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THE  PHYSICS  AND  OPERATION  OF 
ULTRA-SUBMICRON  LENGTH  SEMICONDUCTOR  DEVICES 

ABSTRACT 


This  doounent  simmarizes  activities  under  ONR  Contract:  N00014-86-C-0780, 
under  vdiich  equilibrium  and  nonequilibrium  electron  and  hole  transport  in  micron 
and  subnucron  structures  were  studied  via  a  wide  range  of  numoical  procedures. 
These  included  Monte  Cario  methods,  moments  of  the  Boltzmann  transport 
equation,  Schrodinger’s  equation  and  the  quantum  Liouville  equation  in  the 
coordinate  rq)resentation.  While  all  of  the  studies  have  resulted  in  a  large 
collection  of  publications,  the  basic  thane  of  the  studies  was  the  determination  of 
the  pl^cs  of  de>dce  operation  and  the  influence  of  anall  structure  size  on  this 
operation.  The  most  recent  activities  have  involved  the  quantum  Uouville 
equation  with  emphatis  on  dissipation  and  the  calculation  of  current.  This 
document  includes  a  description  of  quantum  transport  via  the  quantum  Uouwlle 
equation,  as  we  now  understand  h,  as  wdl  as  a  brief  summaiy  of  the  previous 
activities  involving  larger  submicron  devices.  While  the  principle  goal  of  this  study 
was  duddating  the  physics  and  operation  of  nanoscale  dewces,  a  contimiing 
requirement  was  that  all  algorithms  be  menu  driven  and  accessible  to  device 
sdoitists  and  engineers.  The  quantum  transport  algorithm  is  accessible  on  UNIX 
woricstations  and  in  a  PC  V^dows  format. 
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THE  PHYSICS  AND  OPERATION  OF 

ULTRA-SUBMICRON  LENGTH  SEMICONDUCTOR  DEVICES 

1.  PREFACE 

From  it’s  inception  the  study  discussed  bdow,  performed  under  ONR 
Contract  N00014-86-C-0780  has  conconed  hsdf  with  equilibrium  and 
nonequflibrium  electron  and  hole  transport  in  micron  and  submicron  Aructures. 
All  rdevant  equations  and  procedures  were  invoked  and  included  Mome  Carlo 
methods,  moments  of  the  Boltzmarm  transport  equation,  Schrodinger’s  equation 
and  the  quantum  Liouville  equation  in  the  coordinate  representation. 

The  more  classical  problems  emphasized  hot  carrier  phenomena  and 
transients,  while  the  quantum  transport  was  concerned  with  specific  quantum 
phenometui  and  the  best  means  of  studying  it.  Quantum  transport  has  occupied 
most  of  our  actmties  in  the  past  few  years,  and  the  migor  success  in  the  program 
was  the  recent  ability  to  compute  current  self-contistently  within  the  fimnework  of 
a  dissipation  model.  Two  examples  serve  to  illustrate.  This  model  vdien  coupled 
to  earlier  models  now  permits  us  to  deal  vnth  trantients  in  a  sensible  manner  in  that 
the  relaxation  to  an  intermediate  state  is  better  defined. 

The  approach  we  have  taken  is  different  fi'om  those  of  others  because  our 
goals  were  very  general  and  included  the  requirement  that  any  and  all  algorithms 
indude  tools  that  dewce  sdentists  and  en^eers  could  utilize  as  part  of  routine 
device  detign  tasks.  In  otho-  words  one  goal  was  to  include  algorithms  that  would 
be  as  accessible  as  the  standard  drift  and  difiurion  equations. 

The  quantum  transport  equation  we  deal  with  is  the  quantum  Liouville 
equation  in  the  coordinate  representation.  Recall  that  Schrodmger’s  equation  is  a 
coordinate  representation  description.  In  dealing  with  the  quantum  Liouville 
equation  in  the  coordinate  representation  we  broke  new  ground,  particularly  with 
reared  to  donees.  For  example  boundary  conditions  that  workers  typically 
employ  in  solving  the  drift  and  diffusion  equations  were  discarded.  In  it’s  place  it 
was  necessary  to  incorporate  quasi-Fermi  level  conditions  at  the  boundaries  to 
assure  flat  band  contact  conditions.  The  issues  of  Fermi  statistics  was  not  treated 
within  the  fiamework  of  the  differential  equations,  which  would  fomudly  require 
the  introduction  of  the  Dirac  Hamiltonian  into  the  quantum  Uouville  equation. 
Instead  statistics  were  accounted  for  throu^  boundary  conditions. 
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The  calculation  of  current  was  introduced  self-conastently  and  coupled  to 
the  quasi-Fermi  levels.  The  quantum  Liouville  equations  were  also  used  as  a  basis 
for  justi^ng  eailier  and  more  recent  woric  on  the  quantum  potential. 

This  documoit  summarizes  these  studies. 

Mai^  pqjers  were  either  published  or  submitted  for  publication  during  this 
study  and  one  huge  re^new  article  was  initiated.  A  copy  of  each  of  these  is 
induded  v^th  this  report. 
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2.  INTRODUCTION 

Since  the  pioneering  work  of  Tsu  and  Esaki* ,  the  experimental  studies  of 
SoUner  et  al.^  on  double^rarrier  resonant  tunneling  devices,  and  the  superiattice 
detector  woilc  of  Lewne  et  al.^ ,  there  has  been  growing  interest  in  barrier/well 
devices  and  in  the  fundamental  underpinnings  of  quantum  device  operation. 
Further,  following  the  woilc  Datta  et  al.^,  there  has  also  been  rising  interest  in  the 
basic  physics  accompanying  the  Aharonov-Bohm’  effect  in  heterostructures. 
Indeed,  nuyor  advances  in  nuterial  tedinology  has  enabled  device  sdoitists  to 
conjecture  about  new  device  structures  that  both  test  and  illustrate  basic 
fondamental  quantum  physics  issues  of  few  and  mar^  particle  ^sterns.  For 
example  the  issue  of  nonlocality  now  finds  its  way  mto  discussions  of  transport  in 
quantum  devices.  Nonlocality  in  dasacal  phyacs  is  illustrated  by  the  coulomb 
interaction  that  decreases  as  the  square  of  the  distance  between  partides.  In 
quantum  mechanics  there  are  additional  interactions  that  do  not  necessarily  drop 
off  with  distance  and  these  are  discussed  bdow. 

Another  issue  involves  dissipation.  Schrodingar's  equation  as  traditionally 
used  is  dissipationless,  and  if  all  transport  in  subsystems  were  govmied  by 
Sdirodingei's  equation  without  interactions  between  the  subsystems,  all  transport 
would  be  ballistic.  Disapafion  in  quantum  mechanics  is  treated  by  introducing 
additional  systems,  e.g.,  phonons,  and  allo^g  the  additional  system  to  cause  a 
transition  between  states  of  the  origuuil  system. 

A  third  issue,  specific  to  the  treatment  of  dectronic  devices,  is  the 
reservoir.  Traditionally,  the  exanunafion  of  dasacal  devices  involves  the 
qredfication  of  doisities  on  the  bounding  sur&ces,  regarded  as  reservoirs.  Such 
spedfication,  whidi  is  assumed  to  remain  valid  under  bias,  often  involves  the 


'  R.  Tsu  and  L.  Esaki:  "Tiumeling  in  a  Finite  Superiattice,”  Phys.  Lett.,  22, 562  (1973) 
’T.CL.G.  SoUner,  W.D.  Goodhue,  P£.  Tannewald,  CD.  Parker  and  DJ3.  Pedc:  "Resonant 
Tunnding  Through  Quantum  WeUs  at  Frequencies  up  to  2.5  THr;”  Appl.  Pfys.  Letts.,  43,  588 
(1983). 

’BF.  Levine,  KJC.  Oioi,  C.G.  Bethea,  J.  WaUcer  and  RJ.  Malik:  "New  10  micron  Infiared 
Detector  Using  Inteiband  Absorption  in  Resonant  Tunneling  GaAlAs  Superlattioes,"  Appl.  Phys. 
lertK,  50. 1092(1987). 

^S.  Datta,  M.R.  hfeUcdi,  S.  Bandyopadhyi^  and  M.S.  Lundstrom:  AppL  Phys.  Letts.,  48,  487 
(1986). 

’  Y.  Ahaianov  and  D.  Bohm:  Phys.  Rev.,  115, 485  (1959). 
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coiicq)t  of  a  quasi-Fermi  levd,  in  whidi  the  raergy  separation  betwem  the  bottom 
of  the  conduction  band  and  the  Fermi  level  at  the  boundary  remains  unchanged. 

Presently,  our  ability  to  incorporate  these  quantum  mechanical  issues  to 
describe  physical  phenomena  in  ultra  small  devices  and  to  propose  quantum  phase 
based  devices  has  been  evolutionary.  Through  a  coupling  of  experiment,  theory 
ami  numerical  emulation  we  have  been  better  able  to  understand  how  basic 
quantum  mechanical  processes  affect  device  physics.  But  the  'goodness'  of  a 
description  of  quantum  transport  lies  in  the  ability  of  the  theory  to  erqrlain  the 
detailed  experimental  results  obtained  from  such  complex  devices  as,  e.g.,  two 
terminal  resonant  turmeling  diodes  (RTD),  quantum  wdl  superlattice  detectors, 
and  the  more  common  heterostructure  FETs.  However,  the  complerdty  of  the 
RTD  and  the  puzzle  assodated  with  understanding  its  detailed  operational 
principles  has  led  Ferry^  to  describe  it  as  the  fruit  fly  of  quantum  transport  device 
theory.  How  good  is  the  fruit  fly  analog. 

Traditionally,  transport  in  RTDs  and  other  barrier  structures  has  been 
atudyzed  through  implementation  of  the  formula' : 

(1)  y=[fe/(2»)’]J<acv(k)[/,„(£)-A„(£+e,>)]|r(£,^)r 

It  is  the  approximations  assodated  with  this  formulae  that  provide  the  bounds  of 
our  understanding  of  transport  in  quantum  structures.  In  equation  (1)  /j^^is  the 
equilibrium  Fermi-Dirac  distribution  function,  7l(£,(l)is  the  transmission 

coeflSdoit  obtained  fium  solutions  to  the  time  independent  Sdirodinger  equation, 
£is  the  energy  of  the  tunneling  particle  and  #  the  applied  potratial.  As  discussed 
by  Kluksdahl  ^  al.^  a  major  critidsm  of  this  approach  is  that  it  requires  knowledge 
of  the  distribution  function  at  each  tide  of  the  tunneling  interfiu«,  ratho’  that  the 
bulk  like  distribution  far  fixrm  the  tunneling  inter&ce.  Additionally  the  form  of 
equation  (1)  also  implies;  (1)  the  use  of  equilibrium  distribution  functions  to 
describe  a  biased  state,  when  the  biased  resonant  turuiding  diode  is  in  a  non¬ 
equilibrium  state;  (2)  the  nqglect  of  scattering,  although  scattering  would  be 


*D.  K.  Ferry,  Theory  of  Resonant  Tunneling  and  Surface  Superlattices ,  a  diapter  in  Rhymes  of 
Quantum  Electron  Devices ,  F.  C^msso  (ed)  Springer-Veilag,  Berlin  pp77-106  (1989) 

^N.C.  Kltthsdahl,  A.M.  Kriman  and  D.K.  Ferry:  "Self-Consistent  Study  of  the  Resonant 
Tunneling  Diode,*  Phys.  Rev.,  B39. 7720  (1989). 
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required  to  force  a  system  to  a  state  of  equilibrium;  and  (3)  the  concept  of  a  Fami 
levd,  which  clearly  implies  the  presmce  of  strong  carrier-carrier  into'actions, 
particularly  in  the  quantum  well. 

While  the  use  of  equation  (1)  has  been  successful  in  predicting  n^ative 
conductance  in  RTD  its  inadequades  in  explaiiung  experiment  have  been  well 
documented.  These  include:  First:  the  dc  studies  do  not  account  for  the  peak-to- 
valley  ratio  of  resonant  tunneling  devices.  Second,  the  dc  studies  do  not 
adequately  treat  dissipation.  Third,  the  dc  studies  do  not  treat  hysteresis  in  the 
currmt  voltage  characteristics,  observed  experimentally.  Fourth:  the  dc  treatment 
cannot  predict  how  the  dewces  will  be  used  in  applications.  Fifth:  the  dc 
treatment  cannot  treat  the  time  dependent  nature  of  tlw  boundary  conditions  that 
represent  physical  contacts. 

The  above  studies  suffer  from  lack  of  incorporating  the  feature  basic  to 
quantum  mechanical  phenomena:  aU  quantum  mechanical  devices  are  time 
dependent.  Apart  from  dissipation,  there  are  always  reflections  off  boundaries, 
barriers,  wdls,  imperfections  and  contacts.  ffTiai  is  needed  is  a  time  dependent 
large  signal  numerical  studies  of  quantum  feature  size  devices.  This  need  has 
been  discussed  by  Ravaioli  et  al.'  and  Frensley’  and  more  recently  Ferry  and 
Grubin'* .  This  approach  emphasizes  the  details  of  tranrient  behawor,  the  numbm 
of  particles  involved  in  device  operation,  the  temporal  duration  under  which  the 
effective  mass  iqrproximation  is  vafid,  the  Sgnificance  of  the  Fermi-golden  rule, 
and  otho-  short  time  phenomena.  As  currently  practiced,  when  scattoing  is 
presort,  or  when  time  dependent  fields  are  present  and  treated  as  perturbations,  it 
is  supposed  that  the  poturbation  does  not  modify  the  states  of  an  unperturbed 
system,  rather  the  perturbed  system  instead  of  remaining  permanently  in  one  of  the 
unpoturbed  states  is  assumed  to  be  continually  chan^g  from  one  to  another,  i.e., 
undergoing  tranritions  from  one  state  to  another  state.  This  approach  is  at  the 


'U.  Ravaioli,  M.A.  Osman,  W.  Potz,  N.  Kluksdahl  and  D.K.  Feny:  "Investigation  of  Ballistic 
Tianqiort  Through  Resonant-Tunneling  QuaiUum  Wells  Using  Wigner  Function  Approadi," 
/>;b«^ca,134B,36(198S). 

’W.  Fiensley:  "Boundary  Conditions  for  Opmi  Quantum  Systems  Driven  Far  From 
Equilibrium,"  Reviews  ofMotkm  Pl^cs,  62, 74S  (1990). 

"*0.  K.  Ferry  and  R  L.  Gntbin,  ‘Modeling  of  Quantum  Transport  in  Semiconductor  Devices” 
Chap.  In  Solid  State  Pl^cs  (R  Ehrerueich,  ed)  Academic  Press  (1994) 
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heart  of  those  calculations  employing  the  density  matrix*' ,  those  employing  the 
>^gner  distribution  function*^ ,  and  those  emplo^g  Green's  function  techniques*^ . 

In  addition  to  these  fundamental  approaches  there  are  also  derivative 
procedures  that  enjoy  wide  spread  use,  both  for  the  intuitive  nature  of  the 
equations  and  because  of  the  ease  with  ^ch  clasacal  concqits  emei^e.  These 
discussions  include  the  quantum  moment  equations,  see  e.g.,  lafiate  et  al.*^, 
Strosdo*^ ,  and  Grubin  and  Kreskovsky*^ . 

In  the  discussion  that  follows  the  density  matrix  and  quantum  moment 
equations  were  implemented  in  the  study  of  quantum  feature  size  devices.  Further 
we  have  found  inaght  for  multipartide  transport  based  upon  concepts  obtained 
through  a  recasting  of  the  single  particle  Schrodinger  equation.  Adopting  the 
approadi  of  Bohm*^  the  angle  paitide  wave  function  is  written  in  the  form: 

(2)  v=/Jexpj^Y] 

subject  to  the  condition  that  increating  the  phase  by  In,  does  not  change  the 
wave-fimction.  This  wave  function  when  inserted  into  Schrodinger’s  equation 
results  in  two  equations: 


"  See,  e.g..  R  Etuehiekb  and  M.  R  Cohen,  Phys.  Rev.,  115, 786  (1959),  J.  GoMstone  and  K. 
Gottfiied,  II  Nuovo  Qmeato,  13,  849  (1959),  and  nxMC  recently,  W.A.  Frendqr,  Rev.  Mod. 
Phys.,  62,  745,  (1990),  which  include  a  discussion  of  the  density  matrix  in  the  coordinate 
rqiresentation.  Most  recently  a  discussion  by  J.  B.  Krieger  and  G.  J.  lafinate,  Pi^.  Rev.  B35, 
9644  (1987)  and  G.  J.  lafiate  and  J.  B.  Krieger,  Phys.  Rev.  B40, 6144  (1989)  for  a  discussion  d 
the  density  matrix  in  the  momentum  lepiesentatioa. 

'*E.  P,  Wigner,  Phys.  Rev.,  40, 749  (1932) 

”  R.  Lake  and  S.  Datta,  Phys.  Rev.  B45, 6670  (1992) 

'^GJ.  lafiate,  RL.  Grubin  and  DJC.  Feny:  "Utilization  of  Quantum  Disiribinion  Functions  fix 
Ultia-SdbmicnHi  Device  Transport,"  J.  De  Physique,  10,  CJ-3(n  (1981). 

'^M.  A.  Stroscio:  "Kfoment*Equation  Rqnesoitation  of  the  Dissipative  (Quantum  Liouviile 
Equation."  SuperUORcesandhRcrostructures,  2, 83  (1986). 

‘*RL.  Chubin  and  J.P.  Kreskovslty:  "Quaidum  Moment  Balance  Equations  and  Resonant 
Thnneling  Structures,"  Solid  State  Electronics,  32, 1071  (1989). 

'^See,  e.g.,  C.  Philippidis,  D.  Bohm  and  RD.  Kaye:  II  Nuovo  Cimento,  71B,  75  (1982).  More 
recently  see  D.  Bohm  and  B.  J.  Hiley,  The  Unthvided  Universe,  Routledge,  London  (1993) 
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(3)  ^+^+''+6  =  ® 

and 

«> 


where: 


(5) 


Equations  (3H^)  indicate  that  the  Schrodinger  wave  represents  a  partide 
with  a  weU  defined  position  whose  value  is  causally  determined.  The  partide  is 
never  separate  fix)m  the  quantum  force,  -VQ,  that  fundamentally  affects  H.  The 
partide  has  an  equation  of  motion: 

(6)  «^=-V(K+e) 

wludi  means  that  the  forces  acting  on  the  partide  consist  of  the  clasdcal  force, 
-W,  and  the  quantum  force,  -VQ.  It  is  important  to  note  that  the  quantum 
potential  is  dependent  on  the  shape  of  the  real  part  of  the  wave  function  rather 
than  on  its  intensity;  and  does  not  necessarily  &11  off  with  distance.  The  quantum 
force  is  dependent  on  the  momentum  of  the  carrier  through  the  continuity 
equation,  but  does  not  require  a  source  term 

The  quantum  potential  is  defined  in  terms  of  a  angle  particle  wavefiinction. 
And  if  S(r,t)  =  s(r,t)-Et,  where  is  a  constant  independmt  of  poation,  thoi 
undo*  zero  current  conditions,  equation  (3)  is  the  real  part  of  Schrodingers 
equation  whose  solutions  subject  to  a  particular  set  of  conditions  leads  to  a  set  of 
bound  state  eigenvalues.  We  will  come  back  to  this  point  ova*  and  over  again,  in 
the  discusaon  that  follows. 

While  the  above  discussion  is  for  single  partide  wave  functions  we  are 
interested  in  quantal  and  dasacal  distribution  functions,  both  representing  an 
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ensemble  of  particles.  Our  experience,  has  developed  from  iq)proximate 
representations  of  the  Wigner  distribution  function‘d  indicates  that  the 
incorporation  of  the  quantum  potential  for  an  ensemble  of  particles,  where  the 
amplitude  R  is  replaced  by  the  square  root  of  the  self-consistent  density,  p(k),  is  a 
significant  aid  in  interpreting  much  of  the  salient  features  of  quantum  transport 
in  devices.  The  use  of  the  quantum  potential  provides  an  altouative  explanation 
for  the  peaking  of  the  charge  d^isity  at  portions  away  from  the  interfru:e  of  wide 
and  narrowband  gap  structures;  for  real  space  transfer,  for  the  potential 
distribution  associated  with  a  Schottky  barrior,  for  density  variations  associated 
vnth  variations  in  effective  mass,  and  a  host  of  additional  features.  To  get  to  these 
points  we  must  get  through  some  mathematics,  part  of  which  is  exact,  and  part 
approximate.  We  be^  vnth  the  development  of  the  single  particle  density  matrix. 

3.  TfiE  SINGLE  PARTICLE  DENSITY  MATRIX 

While  the  density  matrix  approadi  discussed  below  and  the  Wigner 
approach  are  mathematically  equivalent,  we  have  made  the  choice  of  the  density 
matrix  because  the  equation  of  motion  readily  submits  to  algorithms  developed  by 
the  authors;  the  use  of  wtuch  are  extremdy  short  computational  times  for  steady 
state  solutions.  These  algorithms  are  discussed  bdow.  Tho’e  are  limitations  to 
our  treatment.  The  most  important  is  that  the  equation  of  motion  discussed  bdow 
does  not  include  anti-symmetric  components  and  the  density  matrix  has  not  been 
subject  to  anti-symmetrization‘' .  We  note  that  the  application  of  the  Wigno* 
formulation  to  devices  suffers  from  the  same  limitation.  In  some  of  the  studies 
below,  the  inclusion  of  Fermi  statistics  is  through  the  boundary  conditions,  as  in 
the  Wigner  studies. 

The  structures  that  we  discuss  fidl  under  the  cat^ory  of  open  structured, 
which  can  exchange  partides  with  its  surrounding,  and  which  nu^onatically 
expresses  this  interaction  in  terms  of  boundary  conditions.  The  phenomena  we  are 
interested  in  vrill  be  with  systems  that  are  far  from  equilibrium. 


"a  brief  discussion  d  anti-synun^rization  is  included  in  the  monograph:  "Foundations  d 
Electrodynamics,*  S.R.  De  Groot  and  L.G.  Suttorp,  North-HoUand  Publishing  Company, 
Amsterdam  (1972).  See  also  O^Ckinnell  (get  reference). 
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The  doisity  matrix  is  obtained  from  the  density  operator/?^  (/),  which 
following  Dirac  notation'*  is: 

(7)  Po,(0=Z\^(0>p(i)<m 

i 

vriiere  |i(/)  >  rq)resents  an  dgenstate^ .  The  time  evolution  of  the  doirity 
operator  is  obtained  from  the  time  evolution  operator^'  vduch  has  the 

propertyX/(^^oM^o)^'(0>•  Tlie  time  evolution  operator  is  unitary  and  the 
dependence  of  the  density  operator  on  previous  times  is  given  by: 

(8)  Po,(t)  =  U(tJ,)p^(t,)U(tA)^ 

where  the  symbol  'f  represoits  the  adjoint.  The  time  dependence  of  the  daisHy 
opotdor  is  governed  by  the  time  dqjendence  of  the  time  evolution  operator,  whidi 
is“: 

(9)  )£/(,.,,) 

Assuming  that  the  Hamiltoiuan  H(t)  is  Hermedan,  the  time  dependence  of  the 
denrity  operator  is: 


'^P.  A.  M.  Dirac,  The  Principles  of  Quantum  Mechanics,  Oxford  University  Press,  London 
(1958).  Particular  attention  should  be  paid  to  Section  33,  where  we  note  that  if  P(i)  is  the 
probability  of  the  system  being  in  the /th  state,  it  can  never  be  negative.  If // is  an  eigenvalue  of 
and  I// >  is  an  eigenket  belonging  to  this  eigenvalue,  then  Pg^\pf  ff\ff  >.  As 

discussed  in  sectron  33,  pf  cannot  be  negative. 

"  Note,  later  we  will  be  eiqnessing  our  results  in  the  coordinate  r^resentation.  As  discussed  by 
P.  R.  IfoUand,  The  Quattum  Theory  of  Motion,  Cambridge  University  Press,  Cmnbridge  (1993), 
page  104,  a  mixed  state  may  be  decomposed  in  an  infinite  number  ci  ways,  and  so  we  cannot 
uniquely  deduce  from  it  the  set  of  eigenstates  in  the  ensemble  and  their  reflective  wei^ts.  The 
same  wiU  apply  to  the  Wigner  distribution  fimcdon,  vtadi  is  obtained  finom  it  throu^  a 
transformation,  and  has  the  same  content 
”  Reference  [19],  section  27. 

^See,  e.g.,  A.  h^ssiah  Quantum  Mechanics,  Volume  D,  John  Wiley  &  Sons,  NY,  (1%1), 
particularly  Cliapter  XVn. 
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(10) 


The  density  matrix  in  the  coordinate  represmtation  is  given  by: 


(11) 


^£<^!x:>| 


'  teattirliig 


Notice  that  we  are  ignoring  any  spatial  variation  in  the  effective  mass,  although  we 
will  deal  with  this  later^ .  The  last  term  on  the  right  hand  ade  of  equation  (1 1)  is  a 
generic  representation  of  scattering,  which  we  treat  below  in  a  semiclassical 
maimer.  All  of  the  quantum  features  assodated  with  the  devices  bdow  will  arise 
from  the  streaming  terms. 

The  density  matrix  is  Hermetian,  and  /j(x,x')e<x|Pj^|x'>=P(x',x)*. 
Additional  quantities  rdevant  for  transport  indude  the  current  density  matrix; 

(12)  j(l.x')  =  ^(V.-V,)p(x,.-) 
and  the  energy  density  matrix: 

(13)  £(x.x')  =  -|^(V.-V..)‘/)(x,x') 


The  diagonal  componoits  of  each  represent  the  obsnvables. 
Equation  (1 1)  r^en  coupled  to  Poisson’s  equation : 

(14)  V(xV)l'  =  -*’(p(x)-/>,W) 


^  A  key  study  is  itferenoed  here  G.  T.  EinevoU  and  L.  J.  Sham.  Boundary  conditions  for 
envelope  /motions  at  interfaces  between  dissimilar  materials,  Phys.  Rev.  B49, 10S33  (1994). 
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and  equations  describing  scattering  are  the  rdevant  equations  for  device 
transport.  Note  while  the  above  equations  are  for  electrons,  we  will  also  discuss 
hole  transport;  the  relevant  modifications  to  the  equation  will  be  indicated. 

The  Liou>dlle  equation  in  the  coordinate  representation  is  a  function  of  six 
variables  plus  time.  The  ^  variables  represent  a  coordinate  phase  space  vdiose 
rdation  to  the  standard  phase  space  invohing  portion  and  momenta  may  be 
assessed  through  iq)pIication  of  the  transfonnation^,  which  has  been 
modified  to  include  spin 

(1 5)  p{x,x')  =  j^_^jjdk/^[^k.^ii^jexp[ik .  (X  -  X*)] 

To  date  the  description  of  transport  in  devices  via  the  density  matrix  has 
been  confined  to  cases  where  the  partides  are  fiee  in  two  directions,  >^ch  for 
specifidty  we  take  as  the  y  and  V  directions.  Further  in  the  discussion  bdow  we 
will  deal  with  diagonal  components  along  the  fi’ee  directions,  and  treat  the  density 
matrix  p{x,s^,y  =  y\z^  z')  =  p(x,*')- 

To  determine  the  form  of  the  dendty  matrix,  we  can  picture  a  situation  in 
the  absence  of  disdpation  in  vdiidi  boundary  conditions  permit  the  separation  of 
equation  in  two  Schrodinger  type  equations,  ^th  a  solution  that  is  the  product  of 
two  wave  functions.  More  generally  we  sedc  solutions  of  the  type: 

(i«)  p(x.^ 

i 

for  vdiich  equation  (IS)  is  a  spedal  case.  We  now  condder  several  examples. 

4.  EXAMPLES  OF  THE  EQUIUBRIUM  DENSITY  MATRIX 
For  a  Fomi-Dirac  distribution  function: 

(17)  /.■(k.*)  =  n.exp[(£_£_,)/*,  r] 

and  for  parabolic  bands  the  density  matrix  is: 


”  R  Weyl,  Z  Physik,  46, 1  (1927). 
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,  r  Nc  T  ^  Ir*^  sin[//'^(x-x')/A] 
(IS)  P(.x.r)  =  Jf  J 


Here,  Lim,^p(x,x‘)  =  ti,F,„(fi,),  ^<p{F.)  =  lri3/2)y' 

M  =  (E-Ec)/k,T.  /i.  =(£,-£<•)/*.?■,  Nc  =  r(3/2)/(2;r’X’)  is  the  density 
of  states  and  =h^  l(2mkgT)  is  tfx  square  of  the  thennal  deBroglie 
wavelength. 

There  are  two  limiting  cases  that  submit  to  analytical  expression.  In  the 
high  temperature  limit,  where  Boltzmann  statistics  apply  (the  Boltzmann 
distribution  arises  when  /i,  <  -4) 


(19)  =  K  exp[/i,  -  (x  -  x')*  /4A*] 


This  distribution  is  Gaussian.  For  a  material  sudi  as  gallium  arsenide,  the  thermal 
deBroglie  wavelength  at  room  temperature  is  4.7nm  and  =  4.4xlo23/in3.  For 
a  nominal  density  of  lO^/m^,  s:-1.48.  In  the  low  temperature  limit,  e.g., 

T=OK“: 


(20) 


where  j,(z)  is  a  spherical  Bessel  function.  Eg  =  h^kg  /  2m,  and  kg  =  In 

the  linut  as  z  =>  0,  y,  (z)=^z  /  3.  One  of  the  earliest  applications  involving  equation 
(20)  was  in  a  discussion  by  Bardeen  ^  where  it  was  demonstrated  that  the  electron 
density  profile  a  distance  V  fi-om  an  infinite  barrier  was: 


^See  also  equation  A  5.1.7  in  N.R  Mardi,  SoUtb:  Defective  and  Perfect,  apfxaring  in  The 
Si^le-ParUcle  Den^ty  in  Physics  and  Chemistry,  N.  R  Mardi  and  B.  M.  Ddi,  editors. 
Academic  Press,  London  (1987) 

"  J.  Bardeen,  Phys.  Rev.  49, 653  (1936) 
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Figure  1  (Usplays  the  den^  matrix  corresponding  to  equation  (20)  for  a  density  of 
1023/m3. 


Real  Part  of  the  Density  Matrix,  TbOJ)  K 


Figure  1.  Density  nuxtrix  for  free  particles  weighted  by  a  Fermi  distribution  for 
GaAsaiT=0K.  The  density  is  Kji^ /m\ 

The  osdUation  in  the  density  matrix  along  the  direction  (corrdation 
direction)  normal  to  the  diagonal  is  determined  by  the  argument  of  the  spherical 
Bessd  fiinction.  The  periodidty  dq)ends  on  dentity  as  expressed  by  the  Fermi 
wave  number,  and  suggests  the  possibility  of  a  wavenumber  dqiendoit  resonance. 
The  oscillation  disappears  at  room  t«nperature  where  the  distribution  approadied 
as  Gaussian  as  described  by  equation  (19).  The  progrestive  deoease  in  the 
numbers  of  oscillations  as  the  temperature  increases  is  displayed  in  figure  2,  which 
di^lays  a  cut  of  the  density  matrix  in  a  plane  normal  to  the  diagonal  of  the  doishy 
matrix.  The  effects  of  Fermi  statistics  are  also  more  pronounced  as  the  density  is 
increased  (e.g.,  is  increased)  aiJ  we  expect  this  to  manifest  itsdf  in  the 
osdllatory  character  of  the  denrity  matrix. 


15 


Oemify  Matrix  Along  Cross  Diagonal 


Figure!.  Density  nuOrix  versus  correlation  eBsiance  for  free  particles  wei^ited  by 
a  Fermi  distribution/or  GaAs  €a  T^OK,  77K  and  BOOK..  The  density  is  10^ 

The  denaty  matrix  p{x^x*)  shown  in  figure  1  is  plotted  for  a  raiige  of 
values  of  x  and  x’,  (0<x<200!mit,0<x'<20Qnvif).  The  density  ^observable’ 
p(pe)  m  p(x^x)  is  the  value  of  the  density  matrix  along  the  diagonal  and  b  pbMed 
along  the  physical  coordinate  x.  PictoriaUy,  the  density  is  a  projectiat  of  the 
diagonal  component  of  the  density  matrix  onto  the  x>axis.  The  density  nnatrix 
along  the  cross  diagonal  is  defined  as  e  p{L  -  x,x),  where  L  is  the  length 

of  the  structure;  it  is  shown  as  a  projection  onto  the  x-aris. 

The  above  discusrion  prowdes  an  indication  of  what  the  denrity  matrix 
coordinate  representation  profiles  corresponding  to  standard  classical  equilibrium 
distribution  fimctions  look  like.  It  is  expected  on  phyrical  grounds  that  a  dasdcal 
problon  studied  with  the  clasrical  distribution  fiinction  in  momentum  q)ace  would 
yield  the  same  physical  results  as  that  obtained  with  coordinate  rqrresentation 
density  matrix.  For  example,  classically,  with  the  Boltzmann  distribution,  the 
probability  distribution  is  proportional  to  exp-FCx)/  kfT.  Thus,  ^^dien  a  potential 
energy  change  equal  to  kgTtnlQ  (0.0S9S21  ev  at  room  tenq)erature)  is 
considered,  classical  theory  teaches  that  the  denrity  will  be  reduced  Ity  an  order  of 
magnitude.  Solving  the  equation  of  motion  of  the  density  matrix  for  this  case 
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provides  the  same  result.  If  we  go  to  the  other  extreme  at  T’=OK,  and  recognize 
that  the  Fermi  energy  rdative  to  the  bottom  of  the  conduction  band,  E, -Ec^ 
corre^nding  to  a  denaty  of  10^/m^  is  S4.4  mev,  v«diile  that  corre^ndiiig  to  a 
den^  of  10^  In?  is  1 1.7  mev,  ttoi  introdudng  a  barrier  of  42.7  mev  will  reduce 
the  den^  by  an  order  of  magnitude.  This  is  shown  in  figure  i. 

Ferail  Stetteflcs,  Density  and  Potential  Energy 


Figure  3.  For  GoAs  at  T=0K,  Fermi  statistics,  with  a  step  charge  in  potential 
energy  from  0.0  ev  to  0.0427  ev  (dotted  line),  the  non-self  consistent  spatial 
variation  in  density  (solid  line). 

^)art  fix)m  the  asymptotic  (classical)  values  of  density  &r  fiom  the 
inter&ce  we  point  to  the  local  oscillation  in  denaty  on  dther  side  of  the  interfiice, 
and  make  note  of  the  position  of  the  peak  and  minimum  values  of  den^.  Qasacal 
studies  indicate  that  the  peak  value  of  denaty  occurs  at  the  inta'fiu:e;  while  all 
quantum  mechanical  studies  indicate  thiU  the  peak  is  shifted  away  fiom  the 
interfiu^.  In  a  recent  density  matrix  study”,  devoted  to  Boltzmann  statistics,  it 
was  arudytically  denKxistrated  that  the  den^  could  be  represented  in  eqiulibrium 
as  being  equal  to 


L.  Gnibiii,  T.  R  Goviodan,  J.  P.  Kedcovdgr  and  M.  A.  Strosdo,  Sol  St.  Electron,  36,  1697 
(1993) 
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(22)  f3{x)  =  N,cxp[/i,  -(F(*) +Q(x)/3)] 


In  the  absence  of  the  quantum  potential  the  density  is  detennined  soldy  by 
the  potoitial  en«;gy,  and  so  the  doiaty  for  the  potential  energy  distribution  of 
figure  3  would  be  equal  to  its  left  hand  value  right  up  to  the  potential  barrier,  and  a 
second  (lower)  value  \rithin  the  potential  barrier.  The  finite  value  of  the  quantum 
potoitial  and  its  spatial  variation  is  responsible  for  the  minimum  and  maximum 
values  of  the  density  occurring  away  fi’om  the  interfiux.  This  will  be  discussed  in 
more  detail  below  where  we  also  illustrate  the  value  of  the  quantum  potential. 
The  fiictor  *3*  that  appears  in  equation  (22)  is  discussed  in  detail  in  refisrence  27. 

The  potential  variation  in  figure  3  is  imposed  and  abrupt.  Altemativdy  we 
can  envirion  a  structure  in  which  the  density  dianges  abruptly  at  the  same  point 
(lOOnm).  Then  a  solution  to  the  Liouville  equation  and  Poisson’s  equation  yidd  a 
potential  distribution  whose  values  asymptotically  approach  those  of  figure  3.  The 
potential  distribution  at  the  interfiu:e  is  no  longer  abrupt,  and  the  local  peak  seen  in 
figure  3  is  absent.  Rather,  there  is  a  more  gradual  decrease  in  density  across  the 
interfiu:e,  with  values  that  cannot  be  described  1^  the  classical  distribution,  but 
requite  the  presence  ofthe  quantum  potential.  The  two  dinoensional  density  matrix 
fiar  the  calculations  of  figure  3  are  shown  in  figure  4. 

RmI  Part  of  Um  Owway  Matrix 


Figure  4.  Two  ^menstoncd  density  matrix  from  which  the  results  of  figure  3  are 
obtained. 
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The  ori^  of  the  scales  in  figuie  4  is  closest  to  the  reader  where  the  density 
matrix  has  its  highest  values.  Notice  the  ripples  in  the  density  matrix  closest  to  the 
highest  density  regions.  Ripples  are  also  jmsent  at  the  lower  density  regions  but 
their  period  and  magnitude  are  weaker.  Generally  the  effects  of  Femi  statuttics  are 
more  pronounced  at  hi^ier  denaties,  where  from  equation  (20)  it  is  seen  that  the 
amplitude  of  the  oscillation  increases  and  the  period  decreases,  with  increasing 
d«iaty. 

5.  EQUILIBRIUM  DISTRIBUTIONS  AND  THE  QUANTUM  POTENTIAL 

As  indicated  in  the  earlier  discussion  the  dassical  distribution  fimction 
accounts  incorrectly  for  the  chat^  disbibution  in  the  vicinity  discontinuities  in 
potential  energy  and  cannot  be  used  if  the  goal  is  a  description  of  the  operational 
phyacs  of  devices;  the  quantum  potential  must  be  included.  Additionally,  we  have 
also  used  the  quantum  potential  as  an  aid  in  interpretation.  Several  cases  are 
treated  bdow  vduch  illustrate  the  significance  of  the  quantum  potential.  The 
situation  of  the  resorumt  turmeting  diode  be  treated  separatdy  vdrere  the 
significance  of  the  quantum  potent  is  most  apparent. 

The  first  case  of  interest  is  that  of  a  srig^e  barrier  of  modest  height,  42.7 
mev.  This  value  of  barrier  height  is  the  same  value  as  that  of  the  step  potoitial  of 
figures  (3)  and  (4)  where  the  atymptotic  values  of  density  differed  by  an  order  of 
rruigmtude.  For  the  case  illustrated  in  figure  S,  we  again  consider  a  non  self- 
consistent  calculation,  vnth  a  reference  denaty  of  10^  /m^  T=OK,  Fomi  statistics 
and  a  device  length  of  200nm.  For  the  atuation  whoi  a  vety  wide  42.7  barrier, 
lOOnm  width  and  centrally  placed,  is  considered  it  is  found  that  the  asymptotic 
value  of  density  wHlun  a  central  SO  nm  rp^on  is  equal  to  lO”  /m^  a  result 
expected  firom  the  earlier  discus^rt  There  was  additionally  the  structure  in 
density  at  the  potential  discontinuity  that  was  seen  in  figure  3. 

When  a  narrow  lOnm  wide  barrier  is  considered  the  results  are 
quanthativdy  different.  There  is  a  local  peak  away  fix>m  the  barrier,  but  the 
nuninnun  value  of  density  exceeds  that  associated  vnth  the  wider  barrier.  Of 
interest,  however,  is  the  structure  of  the  quantum  potential,  shown  in  figure  S. 
Rrst  we  note  that  the  magnitudes  of  Q(x)  and  V(x)  are  ^)proximatdy  the  same 
within  the  barrier  region.  The  quantum  potential  is  n^ative  within  the  barrier,  a 
consequoice  of  a  povrive  value  of  curvature  for  the  density  within  the  barrier  (the 
denrity  reaches  a  minirmim  at  x^^lOOnm).  The  quantum  potential  is  poative  in  the 
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r^ons  immediatdy  upstream  and  downstream  of  the  barrier,  vriwre  the  curvature 
of  the  density  is  negative.  The  signs  of  the  quantum  potential  are  consistent  with  a 
density  that  is  below  its  classical  value  immediately  outside  the  barrier,  and  above 
its  classical  value  within  the  barrier  r^on. 


QuMtuiiiaiMl  PolMlial  Enwgy 


OiilMiM(nm) 


Figure  5  Quantum  p<^ntial  (solid  curve)  and  V(x)  (dotted}  for  a  single  barrier 
lOnm  wide. 

The  next  case  of  interest,  which  again  offers  the  quantum  potential  as  a 
vdiicle  for  interpretation  is  the  fiuniliar  sdf-conristoit  charge  distribution 
assodated  wth  a  wide  bandgap/narrow  bandgq>  structure.  Figures  6  through  8 
illustrates  results  using  the  density  matrix  for  a  room  temperature  sdf-consistent 
calculation.  Here  the  device  length  is  200nm  where  for  0<x<100nm,  NIy=l(P/m^ 
for  100nm<r<200nm,  Np^lO^/m^ .  A  barrier  of  300  mev  is  imposed.  While  the 
non  sdf-consistent  calculations  of  figure  3  show  a  reduction  in  charge  dendty 
within  and  near  the  edge  of  the  barrio’,  there  is  nothing  in  figure  3  resembling  the 
extent  of  the  charge  reduction  seat  in  figure  6.  The  contributions  to  this  change 
are  several-fold.  First  the  barrier  of  figure  6  is  an  order  of  magnitude  higher  than 
that  of  figure  3.  Seamtf  the  applied  potential  energy  difference  across  the 
structure  is  chosen  to  ^dd  flat  band  conditions^  arui  thus  equal  to  the  hd^  of  the 
barrier  plus  the  biuh-in  potential.  Third,  the  sdf-consistent  potential  diq)lays 
structure.  What  is  the  origin  of  tlus  structure? 
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0«mlly  aiMl  PolmHal  Enwgy 


Figure  6.  Self-consistent  calculation  of  the  density  and  potential  energy  for  a  300 
mev  heterostructure  diode  ttt  T^SOOK,  with  Fermi  statistics  and  flat  band 
conditions 

In  all  of  the  calculations  with  a  h^ostnicture  barrier,  once  we  pass  the 
peak  density,  there  is  a  progresave  decrease  in  denaty  until  a  minimum  value  of 
density  is  readied  within  the  interior  of  the  heterobarrier.  The  ample  explanation 
based  upon  the  quantum  potential  imhcates,  from  equation  (6)  that  the  net  force, 
under  zero  current  conditions  is  zexo.  But  the  quantum  mechanical  self  force, 
generated  by  variations  in  the  single  particle  density  (from  the  quantum  potential 
as  seen  in  figure  7)  is  alwctys  nonzero.  Here  as  we  move  into  the  wide  band  gap 
r^on  vdiere  the  denaty  is  decreaang  and  qiproaching  a  minimum  value,  the 
curvature  of  the  density  is  positive,  resulting  in  a  n^ative  value  for  the  quantum 
potential  Since  there  is  a  miiumum  value  of  the  density  i^thin  the  wide  band  gq> 
r^^n,  there  is  structure  to  the  quantum  potential  leading  to  a  spatially  dqiendent 
drivirig  force.  This  force  must  be  balanced  by  variations  in  the  sdf  conastent 
potential  as  seen  in  figure  7.  The  self-consistent  potential  vdiich  is  driven  1^ 
Poisson’s  equation  is  now  subject  to  tlw  additional  constraint  imposed  by  the 
quantum  potential  The  details  are  not  govoned  by  equation  (6),  rather  th^  are 
governed  by  the  Liouville  equation;  but  the  qualitative  features  are  rqiresented  by 
equation  (6).  When  examining  the  dasacal  atuation  we  note  that  the  potential 
energy,  is  also  constrained  by  a  diffusive  contribution.  Diffusive  contrilnitions  are 
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also  present  whoi  quantum  transport  is  considered.  The  quantum  potential 
contribution  is  an  additional  contribution  that  is  not  dependoit  upon  the  presence 
of  difiusion. 


PotantW  Energy  and  Quantum  Poianllal 
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Figure  7.  Self-coimstent  calculation  of  the  quantum  potential  and  potent 
energy  for  a  300  mev  heterostructure  diode  at  T-300K,  with  Fermi  statistics  and 
flat  band  amditians 

There  are  sevoal  interesting  additional  points  concerning  the  structure  of 
the  charge  distribution  associated  with  the  calculations  of  figures  (6)  and  (7).  A 
good  iq)proximation  to  the  curvature  of  the  potential  energy  within  the  wide  band 
gap  region  and  near  the  intofitce,  is  to  assume  that  the  rqgion  is  fiue  of  mobile 
carriers,  wherd)9||^^F^(x)  =  (e^ /£)po(^)'  As  a  consequence,  the  higher  the 

heterobarrier,  the  larger  the  vddth  of  the  depletion  zone  on  the  wide  band  gap  side 
of  the  structure.  Under  flat  band  conditions  vdiou  the  net  charge  distribution  is 
zero  there  is  a  corre^nding  increase  in  charge  on  the  narrow  band  gap  side,  and 
this  accumulated  charge  will  increase  vnth  increaring  barrio’  hdght.  Thus  unlike 
the  non-sdf  consistent  calculation  of  figure  (3)  tho-e  is  significant  charge 
accumulation  on  the  iuutow  band  ride  of  the  structure.  The  quantum  potential 
whidi  is  negative  on  the  wide  band  giq)  ride  and  therefore  yields  a  larger  than 
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classical  result  for  the  particle  density,  also  has  the  effect  of  yidding  a  lower  than 
dasacal  result  for  the  d«)sity  ju^  outside  of  the  barrier.  The  small  region  of 
n^ative  quantum  potential  to  the  left  of  the  barrier  is  a  consequoice  of  the 
quantum  potential  defined  in  terms  of  the  square  root  in  density.  An  expansion  of 
the  quantum  potoitial  leads  to  contributions  fi'om  the  square  of  the  first  derivative 
of  dentity  as  wdl  as  the  second  derivative. 

What  is  the  situation  with  multiple  barrier  structures;  the  simplest  being 
the  double  barrier  resonant  tunneling  structure.  The  characteristic  feature  of  the 
multiple  barrier  structures  is  the  existence  of  quati-bound  states  within  between  the 
banian.  The  dentity  between  the  barriers  depends  upon  the  barrier  hdght,  barrier 
configuration,  doping,  etc.  As  discussed  eailier^,  the  value  of  the  quantum 
potoitial  vnthin  the  quantum  well  of  a  double  barrier  structure  is  approximately 
equal  to  the  energy  of  the  lowest  quasi  bound  date,  rdative  to  the  bottom  of  the 
conduction  band.  We  note  that  in  terms  of  the  definition  of  the  quantum  potoitial, 
under  steady  state,  zero  current  conditions,  it  is  direct  to  show  fi-om  Schrodinger’s 
equation  that  Q(x) +V(x)  =  E,  where  £  is  the  oiergy  of  the  quasi>bound  state  (see 
also  reference  20).  We  illustrate  the  quantum  potential  for  a  200nm  structure 
double  barrier  structure.  Th^e  are  two  barriers  5nm  wide,  each  300  mev  high, 
separated  by  5  nm,  placed  in  the  center  of  tiw  structure.  The  background  drying  b 
10^/m^  and  uniform,  except  in  the  interior  40nm  rq^on  where  it  b  reduced  to 
10^/m^.  Figures  8  and  9  show,  respectively  the  dordty  and  donor  distribution,  the 
quantum  potential  and  the  self  contistent  potential  energy. 

With  respect  to  figures  8  and  9,  we  note  that  carriers  in  excess  of 
4xl0^/m^  reside  ^thin  the  quantum  well.  The  quantum  potential  is  n^ative 
within  the  barriers  of  the  structure  corresponding  to  the  curvature  of  the  doitity, 
and  is  positive  vdtlun  the  quantum  well.  But  the  remarkable  feature  is  that  the 
quantum  potential  is  approximately  constant  ^thin  the  quantum  well.  We  have 
found  that  for  the  300mev  barrier,  the  quantum  potential  vidthin  the  well  is 
approximately  84mev  (for  a  200mev  barrier  the  quantum  potential  vdthin  the  wdl 
is  approximately  70mev).  A  key  feature  in  utiliang  the  density  matrix  in  the 
coordinate  representation,  is  that  the  quantum  potential  behaves  like  a  quasi-bound 
state. 


L.  GnAin,  J.  P.  Kieskovsky,  T.  R.  Govindan  and D.  K.  Ferry,  SemJ.  Sci  Tedmohg.  (1994) 
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Density  and  Background 


Figure  8.  Self-consistent  T^SOOK  calculatiat  with  Fermi  statistics  showing  the 
density  and  clonor  distribution /<»•  a  s^nmetric  double  barrier  structure. 


Quantum  and  Poteotfal  Energy 


Figttre  9.  For  figjure  8,  the  quantum  cmd  potential  energy  (Sstributimt 
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Furth«'  evidence  for  use  of  the  quantum  potential  within  the  well  as  a 
measure  of  the  energy  of  the  quasi-bound  state  was  pro\dded  by  supplonental  vnth 
calculations  in  which  the  double  barrier  structure  was  placed  within  a  40nm  wide 
quantum  wdl.  1  he  depth  of  the  quantum  well  was  varied.  As  the  depth  increased 
the  quantum  potential  between  the  barriers  ronained  independent  of  position,  but 
increased  slightly  in  value.  The  situation  when  the  quantum  well  was  ISOmev 
deep,  resulted  in  a  value  of  the  quantum  potential  between  the  barriers  that 
increased  to  94mev.  The  detailed  results  are  difierent  than  that  of  figures  8  and  9 
in  that  the  density  between  the  barriers  has  increased^ . 

6.  DISSIPATION  AND  CALCULATION  OF  CURRENT 

The  calculations  of  the  density  and  potential  profiles  for  the  barrio 
structures  in  both  non  self  consistent  and  seif  conristent  studies  indicate  that  for 
distances  sufiBciently  for  from  the  interfitce  the  results  are  the  same  as  that 
expected  using  the  dissipationless  Boltzmann  (or  Vaslov)  equation.  When  current 
flows,  classical  device  transport  studies  usually  proceed  via  the  drift  and  diffusion 
or  hydrodynamic  equations,  or  through  solutions  to  the  Boltzmann  transport 
equation  and  Monte  Carlo  procedures.  Here,  Ibr  cases  where  the  ends  of  the 
device  are  heavily  doped  N^  re^ons,  boundaiy  conditions  on  the  numerical 
procedures  are  invoked  to  assure  that  the  numbers  of  particles  leaving  and  entering 
the  structure  are  the  same.  An  alternative  approach  that  should  yield  the  same 
results  \rith  respect  to  charge  and  potential  energy  distributions  at  the  boundaries, 
is  to  implement  procedures  recognLdng  that  dissipation  at  the  beginning  and  ends 
of  the  structure  may  be  represented  by  carriers  that  thermalize  to  a  local 
equilibrium.  The  issue  then  is  how  is  to  deal  with  this  situation.  To  date,  very 
iq)proximate  methods  have  been  introduced,  and  a  rational  for  this  ^proach  is 
discussed  below,  but  it  is  emphasized  that  some  procedure  for  dissipatirm  must  be 
invoked  if  tran^rt  in  devices  is  to  be  discussed  sensibfy. 

One  of  the  most  succinct  way  to  express  the  problem  of  dissipation  follows 
that  of  Caldeira  and  Leggett^ .  We  consider  a  system  A  (the  device)  interacting 

^This  inciease  in  density  has  at  least  two  origins:  (i)  the  increased  density  on  either  side  rS  the 
barriers,  and  (ii)  the  lowering  of  the  quasi-bound  state  relative  to  the  Fermi  energy  of  the  entering 
carriers. 

A  O.  Caldeira  and  A  J.  Leggett,  Physica.  121A  587  (1983). 
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with  a  second  system  B  (the  reservoir)  described  by  the  Hamiltonian 
'^^Rmmyeir  mcorporatuig  the  reservoir  and  the  Hamiltonian 

describing  the  interaction  between  the  two  systems.  The  breakup  between  the 
dewce  and  reservoir  is  problem  dependoit.  If  upper  case  letters  denote  the 
coordinates  of  the  reservoir  and  Iowa*  case  letters  the  coordinates  of  the  system  of 

interest  (e.g.,  the  dectron  ^stem)  then  the  quantity  interest  is  the  denaty  matrix 
<  This  quantity  describes  the  behavior  of  the  oitire 

^stem.  We  do  not  need  detailed  information  about  the  reservoir,  rather  we  need 
to  determine  its  influence  on,  in  our  case,  the  electron  system,  which  implies: 

p{x,x\i)  =  J  dlWR*  <  xR|e-""V„^(0)e“""*|x'R'>. 

One  method  that  has  been  invoked  to  deal  with  dissipation  and  boundaries 
and  current  flow  in  devices,  has  been  guided  by  pertuibation  theory  on  the  density 
matrix  .  First  the  equation  of  motion  of  the  denaty  matrix  has  been  rewritten  to 
include  a  scattoing  contribution,  as  shown  by  equation  (11).  Below  we 
concentrate  on  the  modifications  of  the  Liouville  equation  through  the 
incorporation  of  scattering  and  deal  only  with  the  Liouville  equivalent  of  classical 
scattering. 

In  the  Boltzmann  picture,  ignoring  Fermi  statistics,  the  scattering  rate  is: 


(21) 


dk'{/.(k.xMx,k’,k)-/,(k’,i))F(i.k,k’)} 


where  the  subscript  *w’  denotes  a  Wigner  function  and  ir(z,k,k')  represents  the 
standard  transition  probability  per  unit  time.  Utiliang  the  Weyl  transformation: 

(22)  p(r+s,r-s)=^J<fk/,(k,r)exp/2sk 

with  the  following  change  in  coordinates:  x+x’=2r,  x-x'=2s,  the  scattering 
rate  of  density  matrix  (after  manipulation  of  the  variables  of  integration)  is  given 
by: 


’’  R  L.  Gnibin,  T.  R.  Govindan  and  M.  A.  Stroscio,  Semi.  Set.  Technolog.  (1994). 
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(23) 


^^p(r-n.r-i)j  ^ 

-|^^j  jAA*{/>,rXexp/2s  k]»F(r,k\k)Kl-expi2^ 


The  structure  of  the  scattering  tom  within  the  frameworic  of  dasacal 
Boltzmann  scattering  expressed  within  tiw  coordinate  rq)resentation  is  obtained 
from  equation  (23).  For  example  the  second  exponential  term  in  equation  (23)  can 
be  esqjressed  as  an  infinite  series,  in  which  case  the  scattering  term  is  expressed  as 
an  infinite  smes  in  powers  of  s.  The  lead  term  is  ^ven  by: 


(24) 


^^p(r-Hs.r-8)j^^ 


Standard  classical  theoiy^^  teaches  that: 


^;^iere  r(r,|k|)  represents  a  scattering  rate.  Thus: 

(26)  =-/2s|^Jiik/,(k,rXexp/2sk]kr(r,|k|)j 

whidi,  uring  the  inverse  of  the  Weyl  transformation: 

(27)  2/,(k,r)=2’J*/p(r+..r-  s)exp-/2s*k 


can  be  rearranged  as 


”D.  K.  Feny,  Semiconductors,  Macmillan  Pid>lislung  Company,  NY,  (1991). 
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(28) 


-/2sr^J**'/<r+s',r-i')[exp<2(i-i').k]kr(r.|k|)J 


A  significant  sunplification  arises  when  the  ciystal  momentum  in  equation  (28)  is 
replaced  by  a  divergence  of  the  correlation  vector: 


(29) 


<?p(r+s,r-s)  _ 

dt 

Jteat 


-IV.  ■  />(r+s',r-s')[«p/2(f-i')k]r(r,|k|) 


For  the  case  when  the  scattering  rate  is  independent  of  momentum,  the  dissipation 
term  reduces  to: 


(30) 


and  the  UouwUe  equation  in  the  coordinate  represoitadon  is  modified  to  read: 

dt 


The  additional  contribution  due  to  disripation  was  discussed  in  refo'oice 
[27]  and  in  a  study  by  Dekker’^ .  Denrity  matrix  algorithms  incorporating  the 
dissipation  contributions  of  equation  (31)  have  been  implonented  with  some 
results  rq)orted^^  But  because  of  numerical  difficulties  at  higher  bias  levels 
modifications  to  the  scattering  were  introduced  whose  consequences  go  beyond 


“H.Dekker,P^/?ev.  A16, 2126  (1977) 
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the  iq)proxiinations  associated  with  the  expansion  of  equation  (24).  It  is 
worthwhile  dwelling  on  these  nuxlificadons. 

In  modifying  the  scattering  term  in  equation  (31)  it  was  first  recognized 
that  the  dissipation  term  could  be  reexpressed  in  terms  of  a  vdodty  density  matrix; 

P2)  j(i..')  =  [^](V.-V,)p(.,.') 


The  diagonal  elements  of  equation  (32)  yield  the  vdodty  flux  dendty.  In  terms  of 
j(z,z')  equation  (30)  becomes; 


(33) 


dt 


The  scattering  term  in  the  above  equation  was  then  written  in  the  form  of  a 
scatterii^  potential.  The  procedures  for  this  were  as  follow.  First,  the  term 
j(x,z')  was  rewritten  as  j(x,x')Ev(x,x')p(x,x'),  wh^e  v(x)k  v(x,x)  rq>resents 
the  expectation  value  of  the  velodty.  Second,  j(x,x')  was  ^proximated  as 
j(x,x')wv(x)p(x,x').  Higher  order  terms  are  at  least  second  order  in  (x-x*),  and 
retaining  them  would  be  incondstent  vnth  the  approximation  leading  to  equation 
(24).  Third,  quad-Fermi  levels  were  introduced  through  the  definition; 

(34)  £,(.)-£,(i')=-J[;*'W(x")».r(i"), 


For  small  values  of  x-x*  about  x,  equation  (34)  is  approximately  rq)resented  by; 
E^{x)  -  £^(x' )  =  -(x  -  x' )  •  v(x)wr(x).  Undo*  this  approximation  equation  (33) 
becomes; 


(35) 


dt 
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Thus  we  hsve  taken  the  differential  equation  (31)  ^ose  right  hand  side  is  complex 
and  reidaoed  it  by  one  whose  right  hand  side  is  real,  when  the  density  matrix  is 
real.  Side  by  side  calculations  at  low  values  of  l»as  yidd  identical  results. 

While  the  above  discussion  leadii^  to  equation  (35)  q>pears  to  be  modd 
dq)endent,  the  results  implied  by  this  equation  have  greater  generality  than  the 
means  used  to  arrive  at  H. 

The  implementation  of  equation  (35)  permits  us  to  calculate  current  in  a 
direct  manner.  How  is  this  done?  In  all  of  the  calculations  iMth  current  the 
assumption  is  that  the  carriers  at  the  upstream  boundary  are  in  local  equilibrium 
and  that  the  (tistributions  are  dtho*  a  di^laced  Maxwdlian  or  a  disfdaced  Fermi- 
Dirac  distribution.  As  discussed  in  reference  [27]  this  implies  that  at  the  upstream 
boundary,  the  zero  current  quantum  distrilmtion  fimction  p(x,x')  is  rqrlaced  by 
p(^\')tx3^imvibomdoay)'{x-x')IK\.  Since  current  is  introduced  as  a 
boundary  condition  to  the  problem  as  formulated  by  equation  (35)  a  prescription 
is  necessary  for Jincbr^  its  vcdue.  An  auxiliary  condition  was  constructed 

To  compute  a  value  of  current  for  use  in  the  Uouville  equation,  a  criteria 
was  introduced  through  momoott  of  equation  (35)^.  Under  time  indqrendent 
conditions,  the  momentum  balance  equation  yidds  the  condition: 

(36)  2V.£+(V.FMx)-(V.£^Kx)  =  0 

where  E  is  the  kinetic  energy  and  is  given  by  the  equation  (13).  Under  the 
assumptions  of  current  continuity,  i.e.,  /7(z)v(x)  is  independent  of  distance 
(satisfied  for  the  Liouwlle  equation),  and  the  condition  thcd  the  energy  of  the 
entering  and  exiting  carriers  are  equal,  equation  (34)  becomes: 

(37)  £,(1)  -  £,(*•)  = 

vdiere  we  have  restricted  the  considerations  to  one  space  dimension.  7>i%  current 
is  chosen  so  that  Ef{L)-Ep{Q),  is  equal  to  tiw  change  in  applied  energy  oktoss 
the  structure. 


^Thece  moment  equations  are  discussed  in  reference  (27)  and  are  incorporated  into  a  later 
section. 
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We  now  iilustrate  some  of  the  above  considerations.  The  simplest  type  of 
calculation  to  deal  with  is  that  of  a  fiee  particle.  For  this  case  and  with  current 
introduced  as  a  boundary  condition  the  dentity  matrix  is  complex.  The  real  part  is 
symmetric  and  the  imeginary  part  (from  which  current  u  obtained)  is  asymmetric 
about  the  diagonal.  The  calculation  displayed  in  figure  (10)  tiiows  the  real  part 
and  figure  (1 1)  the  imagiiuuy  part  for  a  200nm  with  a  doping  of  10^/m^  subject  to 
a  bias  of  lOmev.  For  this  calculation  and  parameters  appropriate  to  GaAs,  a 
scattering  rate  of  10*'^sec,  yields  a  mobility  of  0.2S8m^Ar-sec.  The  mean  carrier 
velocity  for  this  calculation  is  approximately  1.3xl0Wsec. 


Real  Part  of  tfM  Density  Matrix 


Figure  10.  Real  part  of  the  density  matrix far  a  free  partide  subject  to  a  constant 
force. 

Increating  the  applied  bias  results  in  an  increase  in  the  carrier  velocity  and 
an  increase  in  the  kinetic  energy  of  the  carriers.  This  ino-ease  affects  the  curvature 
of  the  density  matrix  in  the  correlation  direction  and  b  displayed  in  figure  12. 


Comliion  DWmm  (nm) 

Figure  12.  Density  matrix  versus  correlation  instance  when  current  is  flowir^. 
Dashed  line  is  for  a  bias  of  10  nwv  and  a  mean  velodty  cf  1.3x^  m/sec;  s(did 
line  is  for  a  bias  200  mev  and  a  mean  velocity  of  2.6xl&  m/sec. 


All  Mmiooiidiictor  devices  sustain  energy  dq>endefit  scatteriitg,  implying 
that  the  scatterii^  rate  within  one  region  of  the  structure  will  be  different  than  at  a 
diflbrent  fagkm  (ff  the  structure.  To  understand  how  this  is  implemented  in  the 
density  matrix  algorithm  several  illustrative  examples  of  nonuniform  scattering 
were  performed.  These  examples  deal  with  the  generation  of  nomuiiform  fields 
firom  variatkm  in  the  mobility  (w  scattering).  We  will  treat  an  dement  with 
material  parameters  nomitudly  the  same  as  those  assodated  with  figures  10-12. 
However,  here  we  vary  the  scatterir^  rate  within  the  central  Snm  of  the  structure. 
On  the  bads  of  the  definition  of  the  quasi-Fermi  energy,  a  decrease  in  the 
scattering  time,  whidi  results  in  a  decrease  in  mobility,  will  yidd  a  sharp  drop  in 
the  quasi-Fermi  levd.  The  density  cannot  change  as  rapidly,  but  is  constrained  by 
the  Debye  length  and  so  results  in  a  more  gradual  change  in  the  self  consistent 
potential  energy.  The  quasi  Fermi  energy  and  potential  energy  as  well  as  the 
density  are  displayed  in  figure  13  for  a  bias  of  10  mev,  where  the  scattering  time 
within  the  central  Snm  was  10*'^sec,  while  that  at  the  boundaries  are  req)ectively 
10''^sec.  There  are  several  points  to  emphadze.  For  the  calculation  of  figure  13 
the  quad-Fermi  energy  varies  in  an  approximatdy  liiHxur  manner  in  three  squurate 
regions.  In  particular  within  the  exterior  dadding  ipgions  the  quad-Fermi  levd  is 
equal  to  the  potential  enogy  distribution  where  it  assures  tlw  presence  of  local 
charge  neutrality.  The  departure  of  the  potential  energy  fi-om  the  quad-Fermi 
energy  for  this  calculadon  is  in  large  part  a  consequence  of  Dd>ye  length 
condderations.  The  quad-Fermi  energy  ^di  is  an  integral  expresdon  follows  the 
same  dope,  to  the  interior  re^on,  who-e  the  predpitous  change  in  value  is 
consequence  of  the  reduction  in  the  scattering  time. 

Figure  14  displays  the  scattering  rate  used  in  the  calculations  and  the  sdf- 
condstoit  density  distribution.  Of  extreme  dgnificance  here  is  the  formation  of  a 
local  dipole  layer  within  the  interior  of  tlw  structure. 
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Figure  13.  Setf-constOeta  calculations  of  the  potential  energy  and  quasi-Fermi 
energy  for  a  wdform  doped  structure  with  a  variable  scattering  rate  within  the 
center  of  the  structure. 


DMuty  and  Scattaiiiis  Variattam 


I 

I 


Figure  14.  Self-consistent  calculation  of  the  density  for  a  uniform  d(q)ed  structure 
with  the  diq>layed  variable  scattering  rate.. 
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7.  SINGLE  BARRIER  DIODE:  CONSTANT  SCATTERING  RATE 


The  quasi-Fermi  scattering  modd  has  been  appUed  to  a  variety  of 
structures  including  sin^e  and  multiple  barrier  diodes  as  wdl  as  dectron-hole 
tran^KHt.  We  illustrate  single  barrier  calculations  in  figures  IS  through  18  for  a 
structure  voth  a  constant  scattering  rate.  Preliminary  results  for  tlus  type  of 
structure  were  presented  eatiier^^.  The  calculations  are  for  a  200iun  strudure 
containing  a  tingle  300  mev  high,  20iun  unde  barrier  embedded  within  a  30iun  N~ 
r^on,  surrounded  by  uniformly  doped  10^/m^  material.  The  scattering  time  x  is 
constant  and  equal  to  10*^  sec.  The  calculations  are  self-contistent  and  assume 
Fermi-Dirac  boundary  conditions.  The  first  three  figures,  IS  through  17  show 
potential  energy,  density,  and  quasi-Fermi  oiergy  distributions,  respectivdy,  for 
difiSnent  bias  levds. 

From  figure  IS  as  the  collector  boundary  is  made  more  n^ative  with 
reqrect  to  the  emitter,  a  local  'notdi'  potential  wdl  forms  on  the  emitter  tide  of  the 
barrier.  The  potential  energy  decreases  linearly  across  the  barrier,  signifying 
n^^le  charge  within  the  barrier,  followed  by  a  broad  r^on  where  the  potential 
energy  decreases  to  its  value  at  the  collecUH-  boundary. 

Hie  charge  distribution,  figure  16,  di^lays  a  buildup  of  charge  on  the 
emitter  tide  of  the  barrier,  and  a  compensatory  r^on  of  charge  dq>letion  on  the 
collector  tide  of  the  barria*.  At  a  bias  of 400  mev  significant  charge  accumulation 
has  fiirmed  on  the  emitter  tide  of  the  barrier,  foOowed  by  a  broad  n^on  of  charge 
depletion  on  the  collector  tide.  Note  that  as  the  bias  increases  there  is  a 
progrestive  increase  in  charge  witlun  the  interior  of  the  barrier.  Both  remits  are 
consistent  wth  the  low  temperature  exp«imental  findings  of  Eaves  et  al^ . 

Figure  17  displays  the  quati-Fomi  energy  (rdative  to  the  equilibrium  Fermi 
energy).  Becuise  of  the  low  values  of  current  is  approximately  zero  fitrm  the 

emitter  to  within  the  first  half  of  the  barrier’^  and  then  drops  to  a  value 
^)pro]dmatdy  equal  to  the  bias  through  the  remaining  part  of  the  structure. 


’’D.  K.  Ferry  and  K  L.  Giubin,  Procee^ngs  of  the  International  Workdiop  on  Computational 
ElecPvnics,  Univ.  of  Leeds  247,(1993) 

“L.  Eaves.  F.  W.  Sbeaid,  and  G.  A  Toombs,  Phy^cs  of  Quantum  Electron  Devices  (ed.  F. 
Capasso),  107  (1990)  Springer  -Verlag ,  Beriia 

’^In  the  emitter  r^on  the  variation  in  Ef  matches  that  of  VOc).  and  insures  that  p  x  is 
constant  in  the  vicinity  d  the  emitter  boundary. 
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Potantiai  Elwrgy  vMSus  DManc* 
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Figure  15.  Self  consisteni  room  temperature  potential  energy  calculations 
assunung  Fermi-Dirac  boundary  conditions  for  a  single  barrier  structure  urukr 
varying  bias  conditiorvs. 

I  versus  V  for  the  2O11111  barrio'  is  shown  in  figure  18.  Note  that  for  a  broad 
range  of  voltage  the  current  depends  exponoitially  on  voltage;  but  tlwre  is  distinct 
sublinearity  to  the  curve.  In  words,  the  sublinearity  indicates  that  at  a  ^ven  value 
of  voltage  the  current  is  lower  than  expected  on  the  baris  of  a  pure  esqwnential 
rdation.  In  seeking  an  origin  of  tins  sublinearity  we  note  firom  the  aocompaitying 
voltage  distributions  that  not  all  of  the  voltage  fidls  across  the  tunnd  barrier, 
indeed  a  substantial  contribution  fidls  across  the  re^on  inunediatdy  adjacent  to  the 
collector  side  of  the  barrier. 
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DMMlIy  varaiM  Distanc* 


Figure  16.  Self  consistent  density  calculations far figure  15. 


Quasl<F«nni  Enaioy  versus  Distance 


so  100 


200 


Figure  17  Quasi-Fermi  energy  distribution  for  the  calculations  of figure  15. 
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Fig^e  18.  Current-voltage  relation  for  the  calculations  of figures  IS  throng  17. 
8.  RESONANT  TUNNEL  DIODE;  VARIABLE  SCATTERING  RATE. 

To  illustrate  the  calculation  for  resonant  tunneling  structures  we  treat  a 
200nm  structure,  wth  two  S  nm  -  300  mev  barriers  sqiarated  by  a  Snm  weU.  The 
structure  has  a  nominal  doping  of  10^  In?  except  for  a  coitral  SOnm  wide  r^on 
where  the  doping  is  reduced  to  10^ /m^.  The  effective  mass  is  constant  and  equal 
to  that  of  GaAs  (0.067mo  );  Fermi  statistics  are  imposed;  the  ambient  is  77K;  and 
current  is  imposed  through  the  denaty  matrix  equivalent  of  a  di^laced  distribution 
at  the  boundaries.  In  these  computations  only  one  set  of  scattering  rates  was  used, 
although  scattering  was  increased  in  the  vicinity  cf  the  double  barriers. 

The  signature  of  the  RTD  is  it*s  current-voltage  relation  with  the  r^on  of 
negative  dififerential  conductivity,  for  the  structure  conndered  this  is  di^layed  in 
figure  19.  The  curroit  is  num^ically  negligible  until  a  bias  of  approximately  SO 
mev,  with  the  peak  current  occurring  at  260  mev,  followed  by  a  sharp  but  modest 
drop  in  current  at  270  mev.  The  intopretation  of  these  results  is  assisted  by 
figures  (20)  and  (21)  and  the  Bohm  quantum  potential.  As  indicated  earlier,  we 
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have  found,  through  an  extensive  number  of  numerical  Emulations,  that  the  value 
of  yi(icJ+Q(xJ,  between  the  barriers  of  an  RTD  is  a  measure  of  the  poEtion  of  the 
quaE-bound  state. 


Cunnant  vs  Applied  Potential  Energy  (77K) 


Bhs(imV) 


Figure  19.  Current  versus  (nuxgnitude)  voltage  for  the  resonant  turmeUng 
structure. 


Consider  figure  20  which  displays  the  equilibrium  self-conEEent  potential 
for  the  RTD.  Also  shown  is  the  value  of  the  equilibrium  Fermi  enogy 
(approximately  54  mev)  and  the  vEues,  at  five  different  values  of  applied  potoitial 
energy,  of  V(x)+Q(x)  within  the  quantum  wdl.  At  100  mev  the  quaE-bound  state 
is  approEmately  equal  to  the  equilibrium  Fermi  energy  and  Egnificant  currrat 
begins  to  flow.  The  current  continues  to  increase  until  the  bias  equEs  260  mev, 
where  there  is  a  sudden  drop  in  current. 
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■quMMyiiiPalMlMSMtgy,  Diy ■ndiat  VH«HQW 
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Figure  20.  Equilibrium  potential  energy  and  the  bias  depencknce  of  V(x)+Q(x) 
within  the  quantum  well  Legend  denotes  collector  bias. 


ttn  SunmmiiV  Biilar  BhiIv 
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Figure  21.  Blow  up  of  figure  20. 


To  see  what  is  happening  we  blow  up  the  re^on  on  dther  side  of  the 
emitter  barrier,  >^ere  we  display  values  of  V(x)+Q(x)  before  the  emitter  barrier 
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and  within  the  quantum  wdl  (figure  21).  Within  the  quantum  wdl  we  see  the 
quaa  bound  state  decreaang  as  the  bias  on  the  collector  is  increaang.  In  the 
r^on  prior  to  the  onitter  barrier  vdiere  a  *notch’  potential  forms  agnifying  charge 
accumulation,  we  see  the  formation  vtith  increased  bias  of  a  r^on  where 
V(x)+Q(x)  is  rdatively  flat.  Of  significance  here  is  that  for  values  of  bias 
associated  with  the  initial  current  increase  the  value  of  V(x)+Q(x)  vtithin  the 
quantum  well  is  greater  then  its  value  in  the  enutter  rqgion.  The  current  readies  a 
maximum  at  the  crossKiver  where  V(x)+Q(x)  in  the  emitter  region  and  in  the 
quantum  well  are  iqiproximately  equal.  (Impl«nentation  of  an  earlier  algorithm, 
generally  resulted  in  solutions  osdllating  b^een  high  and  low  values  of  curroit 
when  this  condition  was  reached).  While  it  is  tonpting  to  associate  V(x)+Q(x) 
within  the  emitter  re^on  with  a  quasi-bound  state,  this  may  be  premature. 

The  distribution  of  potential  energy  V(x)  as  a  fiinction  of  bias  is  displayed  in 
figure  22,  vdiere  the  notch  potential  is  deepened  with  increasing  bias,  signifying 
increased  diaige  accumulation.  This  is  accompaiued  by  a  smaller  share  of  the 
potential  drop  across  the  enutttf  barrier,  rdative  to  the  collector  barrier  r^on. 
Comparing  the  slopes  of  the  voltage  drop  ao-oss  the  emitter  and  collector  barriers, 
we  see  larger  fi-actions  of  potential  energy  fiill  across  the  collector  barrier. 


Bias  Dependent  V(x) 


OWwio*  (nm) 


Figure  22.  Distribution  of potential  energy  as  a  function  of  appUed  bias. 
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Explicit  m  this  calculation  is  disapation  whidi  is  incorporated  through  the 
quasi-Fermi  levd.  Within  the  inanity  of  the  boundaries  the  quaa-Fermi  levd  is 
parallel  to  the  conduction  band  edge.  Indeed^  for  this  calculation  the  quasi-Fermi 
levd  departs  from  the  conduction  band  edge  only  within  the  vidnity  of  the  barriers. 
The  quad-Fermi  levd  is  displayed  in  figure  23  at  a  bias  of 260  mev,  vdiere  we  see 
that  the  quasi-Fermi  levd  is  rdativdy  flat  until  the  middle  of  tlw  first  barrier  at 
which  point  there  is  a  small  drop  in  value  followed  by  a  flat  re^on  vothin  the 
quantum  well.  There  is  a  strong  drop  of  the  quad  Fermi  levd  within  the  second 
barrier. 

The  charge  distribution  accompanying  these  variations  in  bias  shows 
accumulation  on  the  emitter  dde  of  the  barrier  along  with  charge  accuimilation 
vdthin  the  quantum  well.  The  increase  in  charge  within  the  quantum  wdl  and 
adjacent  to  the  emitter  r^on  is  accompai^g  by  charge  depletion  downstream  of 
the  second  barrier,  vnth  the  result  that  the  net  charge  distribution  throughout  the 
structure  is  zero. 

Variations  in  the  quad  Fermi  levd  were  accompanied  by  variations  in 
dendty  and  current  which  were  all  obtained  in  a  self-condstent  manner. 
Supplemoital  computations  were  p«formed  in  which  the  quad-Fomi  levd  was 
varied  by  altering  the  scattering  rates.  Hie  calculations  were  applied  to  the  post 
threshold  case  with  values  for  the  scattering  rate  chosen  so  to  provide  a  large  drop 
in  current.  Indeed  a  current  drop  by  greater  than  a  frctor  of  three  was  obtained 
followed  by  a  shallow  current  increase  with  increadng  bias.  The  dgnificant 
difference  leading  to  these  changes  was  the  manner  in  which  the  quad-Fermi  levd 
changed.  Rather  than  the  shallow  change  depicted  in  figure  19,  there  was  a  larger 
change  in  the  quad-Fermi  level  across  the  first  barrier,  a  result  dmilar  to  that 
obtained  for  dngle  barriers. 
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Potential  and  Quasi  Faraii  Enargiaa  at  260maV 


OMMiM(nnO 

Figure  23  Potential  and  quasi-Fermi  energy  at  a  bias  of 260  mev. 

The  calculations  obtained  for  figures  19  through  23  were  obtained  fiom  a 
new  solution  algorithm  that  was  constructed  for  the  quantum  liouville  equation 
that  permits  a  more  convenient  spedfication  of  boundary  conditions,  in  particular 
\^ien  the  device  is  under  bias.  The  algorithm  is  based  on  a  reformulation  of  the 
govermng  equations  in  wludi  a  higher  order  differential  equation  in  the  local 
direction  [(x+x’)/2]  is  constructed  fi’om  the  quantum  Liouville  equation.  The 
reformulated  equation  behaves  like  an  elliptical  equation  in  the  local  direction 
rather  than  the  hyperbolic  behavior  of  the  quantum  Liouville  equation.  With 
appropriate  boundary  conditions,  solutions  to  the  two  forms  of  the  quantum 
Liouville  equations  are  equivalent.  However  the  refisrmulated  equation  allows 
construction  of  a  more  robust  algorithm  that  provides  detired  solution  bdiawor  at 
the  contacts  by  boundary  condition  spedfication  at  both  contacts. 
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9  THE  QUANTUM  HYDRODYNAMIC  EQUATIONS 


A  detailed  description  of  the  denaty  matrix  under  equilibrium  and 
nonequilibrium  conditions  was  given  in  the  prewous  sections.  Much  of  the  work 
reported  there  was  a  consequence  of  considerable  effort  at  understanding  the 
nature  of  the  quantum  Uouville  equation  in  the  coordinate  representation.  As  such 
many  of  these  results  were  obtained  at  the  end  of  the  program;  with  some  results 
particulariy  the  RTD  results  obtained  afto-  the  program  was  completed. 

Simultaneous  with  this  effort  was  attempts  at  determining  an  understanding 
of  the  quantum  hydrodynamic  equations,  as  it  was  felt  that  these  equations  being 
approximate  in  nature  would  find  greato-  acceptance  in  the  engineering  community 
as  a  vdiicle  for  the  deagn  of  multidimenrional  devices.  A  discussion  of  these 
equations  is  not  given  here  as  an  ectensive  piq)er  along  with  several  snudler 
supplement  al  studies  have  other  been  published  or  will  be  pubUshed.  These 
papers  form  part  of  this  final  report  and  constitute  part  of  this  section. 
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Uses  of  the  quantum  potential  in  modelling 
hot-carrier  semiconductor  devices 
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Ortraci  TlKeagliaMMaa(anaiplaa.vario«tfonnao(tlMqiianliimpolantlal 
ara  aMatinad  ter  aiadaMaa  and  InterpraOng  Mia  oparaMoo  of  aamteondiietor 
davteaa. 


1.  Inlrodacllaa 

Quantum  cOaeU  occur  in  devka  atnidufca  when  tlie 
laicnl  confiacoMnt  dimeiMioBi  (the  diitance  over  urkich 
atniScant  dhantm  in  deniiiy  occur)  are  oompacalik  to 
the  Ihennal  de  Broglie  warelength.  B  b  poaiible  to 
model  quntam  atnictnrei  teilli  a  <|oaiitom  deuitjr  bom 
the  Schiddiager  equation,  the  LioaeSo  equation  or  the 
Wgnereqnatioa  Howerer.modellmg  of  there  atructuree 
haa  alao  evolved  bom  the  ure  of  a  aet  of  daaaicnl 
hydrodriiamic  equatioiia,  which  include  oorreetioiia  in 
the  form  of  the  quantum  potentiaL  The  moat  bmffiar 
forma  uaed  involve  the  Bohm  potential  [1] 


and  the  Wigner  potential  [2] 


correction  to  JE(x)  The  term  ^mrepreaeata  an  eqnililitiom 
energy  term  to  which  the  ayatem  rdaxca  [3].  To  place 
equationa  (1)  in  a  bmiliar  form  requirca  additional 
manipulation,  and  one  common  form  te  (4) 


^  ^  p9(B^)f3) 

tti  \m/  ttx  dx  ttx 


+  ff5:«  +  55i-o 

dx  T 


Pm) 


>0 


dx 

tt«aj3)+yM) 

dx 


-  r w .  0  (2*) 

\  6m  J\  dx*  J  dx  V 
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For  example,  uaing  the  denaky  matrix  in  the  ooordiBate 
lepreaentatioa,  la  which  dhu^tioa  b  modelled  hy  a 
Fokiccr-Plaiiek  oontribotioii,  the  hydrodynamic  equa- 
tiooa  are  of  the  form  P] 


di  VW  dx 

«55i  ^  ^  ^  _  0 

dt  dx  dx  e 


(la) 


(Ih) 


(Ic) 

di  \2m*/  dx  \mjdx  v  m 


wherep^  b  the  mean  momentum,  £(x)b  thetnean  kinetic 
CMTgy  and  repreaeata  the  energy  Buz  £33  (which  in 
the  cbaaical  care  repreaente  the  tranapoct  of  energy)  The 
Wigner  potenUal  b  uaually  interpreted  aa  a  quantum 


whew  2  b  the  thermal  de  Broglie  wavelength  and  E,  b 
aa  equflihrinm  energy.  In  the  above  equaSona.  the 
quantum  potential  haa  been  reduced  by  a  breor  of  three, 
allfaough  thb  redoctioB  b  a  aubjeefof  aome  debatn  We 
have  taken  thb  view,  with  othcra  £3),  that  thb  b  an 
at^natable  paraaMter,  calcobtibna  beM  iBuatrate  the 
eflect  of  varying  thb  parameter. 

It  b  poaaible  to  determine  the  validity  of  the  above 
equatiom  for  quantum  devicea  by  pciibnning  oompomble 
cakubtiona  with  one  of  the  foU  quantum  tranapoct 
equationa.  In  adtBtion,  it  b  uaefoi  to  aak  whether  the 
quantum  potential  aide  in  the  interpretation  of  reaulta. 
The  full  quantum  treatment  uaed  here  iavolvea  the  denaky 
matrix  (in  the  coordmate  repreaentatida)  £6)  The  httcr 
calcubtiooa  indude  the  appropriate  Fen^  or  Boltzmann 
atatbtica,  aithoogh  the  moment  equationa  diacuiaed 
below  involve  only  Boltzmann  atatbticx.  When  a  ootn- 
parboo  between  the  quantum  hydrodynamic  equationa 
and  the  denaky  matrix  b  made,  Boltzmann  atatbtica  b 
aaaumed. 


H  L  Orubln  M  tt 


2.  HMaiMM  imuMl  diod*  widf  >qiilMbrluiii 


3.  Moom  calciiteUoii 


Tbe  ISn(  example  ii  dial  of  a  OaAVAKJaAt  double* 
barrier  (5  nm  bankn  and  S  urn  umI)  ceaooaol  tuouefliof 
atrudure,  with  a  aomiaal  doping  of  10**cm~*.  Tbit 
ttrudure  hat  a  quatMtomd  date  of  UomV.  U  it 
anikipeted  that  withia  the  quaotum  well  would 
determioe  the  rptatMiouod  date  Flgare  1  dwpltyt  a 
blow-up  of  the  regloa  aunouiiding  the  douUe^tricr 
dructortL  In  figure  1(a),  tbe  aeV-oootident  potential 
energy,  K(x),  the  Bohm  a^  Wigner  gnantum  potentiab 
are  compued  It  it  aeen  that  the  Bohm  potential  it 
ooodant  within  the  wefl  at  the  quati-hooad  date  talue. 
There  it  a  amaD  accumnltUon  of  catriett  whoee  trahie 
peaka  in  the  centre  of  the  wefl.  Figaro  l(fi)  tcpiaoet  F(x) 
with  the  mean  bindic  energy  per  particle  for  thit 
compariaon.  It  it  aeen  that  the  bed  repnaentatioa  of  thb 
energy  it  with  the  Wigner  potentiaL  The  doeenem  of  the 
retulla  auggeda  thd  the  Wigner  potential  it  hapoctant 
for  the  atructure  of  the  energy  wUiin  the  wdL 


OMean  ImI 


Ftoure  1.  (e)  Bteer-ep  d  toudbiluai  potedld  energy 
(biotMO  eunfo).  Betim  pettdlat  (fdt  eurue)  and  WItgnar 
petenttal  (doMad  eurue)  wtlMn  and  aurroundtng  a  loaenad 
lunnitting  drueture.  (6)  Blew  up  at  egutlfiiiluni  energy  per 
pertMo  (brotcen  curue).  Betun  poleniM  (tdl  curve)  and 
Wigner  petendet  (dellad  curvd  wIIMn  and  auteunding  a 
foaenenl  tumietltiHi  structure. 


The  next  example  it  a  oompamiive  cakulation  along  a 
line  perpendicular  to  the  oonduetion  channel  of  a  mmnvt. 
For  thit  device  (figure  2)  the  OaAt  region  deft-hand  aide 
of  the  ttrudure)  it  100  am  long  (doped  to  10‘*cm~*) 
and  it  adjacent  to  a  wide-bandgap  (300  meV)  region 
doped  to  10‘*  cm~*.  The  fird  10  nm  of  the  wide-bandgap 
region  it  undoped.  The  hydrodynamic  equatfont  were 
tolved  uting  tte  Bohm  potential  with  three  dilcreat 
vahtet  for  tbe  multipHeative  eondantt:  3  (at  in  equalioot 
(2)).  land  9.  Figure  2(a)  thowt  the  denaity  obtained  liom 
the  denaily  matrix,  while  figure  2(6)  thowt  deadly 
obtained  from  the  hydrodynamic  equatioot.  Several 
poinit  are  worth  emphadzing.  The  dmeture  of  the 
dentity  it  the  tame  for  the  three  quantum  hydrodynamic 
equatioiit  and  for  the  dcadty  matrix  calculation.  The 
value  of  the  dendty  for  the  hydrodynamic  equatioot  it 
doted  to  that  of  the  dendty  nulrix  idMn  the  eondant 
it  between  1  and  3.  It  dtouM  be  mentioned  that  while 
the  value  of  the  nmltipficatiie  eondant  wat  varied  by 
almod  an  order  of  magnitude  in  there  calcuhthmt,  the 
quantum  potential  abo  underwent  chaagea,  and  wat 
dflerent  far  each  of  the  three  hydrodynamic  equation 


to  -«a  t  tt  to 
OitlUH*  iMii 


Dnlwtcc  iMit 

Ftpsire  X.  (t)  Oendty  watrh  celeutdien  et  the  denaily 
wIMn  e  MOorcT  eonaguratlew.  (fi)  Queanmi  aydrodynamte 
eetedatlew  et  de  demtty  tor  three  dWtreni  vatuee  et  de 
■NilllpIlMtive  laelor;  o/b  (brekan  curve),  (tub  curve) 
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Hem$t  HM*  ikm  ilittifr«fr«^  kaat  been  outftM  by 
m  Utnjtt  yriMtr.  At  retbt.  m  wiB  vUyM  tknyyh  m 
Hmmyyt  Laterc^nf,  ^dlmg  •  for  tuytrior  Imago. 


Um«  of  ttw  quaniwii  potootlal 


OwsoK*  foul 


FHforo  b,  (•)  Doping  p«n  eunw)  and  bairior  oonllguraUoo  (brehon 
oum)  ior  m  giwnlMni  hyd>wdyn«nlo  colculoffon.  (A)  EgulMNlum 
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4.  Qoantam  waR  aat  cglMMMW 

The  floal  eiaaiple  of  the  «w  of  (he  g«a(am  pofentiai 
in  wodeBing  devioee  ie  Ihat  of  a  gnaalwn  nr«l  nr 
eakulaUon.  Foe  thitcalCTili(ion.«  one  iBaenitonal  eiioe 
Ihronth  the  aowoe  legioa  ahoiwing  (he  ptanar  doped 
fegioa  in  (he  wide^nadgep  aurteiial  ie  ahowa  in  Bpire 
3(d)^  la  the  twO'dioMasional  atmcture.  the  gafe  fHa  20  na 
info  the  ttniciore,  which  ie  ofhctwbe  nooiaally  doped. 
There  is  a  large  ooooeatiation  of  canwa  adjaoeaf  to  the 
wide-ban^p  region,  doe  (o  the  pbocment  of  the 
delu-doped  region.  The  potential  energy  diqday*  a  liocai 
variation  with  charge  depletion  in  the  wido^mndgap 
region  (figure  3(h)X  AM  of  the  density  within  the  quantum 


w«l  is  a  coaseqoanoeof  the  plaaar>doped  ngionand  the 
harrier  is  depleted  of  oarriecs.  TUs  hlM  csMatlen  shoos 
the  aignificant  advantages  Coond  in  using  the  quantuin 
potential  to  nwdei  senloondoctor  nr  stmotans. 


S.  Conduslena 

la  the  above  discuMion.  we  have  ilhistiated  the  nse  of 
the  quantum  potential  as  both  an  adSonct  for  interpse* 
tatioo  and  as  a  tool  for  examiniag  transport  with  Rm 
hydrodynamic  equatioiis.  While  calculations  using  (he 
quantum  potential  have  been  perfotmed  [7],  and 
iMusttale  (be  significance  of  its  contribution  to  (he  density 
distribution  within  email  devices,  significant  advances 
will  occur  when  one  can  model  two-dimearional  Bow  in 
MODTET  structures.  Preiiminaty  studia  indicate  (hat 
realistic  simulations  can  be  performed  [8]. 
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Ahsiract— This  paper  (i)  examinet  through  numerical  aohiliou  of  the  coupled  eoordbuM  repreteauaioa 
liouviUe  and  Poi^n  equationa,  the  use  of  the  Bohm  quantum  potential  to  repeeient  the  equilibrium 
distribution  of  density  and  energy  in  quantum  feature  sine  structures;  00  tlwcawea  the  development  of 
the  nonequilibrium  quantum  hydrodynamic  (QHD)  equations  with  diasipaliott  through  the  truncation  of 
the  quantum  dutribution  function:  and  Oti)  compares  sdect  results  of  the  QHD  ^uatioiie  incotporating 
the  Bohm  potential  to  the  exact  LiouviUe  equatkm  eolutioot.  The  broad  oondueion  of  die  study  is  that 
for  structures  of  current  interest  such  as  HEMTs,  only  quantum  murhanical  aohitiont,  or  the  incorpor¬ 
ation  of  the  quantum  potential  as  a  modification  of  the  daseical  equations  wiU  permit  representative 
solutions  of  siKh  critical  features  as  the  ebeet  dtarpe  density. 


INTItODUCnON 

Advances  in  crystal  growth  and  processing  tedi- 
niques  have  assured  the  constniction  of  nanoscak 
devices  with  sharp  interfaces.  Concomitantly,  new 
device  concepts  have  emerged,  including  resonant 
tiimuding  structures,  quantum  wires,  quantum  dots; 
and  varknts  of  classical  structures  with  quantum 
features,  e.g.  HEMTs  and  HBTs.  While  all  devices 
are  governed  by  quantum  mechanics,  many  devices 
induding  HEMTs  and  HBTs  do  not  requite  quantum 
transport  for  a  description  of  their  bask  operation. 
Nevotheiess,  quantum  mechanics  is  required  to 
provide  key  electrical  features.  For  example,  HBTs 
sustain  low  levete  of  current  at  low  bias  levels; 
these  currents  ate  dominantly  turmeUng  currents. 
Thermionk  contributions  to  current  occur  at  high 
bias  levels.  Recently,  device  formulations  utilizing  the 
drift  and  diffusion  equations  and  the  moments  of  the 
Boltznuum  transport  equation  were  generalized  to 
inchide  a  description  of  tunneling  currents  (Anootu 
and  lafiatefl],  Grubin  and  Kreskovsky(2]).  These 
newer  studies  indicated  that  quantum  contributions 
ct  the  type  first  consideted  by  Wignet(3],  could 
be  incorporated  as  modifications  to  the  more  tra¬ 
ditional  approaches  to  studying  transport  of  catrien 
through  devices.  Such  an  apiRoacfa  was  taken  by 
Zhou  and  Fetry(4,S]  in  a  study  of  quantum  contri¬ 
butions  to  tranqxirt  in  MESFETs.  How  wdl  do 
the  quantum  mtxlifications  of  dasskal  transport 
represent  actual  transport?  This  question  is  addmsed 
for  a  fimited  number  of  cases  through  comparison 
of  (i)  quantum  “corrected”  solutions  with  (ii)  exact 
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coordauue  representation  solutions  to  the  quantum 
Uouvilk  equation  for  the  density  operator  p^, 
whose  time  dependence  is  governed  by  the  Hamil¬ 
tonian  H: 

aapjdt~\H,p^l  (1) 

The  relevant  quantity  in  the  LiouviUe  simulations  is 
the  density  matrix  p(x,  x',  t)  v  <x|p,,|x'>  whose  role 
is  similar  to  that  of  the  distribution  ftmctitm  in 
dasskal  physks. 

The  proo^ure  for  assessing  the  quantum  contri¬ 
butions  has  two  parts:  First,  approximate  and  exact 
equiliMum  solutions  to  the  dis^rationless  quantum 
Liouville  equation  for  a  variety  tff  structures,  indud¬ 
ing  a  barrier,  are  compared.  The  approximate  sol¬ 
utions  wfakh  arise  from  a  new  procedure,  with  results 
simiiar  to  that  of  Wignei(3],  are  also  expressed  in 
terms  of  the  Bohm  quantum  potential(fi]: 

G,s-(fi»/2m)td*(p)''*/dxT/(p)'«,  (2) 

't^iose  physical  significance  is  addressed.  Second,  the 
quantum  Liouville  equation  with  Fokker-Pland; 
dissipation  mechanisms  is  introducbdl?];  from  which 
a  new  derivation  of  the  quantum  hydrodynamk 
(QHD)  equations  are  obtain^  NoneqvdUbtium  zero 
current  QHD  and  exact  LiouviUe  solutions  are  com¬ 
pared  for  a  simple  heterostructure  diode  configur¬ 
ation  relevant  to  HEMT  structures.  We  confirm  that 
the  simplest  version  of  the  QHD  equations,  the  drift 
and  diffusion  currmt  density  equation,  and  its  zero 
current  solution  are  modified  as  follows(l,2]; 

/(x,r)  =  pp*,re[(F-l-aG,)/k,r-|-ln(p)l/djr  (3) 
p  -  p,exp - IF(x) -1- aQ,(xm,r,  (4) 
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where  <i  is  a  consunt,  determined  aiuilyticaUy  below 
and  in  [I]  to  be  n  »  1/3.  More  often  a  is  cbosoi  to 
provide  the  best  fit  to  exact  results,  and  is  thus 
determined  from  numerical  simulations  as  discussed 
below.  Any  value  of  a  other  than  n  »  1  is  of  concern, 
in  that  arguments  associated  with  the  single  partide 
Schrodinger  equation,  suggest  a  value  of  unity,  see 
e.g.  [2].  Nevertheless,  we  show  for  coitditions  appro* 
priate  to  Boltzmann  statistics,  that  the  exact  and 
approximate  solutions  for  n  ft  l  are  remarkably  simi* 
lar;  and  that  solutions  without  quantum  contri¬ 
butions  will  not  represent  the  local  diarge 
distribution  in  barrier  dominated  structures.  Finally, 
the  results  are  related  to  earlier  work  on  the  Wigner 
function  for  mixed  and  pure  states(8].  These  latter 
issues  are  addressed  in  ^  appendices,  which  also 
include  a  discussion  of  the  numerical  algorithm. 

THE  EXACT  EQUATION  OF  MOTION  FOB  THE 
DENSrrV  MATBIX 

The  Liouville  equation  in  the  coordinate  represen¬ 
tation  without  dissipation  is: 

ap/di  +  (k/2mi)(v5-vi:.)p 

-  (1/*)[K(X,  t)  -  F(X'.  t)]p  -  0.  (5) 

Solutions  yield  the  time  dependent  density  matrix 
p(X,  X',  r),  whose  diagonal  components  provide  the 
density,  and  whose  values  are  constrained  by  the 
integral:  |  d*  Xp(X,  X)  Af*,  where  A’*  is  the  total 
number  of  electrons.  Assuming  free  particle  con¬ 
ditions  along  the  Y  and  Z  directions,  the  density 
nutrix,  with  h^/lmk^T,  separates  and  we  seek 
PiX,X\ty. 

p(X,X'.r)  =  p(.ir,r.r) 

X  exp  -[[(Y-Yy  +  (.Z-  zy]m  (6) 

Here  2,  is  the  thermal  de  Broglie  wavelength. 
Equation  (5)  separates  and  the  X,X’  portion  b 
rewritten  in  terms  of  center  of  mass  and  nonlocal 
coordinates: 

center  of  mass  coordinates:  {X  +  X')/2  ■=  x; 

nonlocal  coordinates:  {X  —  X')P, »  (.  (7) 

Note:  the  transformation  b  consistent  with  [3],  but  b 
not  area  preserving  (the  Jacobian  b  rrat  unity).  In 
terms  of  these  variables  and  for  free  partide  con¬ 
ditions  along  the  other  directions,  the  governing 
equation  for  p(x  +  C,  x  —  r)  b: 

dp/dt  +(hl2mi)d^pldxdi 

-  (l/«)(K(x  K(x  -  C,  t)]p  -  0.  (8) 

All  results  arise  from  p(x  -I-  (,  x  —  {);  nevertheless, 
we  require  expressions  for  current  and  energy,  whidi 
are  obtained  from  the  diagonal  components  of  the 
following  matrices: 

density:  p{x  -H  (,  x  -  C);  (9a) 


currait  density:  y(x  -I-  (,  x  -  ()  =  (h/(2i»i/)Jdp/dC; 

(9b) 

energy  density:  £(x  +  {,  x  —  ()  =  -(k’/8»i)5’p/dC*. 

(9c) 

The  above  definitions  are  discussed  in  [7],  and  in 
Appendix  B. 

THE  APPBOXIMATE  DENSITY  MATBIX  EQUATION  AND 
EQUIUBBIUM  SOLUTION 

Numerical  solutions  are  obtained  from  eqn  (8).  For 
interpretive  purposes  and  for  developing  the  QHD 
equations,  we  approximate  eqn  (8)  in  two  steps.  Fini, 
we  assume  an  infinitdy  diffmentiable  potential,  in 
which  case  eqn  (8)  becomes: 

dpidt  +  (hpmi)d^pldxdli  -(2/ih) 

X  Y.  fl/(2/  +  1)!K‘"*  '*^‘*'*  '^yidx^*  "P  -=  0,  (10) 

where  the  sum  b  over  0^1  <co.  Second,  we  retain 
only  the  first  two  terms  in  the  expansion: 

dpidt  +  ihf2mi)d^pldxdC 

-  (I/M)(2Cd^'/ax  -I-  (Cl3)d^Vldx*]p  -  0.  (1 1) 

Note:  retaining  only  the  term  linear  in  C,  yields  an 
equation  equivalent  to  the  time  dependent  collision- 
less  Boltzmarm  equation;  demonstrating  that  quan¬ 
tum  effects  arise  from  higher  order  terms  in  the 
expansion  of  (F(x-f-{,r)—  K/x —f,/)}.  For  the  den¬ 
sity  matrix  equivalent  to  the  collisionless  Boltzmarm 
equation,  and  for  dpidt  0; 

p(x  +i:,x  - C.  I)  =  p,exp  - +  mx)]  (12) 

b  an  exact  solution  for  free  partides  (no  collisions)  in 
a  potential  energy  distribution  K(x).  Mote  generally: 
p(X,  X'.  f)  -  p,exp  -  +  PK(x)Jexp  -  ({(K  - 

Yy + (Z  —  Z'f)l4X*f.  For  a  reference  density 
Po,  the  Fermi  energy  » (l//l]ln(pb/N],  where 
ti’*‘7(ml2Kfihy^.  Equation  (12)  b  equivalent  to 
^Pl—P{{p^Ptn)+  F(x)}]  (see  Appendix  B). 

To  obtain  the  quantum  modifications,  we  recog¬ 
nize  that  the  dassical  carrier  density  and  mean 
kinetic  energy  density  undo’  zero  current  conditions 
ate  respectively:  p(x,x)*psexp-{/lK(x)],  and 
£(x,x)mc(x,x)p(x,x)-p(x,x)1^7/2,  where 
e(x,  x)  b  the  mean  kmetlc  energy  per  particle,  and  that 
eqn  (12)  can  be  recast  aK 

p(x  -l-C,  X  - C)  -  p(jc,  x)exp-(2C*/k(x.  x)/2*l.  (13) 

Equations  (13)  and  (12)  have  the  sanw  content  for 
dasdcal  transport  For  quantum  transport  the  mean 
kinetic  energy  per  partide  indudes  m^fications  to 
the  classical  vdue(3].  The  numerical  studies  below 
suggest  that  the  effect  of  the  quantum  correction  b  to 
dt^  pinch  or  stretch  the  d^ity  matrix  along  the 
nonlo^  direction.  Equation  (13)  represrats  both 
contributions.  To  obtain  these  corrections  eqn  (13)  b 
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substituted  into  eqn  (11)  with  the  resulting  equation 
ordered  in  powers  of  { : 

C{4dit(x,  x)p(x,  x)ydx  +  2(dyidx)p(x,  x)} 

-  (8{ x)d€(x,  x)ldx 

-  (X^I24)d^yidx^)p(x,  X)}  =  0.  (14) 

Thus  separately: 

25(t(x,  x)p(x,  x)]/dx  +  {dyidx)p(x,  x)  -  0  (15) 

€(x.  x)d€(x,  x)ldx  -  (X^I24fi)d^yidx^  -  0.  (16) 

Equation  (16)  submits  to  an  immediate  solution: 
t(x,  x) »  f«|l  +  (>l/£,)Hl/12P)a*l'/ax*l‘'»,  when 

€^’=k,,TI2  is  independent  of  position  and  the  inte¬ 
gration  constant  is  chosen  to  retrieve  the  classical 
result  under  uniform  field  and  density  conditions.  If 
the  quantum  corrections  are  small  compared  to 
the  quantum  corrected  energy  is: 

£(x,  x)  s  £(x,  x)p(x,  x) 

-  (ikfc r/2 + (X*i\2)d^yiax^(x.  x).  (i?) 

which  corresponds  to  Wigner's  result  ([3],  eqn  (30)). 
For  the  quantum  corrected  density,  we  solve  a  re¬ 
arranged  eqn  (IS),  using  eqn  (17)  for  energy: 

Kx^md^y/dx^+ay/dx] 

l[(x^l6)(a^y/ax^ + im  4-  a  in  piax  -  o.  (is) 

For  small  quantum  modifications  the  above  inte¬ 
grates  (approximately)  to: 

p(x,x)-p,(exp-/J(K  +  fiw/3)l,  (19) 

where: 

Q^^(h^l4m)pia^yiax^-p(ayiaxYi2]  (20) 

defines  a  Wigner  quantum  potential.  For  smaU 
modifications  eqn  (19)  becomes:  p(x,x)» 
Pf[exp  — py{x)]{l  —PQwfi},  which  conesponds  to 
Wigner’s  equation  (2B). 

Equations  (19)  and  (4)  have  the  same  form 
although  the  modification  is  in  terms  of  potential 
energy  rather  than  density.  To  the  extent  that  the 
above  approximations  are  of  order  h^,  we  replace  the 
potentid  energy  in  eqn  (20)  with  its  classical  density 
equivalent:  ^K(x)— —  lnIp(x,x)/po]'  In  this  case 
Qb^Qw  >nd  eqn  (4)  is  retrieved  with  a  In  terms 
of  denrity,  the  energy  (eqn  (17)]  is  teexpressed  as: 

E(x,  X)  -  [*,  Tf2  -  (h^/2Am)aH\n  p)/a*x]p(x,  x). 

(21) 

The  quantum  corrected  form  of  the  density  matrix 
using  eqns  (13),  (17)  and  (19)  is: 

p(x  +  C.  X  -  ()  -  p,exp  -  I^{F(x)  +  Cw/3} 

+  (C/A)*{1  +  (2V6)/W*K/ax>}].  (22) 

For  small  corrections,  p(x  +  C,  x  —  C)  a  p#exp  —  [((/ 
xy  +  py(x)]{l  -  pQy,l3  -  (CV6)/W'F/ax»}.  which  as 
discussed  in  Appendix  B,  yields  upon  application  of 


the  Weyl  integral,  Wigner’s  form  of  the  quantum 
corrections  (eqn  (25)  of  (3]). 

Equation  (19)  highlights  the  quantum  modifi¬ 
cations.  For  example,  in  the  case  of  a  symmetric 
barrier,  classical  thrary  teaches  that  density  is  deter¬ 
mined  solely  by  the  value  of  the  potential.  Quantum 
theory  is  predicated  upon  continuity  of  the  wave 
functions,  permits  tuimeling,  and  teaches  that  the 
density  within  a  barrier  can  be  higher  than  that 
determined  classically.  At  the  peak  of  the  barrier, 
K.  ~  0,  <  0, 6w  <  0  and  the  dennty  exceeds  its 
classical  value.  Within  a  symmetric  quantum  well, 
at  the  center  of  symmetry,  —  0,  >  0,  ^  >  0 

and  toe  density  can  be  less  than  that  obtained  dassi- 
cally. 

COMPAUSON  OF  EXACT  AND  APPROXIMATE 
EQUIUBRIUM  DISTRIBUTION  FUNCTIONS 

The  extent  to  which  quantum  modification  rep¬ 
resent  quantum  transport  in  structures  under  equi¬ 
librium  was  addressed  in  two  stq>s.  First,  solutions 
were  obtained  for  the  coupled  LiouviUe  equation  (8) 
and  Poisson’s  equation: 

ai€(x)ayiaxy3x  -  -eip(x.  x)  -  pd(x)1.  (23) 

(In  the  discussion  below,  the  permittivity  and  effective 
mass  are  constant,  with  values  are  those  appropriate 
to  GaAs.)  Second,  the  exact  density  computed  from 
the  LiouviUe  equation  was  inserted  into  eqn  (2), 
(2t  was  computed,  and  an  approximate  density 
and  energy  per  particle  was  obtained.  The  results 
of  part  one  and  part  two  were  compared.  In  aU 
calculations  global  charge  neutrality  oocured: 
{dx|p(x,x)  — pi>(x)]»0.  Two  examples  were  con¬ 
sidered,  each  at  300 IC,  where  for  GaAs  the  thermal 
de  Broglie  wavehmgth  is  2  —  45  A.  In  both  calcu¬ 
lations  the  nominal  density  was  10'*/cm’.  (At  these 
densities  galUum  arsenide  calculations  necessitate  the 
use  of  Fermi  statistics.  These  have  been  performed  by 
the  authors[10],  and  demonstrate  two  density  depen¬ 
dent  contributions  to  energy,  one  classical  and  a 
second  quantum  mechanical  in  origin.  At  lower  den¬ 
sities  where  Fermi  statistics  ate  not  an  issue  calcu¬ 
lations  demonstrate  that  the  effects  of  the  quantum 
potential  are  qualitativdy  similar  to  the  results 
discussed  below,  but  the  magnitudes  of  the  density 
derivatives  are  smaller  Oonger  Debye  length)  and  the 
quantum  corrections  are  reduced.) 

Classical  N*N~N*  structures 

The  structure  is  1600  A  long  with  a  nominal  doping 
of  lO^/cm’  and  a  centrally  placed  500  A,  10'*/cm’ 
region.  The  variation  in  background  doping  was 
over  one  grid  point  or  4  A.  Solutions  yield  p{X,  X'), 
which  in  equilibrium  is  real  and  symmetric, 
p{X,  X')  B  p(X\  JO,  Bs  displayed  in  Fig.  1(a).  The 
inset  to  Fig.  1(a)  is  the  free  particle  density  matrix. 
Free  particle  Boltzmann  boundary  conditions  are  as¬ 
sumed;  i.e.  p(X,  X')  po exp  —  (C/Xy.  All  numerical 
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Fig.  1.  (a)  Density  matfia  for  the  N*S~N*  structure  with  free  particle  boundary  conditions,  as  obtained 
firam  the  UouviUe  equation  (Jt).  The  physical  dimension  of  the  structure  is  1600  A,  requiring  that  the 
density  matrix,  which  is  calculated  over  a  square  matrix,  is  of  side  1600  A/^2.  The  center  of  mass  and 

nonlocal  coordinates  are  indkatfed.  The  inset  is  the  free  particle  density  matrix,  (b)  ( - )  Diagonal 

component  of  the  density  matrix  from  (a).  ( — Density  as  obtained  from  eqn  (4)  with  o  >  0;  ( - )  with 

e  1/3.  (c)  Self  consistent  potential  energy  for  the  calculation  of  (a),  (d)  (^ntum  potential  for  the 

calculation  of  (a),  (e)  ( - Energy  per  particle  as  obtained  from  the  diagonal  component  of  eqn  (9c) 

for  the  calculation  of  (a).  ( - ^)  Ene^  per  particle  as  obtained  from  eqn  (24)  with  a  —  1/3;  ( - )  with 

a  =  1/4;  (— )  with  a  »  1/5. 


to  the  exact  solution  occur  for  a  between  1/3  and  1/4. 
The  significance  of  these  results  is  that  the  quantum 
contributions  are  soMy  responsible  for  the  spatial 
variation  in  the  energy  per  partkie,  and  demonstrates 
the  presence  of  quantum  contributions  with  classical 
structures. 

Singk  barrier  diodet 

For  this  calculatkm  the  background  density  is  flat 
and  equal  to  lO'Vcm*;  the  structure  is  2000  A  long 
and  the  grid  spacing  is  uinform  and  equal  to  3.33  A. 
Figure  2  displaya  the.  ic^ts  for  a  iOiOnieV  barrier 
represented  anidytically  by: 

>'*«ur(7f)-500ineVI{l  +  tanh{(x  -4r,)/hJ}/2 

+  {1  -  tanh((jc  -  o,)/h]}/2  ~  1],  (25) 

where  a, -—ISO  A,  02  =150  A,  50  A/3.80. 
Figure  2(a)  displays  eqn  (25),  where  increases 


continuously  (from  near  zero)  to  500  meV,  over 
approx.  75  A. 

Figure  2(b)  displays  p{X,X').  As  in  the  Fig.  1 
calculation,  free  particle  boundary  conditions  are 
assumed.  The  dramatic  “hole"  is  a  consequence  of 
the  barrier.  Figure  2(c)  is  a  line  plot  of  density.  The 
solid  line  is  obtained  from  the  liouvifie  and  Poisson 
equations;  the  dashed  lines  are  from  eqn  (4)  with 
a  =  0  (long  dashed  line)  and  1/3  (short  dashed  line). 
Common  to  each  calculation  is  a  significant  reduction 
of  duuge  within  the  barrier,  as  well  as  charge  ac¬ 
cumulation  on  either  side  of  the  barrier.  At  the  edges 
of  the  barrier  the  solutions  closest  to  the  Liouville 
equation  are  those  for  a  =  1/3.  The  reduction 
of  charge  within  the  barrier  is  a  consequence  of 
the  barrier,  while  the  presence  of  charge  adjacent 
to  but  outside  of  the  barrier  is  a  consequence  of 
self-consistency  in  the  calculation.  Its  magnitude  is 
dependent  on  the  condition  of  global  charge  neu¬ 
trality.  Figure  2(d)  displays  the  potential  energy 
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Fig.  2.  (a)  Sketdi  of  the  hyperbolic  tangent  barrier  centrally  placed  within  the  2000  A  structure,  (b) 
Density  matrix  for  the  single  barrier  structure  with  free  partide  boundary  conditions,  as  obtained  from 
the  Liouville  equation  (8).  The  physical  dimension  of  the  structure  is  2000  A,  requiring  that  the  density 
matrix,  which  is  calculated  over  a  square  matrix,  is  of  «de  UXOkly/t.  The  center  of  mass  and  nonlocal 

coordinates  are  indicated,  (c)  ( - Diagonal  cooponeat  of  the  density  matrix  from  (a).  (— )  Density 

as  obtained  from  eqn  (4)  with  a  - 1/3;  ( - with  a  >0.  (d)  Seif  consistent  potential  energy  for  the 

calculation  of  (a),  (e)  Qiutum  potential  for  the  cakuiation  of  (a).  Inset  includes  (d).  (f)  ( - )  Energy 

per  partide  as  obtained  firom  eqn  (9c).  ( — )  Energy  per  partide  as  obtained  from  eqn  ^4)  with  a.m  |/3. 


THE  APPROXIMATE  NONEQUIUBRIUM  DENSITY 
MATRIX:  THE  CONSTRUCTION  OF  THE  QHD 
EQUATIONS 

The  nonequilibrium  situatioii,  is  considered  within 
the  framework  of  the  QHD  equations,  which  incor¬ 
porate  quantum  contri^tions  as  mor^atUms.  The 
QHD  equations  incude  dissipation  within  the  context 


of  Fokker-Planck  (FP)  scattering.  The  motivation 
for  FP  dissipation  is  simplidty.  When  scattering  is 
treated  as  in  the  Boltzmann  tran^xnt  equation, 
utilization  of  the  Weyl  transformation  results  in  an 
equivalent  scattering  integral,  that  is  approximatdy 
dependent  upon  powers  of  (,  and  derivatives  with 
respect  to  C[I0]-  Under  special  drcumstances  these 
take  the  form  of  FP  dissipation.  The  equation  of 
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Fig.  3.  (a)  200  meV  ttanier  within  the  6000  A  structure  studied  with  the  siinulation.  (b)  Density  nwtiix 
for  the  single  batiier  structure  with  free  partide  boundary  conditions,  as  obtained  firm  the  Liouville 
equation  (8)  for  a  bias  of  —0.2  meV.  The  physical  dimension  of  the  structure  is  6000 A,  requiring  that 
the  density  matrix,  which  is  adoulated  over  a  square  matrix,  is  of  side  6000  A/y2.  The  center  of  mast 

and  nonlocal  coordinates  are  indicated.  (g)( - )Sdf  consistent  diagonal  component  of  the  dcnsi^  matrix 

forabiasof  — 0.2meV.( - Quantum  corrected  solution.  (d)( - )Sctfcontistent  diagonal  component 

of  the  density  matrix  for  a  bim  of  O.OmeV.  ( - )  Quantum  collected  solution,  (e)  ( - )  Sdf  consistent 

potential  energy  for  a  bias  of  —0.2  meV.  ( - )  Quantum  corrected  solution.  (0  ( - )  Seif  consistent 

potential  energy  for  a  bias  of  0.0  meV.  ( - )  Quantum  corrected  sdution. 


cated  by  representing  nonequilibrium  by  a  displaced 
Maxwellian,  exp  —  —  Pifl2m  -F  V}],  and  where 

the  mean  momentum,  p^,  the  density,  and  a  particle 
temperature,  are  to  be  determined.  The  argument  for 
a  displaced  Maxwellian  is  the  assumption  of  rapid 
thermalization.  While  there  is  experimental  evidence 
that  some  quantum  feature  size  devices  sustain  strong 
relaxation  effects,  such  phenomena  is  not  likely  to  be 
universal.  Nevotheless,  as  a  first  stq>  in  dev^piiig 
a  set  of  nonequiUbrium  QHD  equations  we  examiiK 
the  consequences  of  mo^fying  the  quantum  equi¬ 
librium  distributions  to  describe  nonequihbrium 
conditions.  Within  the  context  of  the  coordinate 
rqiresentation,  the  Weyl  transformation  as  discussed 
in  the  Appendix  B,  di^tes  that  the  displaced  equi- 
fibrium  d^ty  matrix  (generally  non-Maxwellian) 
used  below,  is  obtained  through  tte  following  modifi¬ 
cation  of  a  zero  current  denaty  matrix: 

p{x  +C,x  -C)»p(x  -K,  Jc  -  C)exp  +  lUptCIh],  (28) 


where  is  at  most  a  function  of  x.  With  the  current 
incorporated  as  in  eqn  (28)  the  construction  of  the 
QHD  equations  proceeds  in  three  parts:  First,  the 
truncated  density  matrix  is  identified  as  eqn  (28)  with 
the  form  of  the  equiUbrium  contribution  given  by 
eqn  (22)  (the  potential  is  replaced  by  the  dassical 
d^ty  equivalent);  second,  the  relent  transport 
quantities  are  idratified  as  carrier  density,  p(x), 
mean  momentum,  pt{x)  and  dectron. temperature, 
T^{x) »  third,  the  moment  equations,  are 

obtained  from  a  succesrion  of  derivatives,  followed 
by  the  limit  as  (  ■»  0.  In  taking  moments  we  note  that 
mudi  information  contained  in  the  off-diagonal  el¬ 
ements  of  the  density  matrix  is  lost 
With  eqn  (28)  Um  following  quantities  (from 
eqn  (9)]  are  relevant  to  the  moment  equations  (with¬ 
out  the  equipartition  contributions  of  the  T  and  Z 
directions): 

J{x,x)^p(x,x)pJm; 


(29a) 
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E(x,x)~[p\l2m+k^np. 

-{h^llAm)d^i)np)ldx^\  (29b) 

l^\x,  x)  -  [/»3  +  iptmk^,T^ 

-P4(fty*)S^tip)ldx*]p.  (29c) 

Equation  (29c)  is  the  diagonal  component  of 
P'Hx  +  C.  X  -  C)  “  (h/2i)’dV/^’.  «“<*  represents 
the  energy  flux,  (as  typically  appears,  e.g.  in  the  third 
moment  of  the  Boltzmann  transport  equation). 
Equations  (29)  and  their  dependence  on  d^vatives 
of  density  are  valid  only  in  the  limits  discussed 
in  the  above  sections,  and  rqtresent  mod^icatuxu  to 
dassical  situations.  In  this  sense  it  is  important  to 
note  that  the  derbxOion  of  the  quantum  potential 
in  terms  of  Qei  ex|^tly  involved  the  carrier  tempera¬ 
ture.  The  Bohm  potential  Q^,  is  independent  of 
electron  temperature.  The  consequences  of  using 
rather  than  Q^,  in  the  QHD  equations  should  be 
examined. 

The  QHD  equations  are  obtained  by  taking  succes¬ 
sive  derivatives  with  leqiect  to  C,  as  ddfined  by  eqn  (9) 
and  taking  the  Kmit  C  -»0.  The  QHD  particle,  momen¬ 
tum  and  energy  balance  equations,  are  respectively: 

dp/dt  dlppjmydx  ->  0;  (30) 

+  2d£(x,  x)/dx  (dVldx)p  +  ppjx  «  0; 

(31) 

dEjdt  +  l/(2m*)d/x>Vax 

+  (ppJmWISx  +  2£/t  -  (S/m)p  =  0.  (32) 

We  rearrange  eqns  (31)  and  (32X  noting  that  the 
quantum  correction  driving  force  is  implicit  in  E(x,  x) 
utd  P^x,  x).  Using  eqn  (29b)  for  £(x,  x)  and  noting 
that  d{pdH\np)ldx*ydx  ^ —(4mlh*)p^ffdx,  the 
QHD  momentum  equation  is(2]: 

^(PPeV^t  *1  S(ppym)dx  +  d(pkT)ldx 

+  pd(Qsmidx  +  pdV/dx  +  ppjx  =  0,  (33) 

which  differs  from  its  classical  analog  through  the 
presence  of  6i[2].  When  the  first  two  terms  are 
zero,  and  the  dedron  temperature  is  spatially  inde¬ 
pendent,  the  drift  momentum  density  reduces  to: 
ppt  -  -xk^TpdHV  +  Qt/3)lk^T  +  ln(p)l/5x.  Then 
for  a  K  1/3;  and  J  ■■  —eppjm,  eqn  (3)  is  retrieved; 
for  ^daiO,  the  density,  as  given  by  eqn  (4)  is  a 
solution  to  eqn  (33).  Note:  the  form  of  the  scattering 
term  in  eqn  (33)  identifies  the  first  part  of  the  FP 
scattering  as  a  frictional  term  (see  [7]). 

For  the  energy  balance  equation,  using  eqn  (29X 
eqn  (32)  becomes; 

dE/dt  +  d{(ptlm)lE  -b  (pim  -  li^/€idH\n  p)l 

X  dx*)]}/dx  +  {ppelm)V„  +  2£/t  -  (S/i»i)p  -  0. 

(34) 

To  determine  S,  we  note  that  it  generally  depends 
upon  X,  as  does  t.  In  the  context  of  eqn  (34)  we 
require  that  £  relax  to  which  is  the  »  0  value 


given  by  eqn  (21).  This  is  guaranteed  with 
S  K  TmE^/x.  Thus  eqn  (34)  becomes: 

dEldt  -b  d[ipJm)(E  +  pk*  T)\dx 

+  (ppJm)dlQ,l3-i-yydx 

-p(X%TI6)ldH!tap)iex^ 

x{pJm)ldx+Cl/x)[E-E,]~0.  (35) 

The  second  part  of  the  FP  dissipation  involves  a 
relaxation  to  a  non-zero  thermal  energy.  E,  above  is 
the  same  as  used  by  Woolard  et  n/.(13]. 

The  consequences  of  the  above  approxrrrutioos  is 
the  appearance  of  the  qurmtum  potential  with  the 
factor  “l/3“.  (The  situation  for  Fermi  statistics  is  not 
addressed  here.) 

QUANTUM  MOMENT  EQUATION  COMPUTATIONS; 

OOMPAWBON  TO  THE  EXACT  SOUHIONS 

The  development  of  the  QHD  equations,  is  pre¬ 
dicted  on  future  use  in  the  design  and  understanding 
of  multi-dimensional  quantum  feature  size  devices. 
The  degree  to  which  this  is  useful  renaains  to  be 
determined  for  nonequilibrium  phenomeiui,  and  the 
work  of  [4,5]  represents  an  important  beginning, 
Another  relevant  case  is  the  evaluation  of  density 
across  an  abrupt  heterostructure  region,  as  occurs  in 
either  a  heterostructure  diode  or  in  nuxlulation 
doped  FETs.  While  the  dieet  charge  density  can  be 
obtained  from  solutions  to  Schrodinger's  equation, 
the  incorporation  of  sudi  a  calculation  in  a  quantum 
corrected  standard  set  rff  device  simulation  equations 
has  only  recently  been  addressed.  We  consider  this  in 
assessing  solutions  of  the  QHD  equations  against  the 
Liouville  equation  in  the  zero  current  fimit  It  is  rwted 
that  the  use  of  an  abrupt  interface  violates  the 
following  conditions  regard  as  the  basis  for  the 
development  of  the  quantum  nradificalions:  the  po¬ 
tential  is  continuous,  and  the  value  of  Q  is  small  enough 
to  be  regarded  as  a  “correction”.  It  may  be  ooryec- 
tured  that  the  use  of  quantum  potential  has  more 
generality  than  that  uncovered  in  the  above  deri¬ 
vations;  at  this  time  there  is  no  justification  for  this 
claim. 

The  computation  is  for  a  6000  A  structure  with 
constant  10'*/cm*  doping.  The  grid  spacing  for  the 
Liouville  equation  was  constant  and  equal  to  7.SA; 
the  grid  was  nonuniformly  graced  for  the  QHD 
calculation.  A  200  meV  abn^  barrier  is  plaoed 
across  the  right  half  of  the  structure,  as  shown  in 
Fig.  3(a).  The  self-consisteot  i^iaoe  charge  prt^les 
were  conqmted  for  two  values  of  applied  bias: 
F’^fM^O.OeV  and  — 0.2eV.  In  both  computations 
the  quantum  potential  was  finite  within  the  vicinity  of 
the  interftce,  with  structure  similar  to  that  rff  the 
barrier  problem  discussed  in  Fig.  2;  it  was  zero  within 
the  vicinity  of  the  boundaries.  The  two  dimensional 
zero  current  density  rrutrix  for  y  »  — 0.2eV  is 
shown  in  Fig.  3(b). 
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line  ploto  of  Che  density  matrix  for  both  the  exact 
(solid  fine)  and  quantum  corrected  solution  (dashed 
line)  are  diq}layed  in  Fig.  3(c^).  Note  the  accumu¬ 
lation  of  charge  on  the  narrow  haiwtgap  side  of  the 
structure  followed  by  defdedon  (with  noiMiegligible 
values)  of  diarge  <m  the  wide  bandgap  portion  of  the 
structure.  Under  bias  the  edge  of  the  strudure  is 
depleted  of  charge.  The  potential  energy  distribution 
for  the  two  values  of  bias  are  diqdayed  in  IHg.  3(e/) 
where  we  see  that  the  discontinuity  in  potential  is 
equal  to  the  full  200  meV  associated  with  the  barrier. 
Tte  character  of  this  solution  is  similar  to  that  at  the 
edge  of  the  hyperbolic  tangent  barrier  shown  in 
Fig.  2.  In  particular  the  quantum  potential  is  positive 
(negative)  to  the  left  (ri^)  of  the  metallutgical 
boundary.  The  oomparatrve  density  and  potential 
profiles  are  extremely  dose  and  attest  to  the  oonfi- 
denee  of  the  appcoxitnation,  for  this  type  of  structure. 
But  caution  is  in  order!  The  exceOent  agreement  for 
the  calculations  of  Fig.  3,  bur  the  less  certain  agree¬ 
ment  of  Rg.  2,  indicate  that  a  c  '  did  case-by-case 
asseaaement  may  be  necessary.  Nevertheless  it  ap¬ 
pears  that  obtahiing  representative  charge  densities 
necessitates  the  incorporation  of  quantum  effects 
through  such  additions  as  the  Bohm  quantum  poten- 
tiaL  Altetnatively,  realistic  device  simulations  must 
resort  to  a  full  multidimensional  quantum  tranqwrt 
calculation. 


SUMMAnV 

This  study  assessed  the  introduction  of  quantum 
modifications  of  classical  tranqiort,  with  the  results 
indicating  that  quantum  corrective  transport  is  useful 
under  certain  circumstances,  and  that  many  simple 
device  studies,  tud>  as  those  for  HEMTs  would 
benefit  from  its  incorporation.  It  is  likdy  that  such 
corrective  transport  considerations  would  also  be 
valuaUe  under  nonequilibiium  conditions  particu¬ 
larly  in  evaluating  transport  across  beterostructure 
regions.  It  is  important  to  note  that  introduction  of 
the  quantum  potential  in  a  generic  form  of  the  QHD 
equatiotts  is  not  new.  it  has  been  Imked  to  density 
functional  theory,  as  discussed  by  Deb  and  Ghosh(i4] 
who  also  identify  the  force  as  being  of  quantum 
origin.  B^im  and  Hiley(IS],  point  out  that  an  essen¬ 
tial  new  feature  of  the  quantum  potential  is  that  for 
sinile  particle  Scfarodinger  fidds,  only  the>bnN  of  the 
Scfarodinger  field  counts,  not  the  mtensity.  The  force 
arising  from  this  potential  is  not  like  a  mechanical 
force  a  wave  ptndiing  on  a  particle  with  a  pressure 
proportional  to  the  wave  intensity;  rather  the  force 
arises  from  information  content,  e.g.  structure,  rather 
than  value,  of  the  wave(lS].  Bohm  and  Hiley[15] 
distinguish  this  force  from  tte  Maddungfl6|  hydro- 
dynamic  model  in  which  the  particle  is  pudwd  me- 
diankally  by  the  fluid. 
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APPENDIX  A 
Sobitkm  Procedure 

For  convenience  of  aohitioo  and  determining  auitable  forme 
of  boundary  confitioiia,  «pi  (8),  it  rawritten  as  a  coupled 
Brat  order  ayitem  of  equatioiu(i7]: 

u(X,X’)■^■yR|^^n\^Ol^X■^■^p|^X^’~i|  (Al) 

BpIdt-i-^ldX-dufiXI 

+  0  -  y(X\  Dip  -  0.  (A2) 

Etiuation  (Al)  deBnet  «(Jr,  JT', /);  eqn  (A2)  is  an  ahernative 
form  ofeqnff)  after  acnounling  for  fora  partidecooditioiia 
along  the  y  and  Zdiractiona;  rewritten  in  terma  of  u  and  p. 
The  chancletittk  dhactiona  (or  oqna  (Al)  aad  (A2)  are; 

x-fr-t-jrT/Z-cooatant  (A3) 

C-(y-y')/2-conatant  (A4) 

In  tetma  of  the  characteriatic  directiooa  x  and  C.  oqna  (Al) 
and  (A2)  can  be  written  aa: 

u(X.  X')  4-  (A/2m]dp/dx  -  0  (A3) 

ip  fit  +  dufiH  +  p/»I(F(Jr.  0  -  ViX',  /)]p  -  0.  (AQ 

fontabie  boundary  conditiona  for  eqna  (A3)  and  (A6)  are 
the  apeciBcatfon  of  p  and  a  along  the  botmdary  JIT' 0  and 
the  apecification  of  «  along  the  boundary  y  B  L/.yZ  where 
L  ia  the  length  of  the  device.  Along  the  boundary  y  •O,  p 
it  apeetted  aa  the  complex  conjugate  of  p(y,0),  aince  p  is 


hermitma,  and  a  it  computed  form  the  outgoing  character- 
itiie  eqn  (A6).  Along  the  boundary  X  -  L/^2,  p  it  com¬ 
puted  (ram  the  outgoing  characteristic  eqn  (A3). 

Aa  alternative  tyalem  can  be  formulalad  m  terma  of  the 
current  matrix: 

AX.  y ')  +  (A/2ai J9p/aC  -  0  (AT) 

ip  fit  +  djfix + (i7*)(y(y.  1)  -  y(X’.  /Hp  -  o.  (A8) 

Eqiwtioos  (AT)  and  (A8)  have  the  tame  characterittic 
diredioot  x  aial  (  at  equations  (Al)  aad  (A2).  Suitable 
boundary  conditions  for  eqat  (AT)  and  (A8)  are  the  apedB- 
cation  of  p  and  y  along  the  bomdary  y'^O  arid  the 
iqteciftration  of  p  along  the  boundary  y -E  Along  the 
boundary  y  0,  y  is  specified  at  the  conqilex  conjugate  of 
y(x.O)  and  p  it  oonqiM^  from  the  outgoiag  characteriatic 
eqn  (AT).  Along  the  boundary  XmL,J  is  computed  from 
the  outgoing  characteristic  eqn  (AS).  Both  seta  of  the  Brat 
order  system  of  equatioot.  eqna  (A3)  and  (A6)  and  eqnt 
(AT)  and  (AS),  are  useful  in  appUcationt  tinoe  they  allow 
dilimt  forma  of  boundary  conditiont.  Both  seta  of 
equatioot  can  be  solved  by  the  tame  numerical  procedure. 

The  tohitioo  procedure  oontittt  of  solving  the  first  order 
^stm  of  equatioot  at  an  inilial  boundary-value  probta 
starting  fiom  conditiont  along  the  line  y*-0  and  marchiag 
to  the  Ime  y '  L  using  the  method  of  chancterialka.  A 
characterittic  net  for  the  equation  of  motioo  of  the  density 
matrix  can  be  constructed  e  priori  from  grid  poirua  of  a 
uniform  square  grid.  A  diaerete  form  of  eqnt  (A3)  and  (A6) 
on  this  grid  h  ^ig.  4 

IiiL  +  lA/2«lI<*(Ly)-i»(f-l.y-01/Ax-0  (A9) 

ldp/*L  +  +  l.y  - 1)  -  rKLya/AC 

+ p/»jiK(y.  f)  -  V(X\  /)LbL  -  fl.  (AlO) 

where  (  L.  tepreaents  an  average  over  foe  grid  cell.  Depend¬ 
ing  upon  foe  form  of  averaging  choaen,  eqnt  (A9)  aad  (AlO) 
form  a  ^stem  of  2  X  2  blo^  tridiagonal  or  blo^  diaipmal 
algebraic  equatfona  font  can  be  aoM  at  y' from  known 
vahMt  at  JT  ay  - 1.  Thus,  foe  aolutioo  proMdure  can  be 
marched  from  boundary  conditiona  at  y*»0.  in  tiepe  along 
X‘,  toX’^L  Similar  prooedorm  can  be  utflnsid  frir  eqna 
(AT)  and  (AS). 

S^-oontiatenqr  it  included  in  foe  analyaia  by  iterating  foe 
aohition  of  foe  deatHy  matrix  equation  with  foe  solution  of 
PotaMo’s  equation  to  oonvergenoe,  by  aucceaaive  tubati- 
tinion.  For  this  purpoae,  Foiaton’t  equation  it  written  in  foe 
form: 

{i[[{x)i{Eyyixyixy + dfipfito^v 

-  -  {i\fix)ivfixyixy  -  dfpix.  X)  -  p*(xH-.  (Al  1) 

with  K"*'  m  K*  4- AF*.  where  a  it  the  iteration  number.  The 
second  term  on  foe  left  ha«l  aide  of  equation  (Al  I)  serves 
to  acederate  oonvergenoe  of  foe  iteration,  wherein  ipfiV  it 
evaluated  at  X  eillier  numerically  from  previous  itentiont  or 
analytically  aa  ipfiV  >  —p{x,x)fk^T,  for  Boltzmaim  atat- 
iatict.  A  S^roint  ociMeted  Baile  dUfaenoe  approximation  to 
(Al  1)  teauto  in  a  tridiagonai  rystem  of  algebraic  eipiationa 
that  can  be  atdved  eaaily  aad  rArientfy  for  AF,  sdiiA  it  foe 
incremeat  in  V  between  Heratioiis. 

The  first  step  of  the  iteration  procedure  oonaiata  of 
aaauming  a  dittributioo  for  the  self  oooaiatent  potential 
(QfpicaUy,  aero  everyediete)  and  sotving  the  deaaiQr  matrix 
equation  to  obtain  the  de^ty  distribution.  Baaed  on  the 
computed  density,  Foitaoo't  equation  (All)  b  solved  to 
update  foe  adf  oontitleat  potmtial  diriribution.  For  foe 
compotationt  of  this  paper,  foe  analytical  expieaaion  for 
ipfiV  was  utiliaed.  For  cases  rvhete  ipfiV  it  computed 
numerically,  several  Uerationa  (typically,  four  or  five)  ate 
requited  before  a  reliable  estimate  for  foe  gradient  can  be 
computed.  During  these  initial  iterations,  second  term 
in  equation  (All)  is  replaoed  by  a  term  of  foe  form 
— (c(x)AF/(At{Ax}^]  (Ax  b  foe  mesh  spacing,  Ar  «  30)  for 


TraMport  dMOfy  via  LiotfvtUe  eqitalioa 
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pnnvargwwTi  Hm  kna  coidd  be  atiliKd  Cor  all  iiastioM. 
btWMwanwcai»aott«pid.Sol>nioBieMaMii’tw|«iatkw 
ii  Had  H  apdala  dw  mM  cawdilwt  polaHial  baaed  opoa 
whiob  dw  dcHily  aaitm  aquation  it  toived  apm.  The 
itefidoH  aic  rapaUed  aadl  Ite  dcMity  aad  potential  dittri- 
badoH  coBuefgt  to  the  adT  coaiwlrBt  aohitioa.  For  the 
compatttioM  ptHtalad  here,  tia  oidert  of  veadual  re- 
ducdon  WH  obtained  in  leat  than  10  itcrationt. 


APPENDDCB 

MaihH  of  RoaoUs  to  the  Wigm  Formulation 

The  oonnection  betiven  the  deanty  BMtiiz  in  the  coordi¬ 
nate  lepieeentation  and  the  Wgaer  Amotioa  it  throngh  the 
Wcyl-t^  ttaatfomatioH  with  nonnaiiaatioM  peraliar  to 
the  pfoMem  of  inteiatL  For  the  deuity  matrix: 

p(x  +  C.x-C)-2(li(2**)l’ 

X  J"  dV.(p.*)expP/p  {/hJ.  (Bl) 

where  the  factor  of  2  aooountt  for  the  &et  that  each 
ammentianttaie  can  hold  two  electroni.  The  invitaetfaiia- 
fbnaation  it: 

/.(p.w)-2*/2j*  d>0»(x  +  {.x-{)«qif-21p  C/hl.  (B2) 

where  the  foctor  2*  it  a  oonaequenoe  of  the  defittitian  of  the 
nonlocal  cooedinate  (tee  eqa  PS- 1»  (hit  tiaiitforaiatioo  it 
it  atteried  that  the  Wigner  Amctien  aad  ail  neocataiy 
derivativtt  with  wtpect  to  momentum  vaniA  u  p»  ±  ao. 
NotK  (a)  p(x.x)w(l/(2a*N*f?«d'^(p.w);  (b)  aiAtti- 
tudon  of  eqn  (12)  into  eqn  (B2)  yielda  the  ceaaltt: 
/. +  F(x)  -  E,]\. 

The  Wigner  i^untiom  indudipg  FP  tcattering,  at  dttcuned 
by  Sttoadolll]  it: 


+(I/B)(l/2a*) 


djr70»'.Jt) 


X  (F(x,  I)  -  F(x'.  (»expCi(it  -plx'/k] 
-(l/T)divJin+5Vp«/. 


(B3) 


where.  H  m  the  main  text,  all  apatial  variatioH  are  along 
the  xdinction,  Boknnaan  ttatittict  apply,  and  momentum 
variatioH  in  all  three  dimfiitioM  are  allowed.  The  oo- 
eCcientt  tend  5  are  choaw  Min  the  dentity  matrix  ttudiet. 
For  traatpoft  in  one  tpooe  dmention  it  it  direct  to  demoa- 
atiale  that  the  iniepal  in  eqn  (B3)  redooet  m  the  cUarical 
cate  to  (dK/dx)(4f/dp).  To  aeoond  order  in  h,  the  Wigner 
eqmdon: 

+  (PimWIi*  -  iVIdxdfIdp 


+(ii*/34)(d»F/at*)aW-(iA)«v,(in+H^  (W) 

The  left  hard  tide  of  eqn  (B2)  hM  bem  ditcumed  in  ((.iq. 
AppBcaticnofttetiiMformation.  eqn  (Bl).  yielda  eqn  (11). 

b  the  afaeenoe  of  djaripation  be  approximate  Wigner 
dittiibution  Amctioo  to  aeoood  order  in  h  it  (WignerPL  tee, 
c«cqn(2S)|: 

/.-exp-#lp*/2«  +  V{x%l 
X  ((a>K/dx»  -  HdVIdxfm  -  Hp^ftmyd^VIdx^}  (B5) 


which  upon  application  of  eqn  (Bl)  yieidt  the  equation 
fcdowiiig  eqn  (22). 

The  idatiM  betwem  the  dentity  matrix  and  the  Wigner 
function  exteadt  to  obeervaUet.  permit^  a  coodte  ded- 
nition  of  the  eteocieted  matrioet.  Defining  aortni  density, 
energy  dentity  and  third  moment  matricies  letpectively,  at: 

Kx  +  C.x-C)-2ll/(2ii*)H 

X  J  d’p(p/ni)f.(p.  x)exp(2ip  •  (/h) 

£(x+(.x-C)-2(l/(2x*)l» 

x|"  d>p(p  i/2ml/;(p,x)exp(2r>  (/fil 

pO|(x  +  C.x-{)-2(I/(2x*)p 

X  j"  d»/Wp  p)K.(p.x)exri2ip  C/hJ. 

it  it  direct  to  demonttrate  the  validity  of  eqnt  ^)  aixl  (%). 
The  derivatiwe  dcfinhiott  of  the  third  moment  [aee  diaeuttion 
foDowiag  eqn  Q9cH  follow  directly  from  the  above.  Note: 
If  be  ditttibinion  fimctioo  in  eqn  (Bl)  it  for  aero  curmit, 
and  a  finite  Gurreat  fonction  it  obtahied  ftom^(p  — 1^.  x), 
ben  the  aero  cuncat  and  finite  current  dea^  matrix  are 
rriated  m  eqn  Qt). 


(B6) 

(B7) 

(B8) 


APPENDIX  C 

Pure  State  Results  and  a  Comparison  to  Itrfrau, 
Grutm  and  Feny[8J 

The  pure  ttate  retuht  for  cuncat,  energy  dentity  and  third 
moment  may  be  obtained  m  foUowt.  Exptcat  the  wave 
function  M  P(x.l)«p(x.t)'^expi9(x,i),  wibp(x,()> 
hdOlix.  Then  Schrodinger’t  equation,  Mi’ldt  > 
--(h*/3m)9*P/9x*-i-  y(x,0'F,  which  k  ooaqgfex  it  rewrit¬ 
ten  H  two  real  partial  deferential  eqHtiont: 

dpfdt  ^  (^1) 

UppVit  My  q.  pdVIdx  m 0,  (C2) 

where 

jlT-ttmoix  lUilm  (O) 

£(*.  »V-  T  ■  P  wax  V  (C4) 

While  the  content  of  Scfarodinger’t  eqwtioo  it  contained  in 
etpit  (Cl)  and  (€2),  an  exptcation  for  the  time  dcprndwice 
of  the  energy  may  be  obtained  through  the  time  derivative 
of  eqn  (C4)  and  judkiow  ute  of  eqnt  (Cl)  and  (C2).  We 
find,  with; 

-lp’-(**/*){3ipa'(inpwax'+awax^}]p  (cs) 
hR.^  .  ,_Mf  +(l/2«*)ai^nn»rH-|/^* 

■i-(pplm)dVI»x~0.  (06) 

Note  the  difietenoet  between  the  pure  atate  definitiou 
{eipia  (O),  (C4),  aad  (CS)].  and  that  of  equ  (29).  For 
be  pure  ttate:  there  b  the  ebaence  of  a  temperatme  depen¬ 
dence,  the  factor  of  3  b  absent,  and  there  b  a  vdodty 
correction. 
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Abstnet 

TTiis  study  describes  the  evolution  and  imptementation  of  a  set  of  quantum  balance  equatioiii  for  cauninii^ 
transmit  im  mesoosoopk  structures. 

Keyword: 

Wigtier  functions,  quantum  potendab,  rpiantum  balance  equatioos. 

Introduction 

TUs  study  describes  the  evohition  and  implenientation  of  a  set  of  quantum  balance  equations  for  examiniiig 
mnsport  im  mesoosoopic  stnictuies.  The  study  is  motivated  fay  a  perceived  need  for  an  imnitiw^  accessible  set 
of  multi-dimensional  quantum  tiansport  equations,  that  perron  the  setf-eoniistent  calculation  of  pattide 
current  and  cunentdensiM.  The  yrol  is  the  devekyment  of  a  set  of  quantum  hydrodynantic  equations  that 
reduce  to  the  siqgie  paitiae  equations  [1]  lor  a  pure  state,  and  the  diusical  faydrodynamic  equmioos  12]  as  a -> 

0.  As  discussed  6elow,  these  goab  have  been  partially  met 

Pure  State  and  dasticti  Moment  Equations 


andadatticaIpotentulU(x,tXare,vritht(x,t)  «  Jpeq)(iS(x4yfll.andp,l  «  flS/flc 

*t^+*x(p<ip/hi)  “  0 

«t(^Pd)+»s<^Pd*/to)+^<x(U+Q)  -  0 
Q  - -(a«/2nUp)«,*Jp 


I  one  dimension. 


NR^Iecanse  of  the  one 
satitSed(3].)  Thedas 
•  are]^ 


represented  By  the  qMntumpoteaMl,Q  PL  Moterbecanse 
transport  and  Spatially  independent  effective  mast  are  CT; 
•^+*s(P<K^) 

*tCPPd)+*a(^Pd*^“)+^«aU-Hj|(pkT)  -  dAPd,Go]l 

-  flW^| 

W  -  3pkT/2+/>pd*/2m 


oitiire  of  the  tfinspof 


It  is  vrartlnriiile  emphasizing  that  tile 
with  natial  variations  in  only  one  dii 
The  derivative  notation  in  the  above 


equations  involve  tineedimenshwal  momentum  yace  imitation, 
and  that  dens^  and  momentum  now  rqnesent  patide  density, 
ms  is  fl  j  etc 


Structure  of  the  Quantum  Medianical  Equations 


this  we  turn  to  an  qyradmate 


function  obtained  Iqr  Wigner  [5],  and  diacoated  more  reooidy  by  Anemia  I 

-  eq)^{p*/2m-HJ]{l-2»(ax*U-d(«xU)*/2yV«»(l-dRi*/to)»x*U/3} 

In  equation  (8X^  >  l/kT.a  a*d*/2ni,andp*  ■■PK*'*‘Fy*'''Pii*- 

The  construction  of  the  noneonflifacinm  distribution  ftmetion  involves  re|^ 
in  equation  (8)  by  oorreqiondmg  density  eiytessions.  The  carrier  density  I 


hatr 

h  0(x>t).  Howfoodbll 
mmn,  ditcusted  in  [4]. 

r41  from  the  zero  cii^ 
ntiy  by  Anemia  udufii 


tnlnm  contribotions  win 
statement?  Toesamine 


I  distribution 


/  -  2(l/2aa)*/f^zy)d*p  -  Nezif^UIl-2<»(«,*U^(dxU*y2V3J  (9) 

whereN  <■  .After  demonstrating  that  4.11  ~  -(dj^VC^p)  +  C>(ttXandfl.*U  >■ 

**z{(*xPyp)/^l’  ^  d  direct  matter  to  sKow: 


ma»2~c 
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«o  -  (p/N)«p(^*/2niJ(l+(<iO^Xl-^*AnXax{(«j^V>})+0(«*)l  («» 

^l<««:1rflcacq^Mtion(9)ii^all^lilatedfalU>e<^u«l^OB(4XwMl^J|lMrelelll^ngtllecqll^librillmpoCea^il^ 
equatkM  (8),  io  onfer  a,  it  retrieved. 

To  see  what  equeiioo  (10)  offien,  eoorider  the  iteadjrfWe  tmall  rifoal  Wigner  fuactioa  within  a  leiHatioa  time 
aptwoMinatioB,  and  to  second  Older  in  a  (d|, 

«w-‘o-'w((libeto)«A-<*al^I*+(»*'24XVU^  (11) 

te«ri^cj^riMi^]^hyo(llX  the  foBowiog  key  lesalts emerge:  /^*p«/i^*|i,andforj( 

it  -  /4{pax(U+(W)+a3t(W»}  (12) 

which  was  first  obtaiiied  by  Ancona  and  lafrate  (6).  Heic^p*  er/m.  While  this  result  is  ooaatfenf  with  the 
general  phfloio^  of  the  Introdncioiyiintagyapoflf  this  section  the  factor  of  3  on  the  < 

SlieeattorSl^i . —  .  .  ...  ... 

Wigner- 


(see  also 

1  to  Older  a  and  for  a  d 
p^t).  The  WB  equadon  of  modon  ic 

«t*w+(Pa^)aaV(«aU)*p«w+(»*/24X<i»U)«p**wl  “  <w^ 

and  the  first  three  momeiit  equatoo  comspondiiig  to  diat  of  equadoos  (4)  throng  (7)  ate: 

*|P+*x(PdP^)  ■  0 

atO»Pd)+«xO»Pd*/i“)+P«a(U+<yj)+a3r(pkT)  -  appd,con 
a,W +$jfp^w/nii) + a;,(pt]akT/m) +(pp(|^)aj((U  -t-QO) 
-(pa*/l2in)Sx{(a^)/^}ax(P(l/lm)  «•  aW^ 

W  -  3/»kT/2+r»pj*/2ni^a*/24m)ax{(«*(»)/p} 


(M) 

(15) 

(16) 
(17) 


a  <(Pl>d)*  ^  ”  “all  -  (2«»/3dX»x((*i^Va VpI 

<a(pi»d)*  >  ■  0 


(18) 

(19) 


ich  reference  (81  equations  (lOaX  and  (10b)sdiea  combined  with  (lOa)  yield  egoadons  (14)  and  (15)  of 
y,  where  the  colusian  integrals  are  treimfeaeilcalfy.  The  energy  babmeqnttion  is  treated 
'.  Here,  reference  (71  eqiHnan(10bXtdiennBildpled  by  aPd/hiS^added  to  (iDcXtridcfa  is  nailtDled  by 


fiom  which  reference  I 

this  study, '  ..  ... 

Here,  reference  (7]  eoudian  flObX  when  noltipled  by  ap<^  S’added  to  (lOcX  tddefa  is  nailtqtled  by 
l/2m;thecaatmiiityeqiiatioi:istnchideahithisprooedare.  TheRwtofdiisnianhailaDoaitaoae 
dimensional  pbise  space  version  of  the  energy  balance  equation  of  tUs  study,  W  is  T^jiaced  by: 


W-  /.kT/2+r»p4*/2nKrra*/24m)a,{(ajeya} 
From  die  pohtt  of  view  of  device 


(20) 


defined.  Wihh|gE>U4-Q0-fkUii^/pJtriierearl>sn 
can  be  wihtea  at  j  "pa^B  (6),  and  efuttuas  (IS)  arid  (16)  can  be  reeqirwed  as: 

«t(aPd)+fx(<’Pd*^)+P**^  “  *aPd/»ll 

*t^+**0>dW^)+*x(PdPWViBi)+^P<^)*xB 

-(a«*/i2m)ax{(axaya}ax(pdfei)  -  »Wooll 


(12) 


(21) 


(22) 


hi  the  above  form  it  I 
it  introduced  as  a  tia 


tare  governed  by  an  energy  E.  However,E 


Quantum  moment  balance  equations 
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Anticipated  Solutions  of  the  Quantum  Balance  Equations 


For  bound  stttes,J  -  ^  and  equation  (23)  is  an  emenvalue  problem,  one  that  in  the  case  of  a  resonant  tunnel 
structure  leads  to  q^-bound  states.  Further  under  zero  current  conditions,  with  Ep  tepresenting  the 
eigenvalues,  Q  -  ^-U,  and  the  values  of  Q  are  spatially  dependent  and,  in  some  cases,  are  appradmale^ 
e^ial  to  the  bound  nates.  This  result  will  be  prommantfy  displayed  in  the  discusskm  below. 

The  single  particle  Sduodii^  picture  is  lunhed,  in  that  being  dissipationlets  it  does  not  permit  a  direct 
transition  transition  to  a  niu%>tttide  problem  when  oontaeti  are  considered.  For  eiainple,  in^  case  of 
multiparticle  transport  with  electrons  moving balHstically  within  the  N~  region  of  anN'*^ FFTil nmctare,  die 
mean  carrier  energy  and  vetod^  increase  from  the  cathode  to  the  anode.  Conservation  of  mult^wrtklecnnent 
requires  that  increases  in  velodqr  ate  acconmanied  fay  decreases  in  partide  density.  Thusintheabaoaecof 
dissipation  there  will  necessarily  be  ctoge  depletion  at  the  downstream  anode,  unless  disi^iation  Is  present  in 
the  interior  the  device.  If  the  assumption  is  made,  that  the  pltysical  contact  are  boundanes  wherein 
numbers  of  carriers  at  the  cathode  and  anode  are  equal,  then  scattering  within  the  Intarier  of  the  atnettne  la 
eonoeptoally  necessary.  For  the  Itydrodynamic  fonnulation  of  the  sin|pe  partide  Sdgodinger*t  equation,  diete 
is  no  meaning  to  introducing  N***  cathooe  and  anode  r^ons,  since  we  are  dealing  with  a  sin^paitide. 

In  that  dissmation  is  an  essential  feature  of  transport  in  devices,  die  quantum  balance  equadoos  repteaeated  by 
equations  (13)  through  0'^  fonn  «  staring  p<wnt  for  the  simulations  to  be  discussed  behm.  Todat^our 
simulatioiis  indude  tnenrst  two  moment  equations,  and  POissoa’sequarion.  These  have  been  solved  for  a 
qiatially  dependent  effect  mass,  and  Cor  Fermi  stadkies.  ffete,  since  we  have  neither  generated  a  set  of  WB 
moment  etyudons  for  a  spatially  dependent  effective  mass  nor  have  we  obtained  a  di^plaoed  Wigner  fimedon 
that  sadsfiM  Fermi  stadsnes,  we  haw  instead  patdied  on  these  contributions.  Fbrtber,  we  have  treated  die 
factor  of  *3*  associated  with  equations  (14)  ana  (IS)  as  an  adjustable  parameter  that  reflects  the  stadsdcal 
distribudon  used  as  a  basis  for  the  catolation,  as  discussed  m  [6],  and  have  replaced  it  by  unity.  Induscase 
withva  P(|/ni,  the  continuity  equation  is  unchanged,  while  die  momentum  balance  equation  reads:: 


a  ,pv + «x(^v* ) + (^/m){(a  j(U +Q)1 + [(pv*/2)-(NkTF,/,/p)]aj^  nm} + (2/3)a  jNkTF,/,  +pvr  -0  (24) 


Equation  (25)  is  oou|ded  to  the  equation  of  continuity  and  Poisson’s  equation,  with  U(x)  representing  the 
omduedon  band  energy.  The  heterostructure  is  rqnesented  by  the  Anderson  rule:  U  »  z>x(xX  where  x(iO  is 
a  position  dependent  mectronaffloity.  z  is  obtained  from  Foiaon’s  equation:  v<vz  »  e*> -p,1,vriieiec^) 
is  a  position  dqiendentpermitivityiUd^,  is  a  position  dependeiu  doping  level  Fo{oooaictionband 
variations  between  GaAs  and  Aljuaj.,tAs,  the  Kdlowiiigtmtionships  were  used:  m  •>  0.057  0.063X,  nEc - 
a697x. 

Calculations 

The  calculations  discussed  below  ate  for  the  structure  shown  in  figure  1,  with  resonant  tunnelling  barriers 
located  tymmetrically  at  the  center  of  the  structure.  The  structure  and  dimensions  of  the  banien  ate  dimlayed 
in  figure  2,  whidi  shm  the  current  voltage  duracterisda  this  device  at  77K.  There  is  a  weak  r^kn  of 

n^ative  conductance.  The  conhictioa  bind  pn^  at  different  Was  levels,  figure  3,  shows  the  emected  tilt  as 
the  Uks  is  increased.  The  distribation  of  energy  is  sudidiat  at  0.1v,imprcadmatley  20%  of  the  voltage  drop  frlis 
acroH  the  iqistream  accutmikdon  layer,  30%  vritUn  die  confines  of  the  barrier,  and  50%  across  downstream 
fimn  the  second  barrier.  The  dunge  dmribotion,  figure  4,  shows  a  region  of  charge  accumulation  tmstream  of 
the  barrier  that  increases  with  hiaearing  Was,  as  does  the  duuge  hi  the  well  Wmle  different  bounoaiy 
condithnis  have  not  been  studied  these  results  should  be  eatremely  sensitive  to  the  boundary  conditions  at  the 
cathode;  as  should  the  effects  of  inoMporadng  die  energy  balance  equation.  Itisnotdearuiatindudingtbe 
latter  vwllrqirot^  the  diMgedqiletion  at  resonance  seen  Ity  several  other  studies. 

In  all  our  sitnulations,  we  have  noticed  the  formation  a  deletion  layer  downstream  of  the  second  barrier 
once  we  pass  the  vallCT  current  of  the  1«V  characteristic  This  depletion  layer  is  a  ^ledficsin^  particle 
quantum  effect  The  oqiletion  layer  keqis  on  extending  for  biases  greatetthan  the  voltage  at  the  valley  of  the 
I-V  curve  until  the  dqtiraon  layer  toudies  die  heavily  doped  tegioa(210i^*^).  Then,  the  electroa  density 
downstream  of  the  seomd  barrier  gradually  increases  and  the  deletion  region  msappears. 

The  quantum  potential  is  displayed  in  figure  5.  If  we  concentrate  on  its  value  in  the  weU,  the  most  dramatic 
point  to  note  is  t^  as  the  bim  B  raised  m  value  of  the  quantum  potential  tends  to  duster  around  a  narrow 
range,  increasing  in  magnitude  from  the  upstream  barrier  to  the  wiwnstream  barrier.  Within  the  barrier  the 


10.  ELECTRON  AND  HOLE  TRANSPORT 


Multispedes  transport  is  part  of  the  density  matrix  algorithm  with 
^>plications  to  heterostructure  barriers  and  diodes.  While  extensive  woilr  on 
dectron  and  hole  transport  in  underway  as  part  of  the  ULTRA  program,  some  of 
this  type  of  actmty  has  been  published  and  forms  part  of  this  document,  and  is 
incorporated  into  this  section. 
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A.  INTRODUCTION 

When  a  device  physicist  studies  high  field  transport  there  are  several  issues  to  address.  First, 
what  is  the  momentum  space  kinetics,  and  second  what  are  the  high  field  dynamics  in  devices? 
For  much  of  the  early  phases  of  device  studies  in  which  length  scales  were  of  Um  order  of  tens 
and  hundreds  of  microns,  a  satisfying  picture  often  emerged  in  which  it  was  assumed  that  the 
^ce  charge  distribution  did  not  affect,  in  any  significant  way,  the  kinetics  of  transport  In  this 
picture  the  semiconductor  was  represented  as  a  material  with  a  specified  field  depenttent  velocity, 
which  in  the  case  of  GaAs  sustained  a  region  of  negative  differential  mobility  associated  widi  k- 
s^ce  transfer.  As  the  structure  lengA  scales  decreased,  it  was  clear  that  the  separation  of  the 
l^etics  from  the  space  charge  contributions  was  no  longer  possible,  and  a  more  complicated 
pictim  emerged  that  ^uiied  the  use  of  advanced  numerical  algorithms  for  the  study  of  device 
physics.  These  algorithms,  however,  depended  in  a  detailed  way  on  uncertain  parameters  diat 
were  used  in  the  k-space  calculation.  This  was  the  best  we  could  do;  and  today  when  hi^  field 
transport  is  considet^  it  is  most  often  examined  without  reflect  to  device  configuration.  In  the 
discussion  below,  we  take  this  viewpoint 

>^th  high  field  transport  considered  within  a  device  context  the  ternary  alloy  AlGaAs  can  be 
regarded  as  an  enabling  material.  For  while  AlGaAs  for  a  range  of  aluminium  mole  fraction 
possesses  a  region  of  negative  differently  mobility  (NDM),  other  samples  possess  superior  NDM 
r^ons.  A  major  interest  in  AlGaAs  lies  in  the  fact  that  within  the  And^on  rule,  its  dectron 
amnity  is  significantly  different  ftem  such  materials  as  gallium  arsenide  and  indium  gallium 
ars^de.  This,  in  ite  sinqilest  version,  was  responsible  for  die  presence  of  barrier  structure 
devices,  and  the  earliest  AlGaAs/GaAs  devices  were  among  the  fi^  Ixmd  structure  engineered 
devices.  While  a  host  of  AlGaAs/GaAs  devices  have  emerged  as  a  result  of  barrier  engineering 
we  will  focu5  attention  on  only  several  v^ions  of  this  enabling  technology.  The  structures  we 
will  focus  on  are  Gunn  diodes,  AlGaAs/GaAs  HBTs,  AlGaAs/GaAs  MODFETs  and 
AlGaAs/GaAs  BICFETs.  The  emphasis  will  be  on  transport  The  following  will  brirfly  review 
the  high  field  k-space  transport  properties  of  AlGaAs,  and  then  turn  to  high  field  transport  in  the 
devices  mentioned  above. 


B .  HIGH  FIELD  BULK  PROPERTIES 

Electron  transit  on  the  AlGaAs  alloy  depends  in  a  detailed  maimer  on  the  nimbers  of  carriers  in 
die  r,  L  and  X  portions  of  the  conduction  band.  For  moderate  values  of  the  iiiole  composition  of 
aluminium,  and  rt  low  values  of  electric  field  the  electrons  are  dominandy  in  the  F  valley.  As  the 
field  inoeases  and  LO  phonon  scattering  no  longer  effectively  removes  ^  excess  carrier  enngy, 
a  CCTtain  fraction  of  electrons,  with  the  assistance  of  phonons,  transfers  to  the  subsidiary  L  vriley, 
of  which  th«'e  are  four  ^uivalent  valleys.  The  rate  at  which  these  electrons  are  transferred 
determines  whether  negative  diffnential  mobility  wll  occur.  At  further  increases  in  field  the 
electrons  can  also  transfer  to  the  next  higho*  valley,  the  X  valley,  of  which  diere  are  three 
equivalent  valleys.  Transfer  between  any  two  valleys  including  equivalent  valleys  occurs.  The 
specific  properties  needed  to  determine  these  transfer  rates  are  listed  in  TABLE  1,  where  the 
denrity  of  states  effective  mass  in  the  table  is  m^  s  (mi^mpe,pe„dicai„)’^ .  / 
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TABLE  1  Properties  of  AljGaj., As 


Patameter  I  GaAs  AlAs  |  AljOai.^As 


Static  dielectric  constant  e. 

13.18 

10.06 

13,18 -3.12X 

High  frequency  dielectric  constant 

10.89 

8.16 

10.89  -  2.73x 

1 

Band  gap  energy  Eg  (eV) 

1.424 

2.168 

1.424  +  1.247x(0<x<0.45) 

4.07  -  l.lx  (0<x<0.45) 

Indirect  band  gap  (L) 

1.708 

2.35 

1.708  -  0.642X 

Indirect  band  gap  (X) 

1.900 

2.168 

1.900  +  0.125X +  0.143x2 

Electron  affinity  Xe  (eV) 

Density  of  states  electron  mass  nv 

4.07 

3.5 

3.64-0.14x(0<x<0.45) 

1.900  +  0.125X  +  0.143x2  (0.45<x) 

r  valley 

0.067 

0.150 

0.067 +  0.083X 

X  valley 

0.49 

0.34 

0.49-0.15X 

L  valley 

1 

0.22 

0.26 

0.22  +  0.04X 

Acoustic  deformation  potential  E  (eV) 

6.7 

5.5 

6.7  -  1.2x 

Phenomenological  def.  potential  Em  (eV) 
Intervalley  def.  potential  D(ij)  (eV/cm) 

Dcrpo 

DCTJL) 

D(XJL) 

3.6 

{0.5-l.l)xl0’ 
(0.15- 1.0)  xio’ 
(034-l.l)xl0’ 

2.9 

3.6-0.7X 

D(X,X) 

D(UL) 

(0.27- 1,1)  xio’ 
IxlO’ 

1.47  xio’ 

LO  phonon  energy  (eV) 

Intervalley  phonon  ener^  (eV) 

0.033 

0.050 

0.033 +  0.017X 

{r=>x) 

0.0299 

0.0299  +  0.0175X 

(n=>L) 

0.03 

0.03  +  0.0134X 

(X=>L) 

0.0293 

0.0293 +  0.0181X 

(X=>X) 

0.0299 

0.0299  +  0.0175X 

(I^L) 

0.029 

0.029 +  0.0150X 

A  number  of  studies  have  been  pofonned  which  show  the  degree  to  which  the  percentage  of 
aluminium  affects  the  region  of  NDM.  The  parameters  used  in  these  studies  sustain  a  certain 
degree  of  uncertainty  in  that  such  scattering  contributions  as  the  deformation  potential 
contribution,  the  relevant  optical  phonon  frequencies,  etc.,  are  all  dependent  upon  die  aluminium 
mole  fraction.  A  representative  c^culation  was  performed  in  1988  [1].  This  was  a  Monte  Cario 
calculation  of  the  field  dependent  velocity  of  AIq  32Gao  ^gAs.  Comparison  was  not  made  to 
experiment,  as  no  experimental  drift  velocity  m^urements  were  available.  However,  the 
calculation  is  of  sufficiently  general  nature  to  reproduce  its  results  here.  In  ref  [1]  a  set  of 
parameters  was  chosen  from  Ae  literature,  about  which  a  parameter  variation  was  made.  The 
parameters  are  ^proximately  represented  by  inserting  32%  into  the  expressions  of  TABLE  1, 
from  which  the  velocity  field  curve  of  FIGURE  1  was  obtained. 
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FIGURE  1  Field  dependent  velocity  for  AIq  32630  ^gAs,  from  [1]. 

Several  points  are  relevant:  first  the  peak  velocity  is  slightly  in  excess  of  9  x  10^  cm/s,  which  is 
considerably  below  that  of  gallium  arsenide.  The  saturated  drift  velocity  is  however 
approximately  equal  to  that  of  gallium  arsenide.  The  study  in  [1]  undertook  a  number  of 
parameter  variations.  For  example  when  the  F-L  separation  was  reduced  there  was  as  expected  a 
largo-  fraction  of  carriers  in  the  L  valley  in  equilibrium,  with  a  consequent  reduction  in  the 
mobility.  Indeed  for  an  energy  separation  of  91  meV  the  NDM  region  is  absent  For  an  increase 
in  the  F-L  separation  to  140  meV  there  is  a  reduced  population  in  the  L  valley,  the  mobility  is 
higher  and  a  region  of  negative  differential  mobility  is  present  with  an  increase  in  the  peak  velocity 
to  nearly  1.2  x  10^  cm/s.  Increasing  the  effective  masses  of  all  the  valleys  tends  to  reduce  the 
peak  velocity  as  well  as  the  saturated  drift  velocity.  The  LX)  phonon  scattering  rate  increases  as  the 
optical  phonon  frequency  increases.  Increase  in  the  scattering  rate  results  in  fewer  electrons 
h^ted  to  sufficient  values  for  intervalley  transfer  to  occur.  It  was  found  that  for  an  optical  phonon 
energy  increase  from  38  to  45  meV  the  peak  velocity  increased  by  approximately  15%.  A 
decrease  in  the  optical  phonon  energy  to  30  meV  also  resulted  in  a  near  15%  decrease  in  peak 
velocity.  It  is  important  to  note  that  for  percentages  of  aluminium  greato-  than  50%,  the  region  of 
NDM  disappears  (see  the  general  discussion  in  [2]).  Additionally,  in  the  vicinity  of  40% 
aluminium  Ae  band  structure  eno-gy  minimum  order  changes  from  Ae  F-L-X  ordering  to  the 
X-L-F  ordering. 

The  phrase  heated,  in  the  above  paragraphs,  implies  nonequilibrium  electrons,  and  an  electron 
tempmture  model  is  the  one  most  often  invoked  to  deal  with  this  description.  Th^  are  a  variety 
of  means  by  which  the  electron  temperature  model  is  invoked,  the  most  common  being  usually 
predicated  on  a  displaced  Maxwellian  for  a  distribution  function  and  solving  three  sets  of 
equations:  a  carrier  balance  equation,  a  momentum  balance  equation  and  an  energy  balance 
equation.  This  model  is  briefly  considered  as  the  language  associated  with  it  is  invoked  in  the 
device  discussion. 

Under  uniform  field  conditions  these  equations  for  two  levels  of  transfer  (e.g.  F  and  L)  rqircsent 
particle,  momentum  and  energy  conservation.  For  particle  conservation: 
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(la) 

(lb) 


da|n,^t  =  -a,n,r,  +  ajOjFj 

3a2n2^^t  “  a|n|F|  ■  a2n2F2 

For  Eqn  (la)  there  are  n,  electrons  in  each  lower  energy  valley  and  n2  electrons  in  each  higher 
energy  valley.  Eqn  (la)  indicates  that  there  are  aiOf  electrons  scattered  out,  distributed  equally  to 
the  a  2  higher  energy  va^eys;  and  there  are  a2n2  electrons  scattered  from  the  higher  energy  valleys 
into  die  lower  energy  valleys.  The  respective  cairi^  scattering  rates  arc  designated  Fi  and  F2. 
For  uniform  fields  and  steady  state,  a  condition  under  which  the  velocity  field  curve  is  generated, 
the  following  condition  holds;  ajniFi  =  a2n2F2. 

For  momentum  conservation,  the  second  set  of  equations  describes  the  rate  of  change  of 
momentum  (or  velocity)  in  the  individual  valleys  under  application  of  an  implied  Held  and 
scattering  events.  Under  uniform  fields  and  for  the  low  energy  carriers  this  equation  is: 

8(n,pi)/9t  =  -n,cF„  -  n,p,F3  (2) 

where  the  momentum  is  designated  Pj,  s  m,Vt,  and  F3  is  the  momentum  scattering  rate  for  the 
low  energy  carriers.  Under  steady  state  conditions:  nimjV]  =  -(l/F3)n]eF„.  Similar  equations 
can  be  written  for  the  high  energy  valley  carders  and  for  holes. 

For  energy  conservation,  there  are  various  forms  in  which  the  lower  and  higher  energy  valley 
energy  equations  can  be  described.  We  cast  the  energy  equations  in  terms  of  the  species  *  1  ’  and 
species  ‘2’  electron  temperatures: 

3(ajn|Tj)^t  =  (mj VjV(3kg)]  {a|nj(2F3-F|)  +  ajn2r2}  -  ainiTjF5  +  a2n2T2F^  (3) 

In  Eqn  (3)  kg  denotes  tiie  Boltzmann  constant;  r5  denotes  energy  relaxation  within  die  species  *  1  * 
valley  plus  eneigy  exchange  with  the  species  *2’  valley;  Fg  denotes  return  energy  between  qiecies 
*2*  and  species  *r  valleys.  The  above  analysis  requires  calculations  of  die  scattering  rates.  These 
are  taken  from  scattering  integ^s.  For  a  review  see  [2],  in  which  the  sigriificance  of  the  above 
description  is  dwelled  upon.  In  particular,  with  the  electric  field  as  a  ^ving  force  the  three 
parameters  density,  momentum 
and  temperature  are  determined, 
as  a  function  of  field.  It  is 
generally  assumed  that  under 
equilibrium  conditions  the 
distribution  of  carriers  in  each 
of  the  valleys  is  determined  by 
the  density  of  states  of  each 
vall^  and  the  energy  sqiaration 
of  each  valley. 

Adachi  [3]  has  performed  some 
of  the  above  electron 
temperature  calculations.  In 
particular,  an  estimate  of  the 
inoease  of  electron  ten^peratuie 
as  a  function  of  electric  field  for 
polar  phonon  scattering  is 
shown  in  FIGURE  2.  It  is 
clear  from  FIGURE  2  that  as 
the  alloy  composition  is 
increased  carrier  heating  is  less 
severe,  at  any  given  value  of 


vf 


ELECTRIC  RELD  E  ( V/cm  ) 

FIGURE  2  Electron  temperature  as  a  function  of  electric  field  for 
various  alloy  compositions,  under  the  condition  of  polar  opfical 
scattering.  From  [3]. 
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field.  When  all  of  the  scattering 
mechanisms  are  included  the 
electron  temperature  increases 
much  more  rapidly  with  field, 
indicating  that  at  high  values  of 
electric  field  there  is  a  considerable 
amount  of  transfer  from  the  low 
energy  regions  to  the  high  energy 
regions  [2]. 

Polar  optical  scattering  also 
dominates  the  low  field  mobility  of 
many  of  these  alloys,  and  in 
FIGURE  3  the  dependence  of  the 
polar  optical  scattering  limited 
mobility  for  GaAs  and  indicated 
alloys  is  shown,  with  the  apparent 
degradation  as  the  alloy 
composition  is  increased. 


C.  PERPENDICULAR  TRANSPORT  -  HOT  ELECTRON  INJECTOR 
CATHODES 

While  most  of  the  discussion  of  high  field  transport  in  AlGaAs  will  be  concerned  with  transport 
parallel  to  the  interface,  the  following  discussion  on  perpendicular  transport  will  provide  an 
indication  of  how  high  fidd  transport  in  devices  is  affect^  by  such  materials  as  AlGaAs. 

Consider  FIGURE  4,  which  within  the  context  of  the  Anderson  model  is  a  generic  conduction 
band  diagram  of  a  thrw  le^on  stractuie:  HEA-LEA-HEA  (H:  hi^,  L:  low,  EA:  electron  affinity). 
In  die  case  of  hot  electron  injector  cathodes  [4],  a  recent  application  choice  for  Gunn  diodes,  die 
HEA  material  was  GaAs  while  the  LEA  material  was  AlGaAs.  As  configured  the  structure 
consisted  of  an  n^  heavily  doped  cathode,  a  lineariy  graded  AlGaAs  region  extending  ower  SOO  A 
with  a  height  of  300  meV,  an  n^  spiked  doping  layer,  tqiproximately  50  A  long,  followed  by  a 
drift  region  and  an  n^  anode.  The  design  parameters  of  the  structure  are  consistent  with  die 
requirement  that  for  Gunn  oscillations  to  occur  electrons  must  transfer  firom  a  low  energy  region 
to  a  higher  energy  with  a  consequent  local  region  of  negative  differential  mobility.  Early  theories 
recognised  that  satisfactory  control  of  the  oscillations  n^uired  the  presence  of  local  regions  to 
force  at  least  one  parameter,  e.g.  a  bounding  electric  field,  to  be  insensitive  to  bias  and 
temperature  conditions  such  that  electrons  wodd  enter  the  drift  region  with  a  distribution  of 
energies  consistent  with  a  sufficient  number  of  carriers  in  the  L  valley,  in  the  case  of  e.g.  GaAs. 

The  equilibrium  band  structure  for  this  case  (ignoring  subsidiary  valleys  and  inyddng  Boltzmann 
statistics)  is  represoited  in  FIGURE  4  for  an  undoped  AlGaAs  launcher  One  Concqpt  bdiind  this 
design  is  that  carriers  enter  the  drift  region  with  a  non-zero  velocity  whose  value  is  estimated  from 
the  conversion  of  carrier  potential  energy  to  kinetic  energy.  It  is  assumed  that  a  large  fraction  of 
the  electrons  that  enter  the  device  are  F  valley  carriers,  as  are  those  that  are  in  the  AlGaAs  region. 
Those  carriers  that  pass  through  the  spiked  region  are  F  valley  carriers,  but  when  they  get  to  the 
boundary  of  the  drift  region  thoe  is  enough  energy  to  place  a  considerable  number  of  carriers  in 
the  L  valleys.  The  design  appears  to  provide  improved  performance  of  the  transferred  electron 
oscillator,  but  the  presence  of  the  wide  band  gap  AlGaAs  next  to  the  GaAs  introduces,  as  in  all 
structures  of  this  type,  a  triangular  potential  well  and  a  force  that  tends  to  confine  carriers.  Indeed 
the  presence  of  such  a  force  i  carriers  in  the  continuum  is  of  the  order  of  10  kV/cm  and  in  the 
wrong  direction!  Under  bias  of  course  there  may,  depending  on  conditions,  be  a  net  force  acting 
on  the  carriers  pulling  them  into  the  transit  region,  but  then  the  operating  conditions  would  display 


FIGURE  3  Polar  optical  scattering  limited  mobility  as  a 
function  of  alloy  composition  using  the  electron  temperature 
model.  From  [31. 


-201- 


7J  High  field  transport  and  the  Gunn  effect  in  AIGaAs/GaAs  structures 

a  significant  dependence  on  bias.  The  origin  of  this  dilemma  lies  in  the  incomplete  manner  in 
which  transport  problems  are  addressed.  This  incomplete  feature  is  highlighted  here  because  the 
contribution  we  are  about  to  focus  on  is  present  in  all  structures  where  there  are  large  carrier 
density  gradients,  often  associated  with  heterojunction  interfaces. 


FIGURE  4  Equilibrium  density  and  potential  energy  for  a  hot  electron 
injection  launcher  with  an  n*^  spike  of  lO’^/cm^.  Rom  [S]. 

The  calculation  displayed  in  FIGURE  4  was  obtained  through  a  solution  to  the  Liouville  equation 
in  the  coordinate  representation  [S],  and  contains  the  relevant  quantum  features.  However, 
through  an  analytical  expansion  of  the  Liouville  equation,  in  which  classical  tran^rt  dominates 
and  quantum  contributions  are  treated  as  corrections,  and  Boltzmann  statistics  prevail,  it  has  been 
demonstrated  that  the  net  driving  force  on  carriers  is  given  1^: 

F  =  -V(V  +  Q/3)  (4) 

Q  =  -(h2/2m)(p-'/2)(d2  p«/2/92  >.)  (5) 

We  have  calculated  Q,  often  referred  to  as  the  quantum  potential,  from  the  structure  of  FIGURE  1, 
and  find  that  in  equilibrium  thore  is  a  net  force  acting  on  the  carriers  that  is  of  the  order  of 
10  keV/cm,  acting  in  such  a  direction  as  to  move  carriers  into  the  drift  region.  It  is  perhaps 
important  to  emphasise  that  the  configuration  of  FIGURE  4  is  representative  of  p^pendicular 
transport;  and  for  the  specific  situation  of  the  hot  electron  launcher  it  is  expected  that  die  details  of 
such  things  as  the  field  dependent  velocity  may  be  of  secondary  importtuice  to  the  feature  of 
providing  a  heterostructure  offset  region. 


D.  PERPENDICULAR  TRANSPORT  -  AIGaAs/GaAs  HETEROSTRUCTURE 
BIPOLAR  TRANSISTORS 

The  second  example  of  perpendicular  transport  is  that  of  the  heterostructure  bipolar  transistor 
(HBT).  In  this  case  the  band  structure  of  the  device  falls  into  the  same  generic  category  as  that  of 
Ae  hot  electron  injector.  Typically  the  HEA  region  is  a  heavily  doped  n*  gallium  arsenide  region. 


-202- 


High  fittd  transport  and  the  Gunn  effect  in  AlGaAsiGoAs  structures 


followed  by  an  LEA  region  which  is  a  heavily  do]^  n*  AlGaAs  region  (whose  properties  are 
often  igpplication  speciEc),  followed  by  HEA  material.  The  HEA  material  is  often  a  standard  p^ 
GaAs  b^  followed  by  a  low  doped  GaAs  collector  region.  Most  of  the  attention  associated  with 
this  structure  involves  compositional  grading  of  the  AlGaAs  emitter.  Here  there  are  several 
aspects  to  consider.  First,  the  electrons  must  get  from  the  heavily  doped  GaAs  region  to  the  wide 
band  gap  AlGaAs.  If  the  first  n*(HEA)n*(LEA)  interface  is  not  obliterated  in  the  device 
processing  steps,  then  the  interface  is  expected  to  look  like  that  shown  in  FIGURE  S,  where  the 
wide  band  gap  material  is  at  the  left.  Note  that  in  equilibrium  the  density  in  the  wide  band  gap 
material  approaches  background,  which  for  this  case  is  lO'Vcm^,  within  3(X)  A,  and  is  relatively 
insensitive  to  the  density  of  the  adjacent  material.  Similar  remarks  apply  to  the  potential  energy, 
which  is  also  shown. 


-«00  -300  0  300  coo  -COO  -300  O  300  COO 

OWone*  (AngcfracM)  Along  Th«  Otngoool  OMnneo  (AngolranM)  Along  Iho  Otagonol 


FIGURES  Quantum  mechanical  calculation  of  die  distribution  of  charge  and  potential  energy 
for  a  wide  band  gap/harrow  band  gap  stnicture  with  varying  doping  distributions.  From  [5], 

In  the  design  of  the  HBT,  if  the  upstream  emioer  interface  is  graded,  as  in  die  hot  electron  injector, 
the  necessity  of  requiring  tunnelling  mechanisms  to  move  die  electrons  from  the  gallium  arsenide 
region  to  the  aluimnium  gallium  arsenide  region  is  minimised.  But  most  of  die  effort  in  designing 
the  HBT  is  concerned  with  the  compositional  grading  of  the  AlGaAs  in  die  vicinity  of  the  ba^ 

The  quantum  mechanical  calculation  of  the  equilibrium  potential  energy  profile  and  electron  and 
hole  density,  for  an  emitter  that  incoiporates  an  abrupt  1000  A  long,  300  meV  barrier  adjacent  to 
the  base,  is  displayed  schematically  in  FIGURES  6  and  7,  respectively.  For  this  calculation  an 
acceptor  doping  of  lO'^/cm^,  and  a  wide  band  gap  emitter  doping  of  lO'^Airf,  was  assumed. 
This  latter  is  generally  an  order  of  magnitude  below  that  of  the  usual  design  of  the  HBT,  but  is 
sufficient  to  illustrate  the  features  of  the  role  of  the  het^structure  in  the  device  -  it  prevents  the 
diffusion  of  holes  from  the  base  to  the  enutter  while  enhancing  the  injection  of  electrons  into  the 
base.  The  latter  is  represented  by  die  dip  in  the  potential  at  die  n-hetero-p  int«face.  Thoe  is  a 
diffusion  of  mobile  holes  from  the  base  to  the  collector  re^on,  that  follows  the  standard  results 
when  recombination  is  ignored.  Hi^er  doping  in  the  emitter  v^l  pull  down  the  potential  within 
the  centre  of  the  heterostructure  region  to  a  near  zero  value,  which  is  consistent  with  the  higher 
emitter  background  doping.  While  quantum  effects  associated  ^di  electron  injection  into  the  base 
are  apparent,  on  the  emitter  side  the  hole  density  goes  from  its  peak  value  within  the  base  to  a 
negligible  value  at  a  distance  of  approximately  70  A  into  the  emitter  base;  quantum  effects  are 
likely  here.  The  details  indicate  that  the  band  structure  of  the  base  is  dominated  by  near  charge 
neutrality  within  the  base  away  ffom  the  heterointerface.  Notice  that  the  decrease  of  £e 
conduction  band  energy  on  the  emitter  side  of  the  barrier  corresponds  to  the  increase  in  the  space 
charge  on  the  knitter  side  of  the  barrier.  This  is  similar  to  the  space  charge  distribution  in 
FIGURE  5. 
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FIGURE  6  Quantum  mechanical  calculation  of  the  equilibrium 
conduction  and  valence  band  energies  for  an  HBT  with  an  jdxupt 
electron  barrier  in  the  emittec  From  [S]. 
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FIGURE  7  For  the  calculation  of  FIGURE  6,  the  eqi»’' '  ^^'im 
electron  (solid)  and  hole  (dashed)  distribution  for  tf 
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E.  PERPENDICULAR  TRANSPORT  -  THE  BIPOLAR  INVERSION  FIELD 
EFFECT  TRANSISTOR 

The  bipolar  inversion  field  effect  transistor  [6]  (BICFEp  has  the  same  generic  ordering  as  that  of 
the  two  previous  devices,  but  is  of  interest  because  in  some  of  its  present  configurations  it 
involves  the  strategic  placement  of  a  planar  doped  (or  delta  doped)  layer  of  acceptors  or  donors. 
In  particular,  the  configuration  of  the  HBT  shown  in  the  previous  figures  can  be  altered  to  that  of 
the  BICFET  by  replacing  the  heavily  dop^  p  base  with  a  planar  doped  p  region,  approximately 
30  A  long.  Unlike  ordinary  inversion  regions  which  arise  at  suitable  values  of  the  local  potentisil 
energy,  tlw  local  region  of  holes  generated  by  planar  doping  is  drought  to  contain  many  of  the  bias 
dqrendent  characteristics  of  the  more  common  inversion  layer  (hence  the  term  inversion  in  the 
name  of  the  device). 

For  the  BICFET  the  operating  voltages  are  different  because  the  critical  device  leng^s  are  of  a 
nanostructure  scale.  For  the  configuration  of  an  HBT,  and  on  the  basis  of  earlier  studies  [7],  it  is 
anticipated  that  the  presence  of  the  planar  doped  layer  would  lead  to  potential  contours  that  would 
lie  parallel  to  the  interface  everywhere  except  at  the  vicinity  of  the  metallisation  or  contact  regions. 
In  the  vicinity  of  the  contacts  they  would  spread  fix>m  a  small  region  at  the  emitter  planar  doped 
edge  into  the  contact  regions.  For  the  case  in  which  the  planar  doped  barrier  is  introduced  as  a 
replacement  for  the  base  in  the  HBT.  it  is  anticipated  that  such  terms  as  the  base  transit  time  would 
improve,  simply  because  of  a  reduction  in  the  base  dimensions.  The  presence  of  the  wide  band 
gtqp  matoial  such  as  aluminium  gallium  arsenide  is  crucial  for  die  operation  of  the  device  because 
it  eliminates  the  possibility  of  a  remote  migration  of  holes  toward  the  emitter  and  confines  them  to 
the  collector  region. 


F.  PARALLEL  TRANSPORT  -  GaAs/AIGaAs  MODFETS 


While  much  recent  activity  has  concentrated  on  tite  structures  discussed  above  paiticularly  with 
respect  to  the  analog  and  ^gital  propeities  of  die  devices,  the  vertical  devices,  as  the  above  are 
referred  to,  do  not  reflect  the  transj^rt  properties  of  the  ternary  AlGaAs.  Rather  they  reflect 
primarily  ^  band  structure  of  the  material  in  concert  with  the  lower  band  gap  matoial.  To 
examine  the  role  of  high  field  transport  in  devices  and  the  Gunn  effect,  we  ne^  to  examine  the 
MODFCT  as  a  generic  device  (see  ref  [8]  for  a  review). 


HGURE  8  is  a  sketch  of  a 
MODFET,  and  while  these 
structures  can  be  very 
complicated,  the  essential 
feature  of  the  device  is  that  the 
wide  band  gap  material  is 
doped,  and  a  two  dimensional 
electron  gas  forms  within  the 
narrow  band  gap  material,  near 
the  heterostructure  interface. 
The  structure  shown 
incorporates  AlGaAs/GaAs. 
Other  structures  incorporate 
AlGaAs/InGaAs/GaAs,  which 
includes  the  possibility  of 
transport  in  a  quantum  well. 

Since  the  offset  voltage  of  the 
AlGaAs/GaAs  system  is 
dependent  upwi  the  mole 
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RGURE  8  Schematic  of  an  n-channel  MODFET. 
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fraction  of  aluminium,  and  since  the  distribution  and  number  density  of  the  two  dimensional 
'  electron  gas  are  dependent  upon  the  barrier  height  between  the  dissimilar  materials,  as  shown  in 
nOURE  9,  the  mole  fraction  becomes  an  important  design  feature.  For  the  FIGURE  9 
calculations  it  is  noted  that  the  wide  band  gap  material  is  uniformly  doped,  and  that  the  entire 
offset  is  at  the  heterostructure  interface.  The  density  distribution  displays  a  decrease  within  the 
wide  band  gap  material  where  a  minimum  is  reached.  The  maximum  value  of  charge  density 
occurs  to  the  left  of  the  interface  and  within  the  narrow  band  gap  portion  of  the  structure.  For  the 
100  meV  offset  calculation  the  peak  density  approaches  7  x  10*  vem^,  for  an  approxinute  sheet 
carrier  concentration  of  3  x  10  "/cm^.  For  the  200  meV  offset  calculation  the  sheet  carrier 
concoitration  increases  by  approximately  40%. 
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FIGURE  9a  Quantum  mechanical  calculation  of  the  dependence  of  two  dimensional 
electron  gas  on  the  offset  voltage.  Potential  and  electron  distribution  for  an 
n'ht'(LEA)n'^(HEA)  stnictiue  with  an  r^set  of  100  meV.  From  (S). 
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RGURE  9b  As  in  FIGURE  9a,  but  for  an  offset  of  200  meV. 
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An  examination  of  FIGURE  9  provides  some  important  features  of  transport  in  the  MODF^'. 
First  we  are  interested  in  transport  in  the  GaAs  region.  In  this  region,  there  is  a  high 
concentration  of  electrons  whose  origins  are  in  the  wide  band  gap  material.  Thus  id^ly,  there  is 
no  impurity  scattering  or  alloy  scattering  to  worry  about.  But  if  we  look  at  the  distribution  of 
charge  in  the  aluminium  gallium  arsenide,  it  is  seen  that  only  a  marginal  amount  of  current  flows 
in  the  structure  near  the  interface,  but  a  substantial  amount  of  current  can  flow  at  regions  away 
from  the  interface.  Using  the  parameters  discussed  earlier  it  is  known  that  the  alloy  and  density 
contribution  and  the  scattering  contributions  in  the  aluminium  gallium  arsenide  indicate  that  any 
contributions  to  the  cuircnt  from  the  wide  band  gap  material  are  undesirable.  This  means  that 
design  efforts  must  be  introduced  to  minimise  the  contributions  of  transport  within  the  AlGaAs 
region.  One  prominent 
means  is  to  introduce  a 
region  of  planar 
doping  into  the  wide 
band  gap  region.  The 
situation  is  displayed 
in  FIGURE  10  for 
an  undoped  400 
meV  barrier,  that 
incorporates  a  30  A 
wide  planar  doped 
region  with  a  doping 
of  10‘Vcm3.  The 
important  feature  to 
note  is  that  there  is 
very  little  mobile 
charge  within  the  wide 
band  gap  material,  and 
thus  the  parasitic 
current  is  minimised. 

The  details  indicate  that 
the  petdt  in  charge  falls 
within  the  nairow  band 
gap  material,  and  that 
the  presence  of 
exposed  donors, 
associated  with  the 
planar  doping  within 
the  barrier,  results  in  a 
potential  energy 
minimum  in  the 
vicinity  of  the  donors. 

But  the  significant 
reduction  in  charge 
density  implies  a 
reduction  in  parasitic 
current 


In  addition  to  the 
above  issues  several 
odiers  emerge.  For  the 
materials  of  interest 
there  are  three  bands  to 
consider,  the  F,  L  and 
X  valleys  for  each  of 


RGURE  10  Quantum  mechanical  calculation  of  the  dependence  of  two 
dimensional  electron  gas  fonn  in  the  presence  of  a  planar  doped  donor  region. 
Potential  and  electron  distribution  for  an  n'^(HEA)n'(IJEA)n'(HEA)  structure 
with  an  offset  of  400  meV. 
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the  materials.  Here,  generally  only  the  offset  voltage  between  the  gamma  valleys  of  the  two 
materials  are  consider^  For  concreteness  consider  AE(x,i  J)  where  ij  =  (r,L,X)  and  x  denotes 
the  composition  of  aluminium.  Then  while  the  AE(x,r,  O  increases  as  x  increases  from  x  =  0,  for 
sufficiently  large  x  the  subsidiary  valleys  begin  to  influence  the  statistics.  The  commonly  used 
value  for  x  in  MODFETs  was  30%,  where  the  ofrset  voltage  is  approximately  260  meV.  At  this 
value  the  commonly  known  DX  centre  comes  into  play  and  leads  to  instabilities  particularly  in 
digital  circuits,  as  a  result  of  which  much  effort  has  been  confined  to  MODFETs  with  lower 
aluminium  composition. 


G.  PARALLEL  NONEQUILIBRIUM  TRANSPORT  -  HIGH  FIELD  EFFECTS 
AND  MODFETS 

We  now  examine  high  field  transport  in  these  structures.  There  are  several  ways  in  which  these 
problems  can  be  stuped.  The  drift  and  diffusion  equations  can  be  invoked,  requiring  the  presence 
of  field  dependent  drift  velocity  for  the  carriers,  a  nonequililMium  formulation  transport  can  be 
examined.  For  the  case  of  nonequilibrium  transport,  the  conunon  approach  has  been  to  invoke 
either  Monte  Carlo  or  balance  equation  procedures.  Only  the  Monte  Carlo  approach  will  be 
discussed  below. 

A  particularly  interesting  set  of  results  has  been  discussed  in  [9]  and  [10].  Both  of  these 
calculations  include  real-space  transfer  [11],  permitting  electrons  from  the  narrow  band  gap 
material,  with  a  sufficiendy  large  energy,  to  transfo'  back  into  die  wide  band  gap  material.  Typical 
material  parameters  used  in  these  studies  are  represented  in  TABLE  1. 

The  structure  studied  in  [9]  consisted  of  a  souixre  drain  spacing  of  0.75  microns  and  a  gate 
contact,  0.5  microns  long,  located  0.1  microns  from  the  source  boundary.  The  diickness  of  fhe 
AlGaAs  was  400  A  with  a  uniform  donor  coiurentration  of  10 ‘‘Am’.  The  composition  of 
aluminium  was  22%  with  a  barria  near  190  meV.  The  gallium  arsenide  layer  was  u^oped  and 
0.2  microns  thick.  While  calculations  were  performed  to  examine  the  switching  speed  in  going 
from  one  voltage  state  to  a  state  in  which  the  gate  and  drain  voltages  were  respectively  0.5  V  and 
1.0  V,  we  will  concern  ourselves  only  widi  the  steady  state. 

What  is  to  be  expected?  For  these  calculations  in  which  only  one  volt  falls  across  the  source  to 
drain  region,  the  average  electric  Held  is  approximately  15  kV/cm,  and  is  nonuniformly  distributed 
and  determined  by  Poisson's  equation.  While  the  electrons  acquire  energy  from  the  field  as  they 
traverse  the  structure  there  are  phonon  losses  within  the  gallium  arsenide  layer,  and  the  average 
energy  of  the  carriers  increases  to  only  350  meV  at  the  end  of  the  channel.  Witiiin  tite  channel  ai^ 
for  approximately  75%  of  the  channel  length  most  of  the  carriers  in  the  GaAs  are  F  valley  carriers, 
with  most  of  the  transfer  to  the  subsidiary  valleys  occurring  near  the  downstream  portion  of  tite 
structure  as  shown  in  FIGURE  11a.  The  interesting  feature  of  this  result  is  that  most  of  the 
carriers  in  the  structure  are  gamma  valley  carriers,  then  the  electron  velocity  within  the  chatmel 
should  be  dominated  by  the  low  field  mobility  of  the  materiaL  But  a  study  of  tiie  mean  carrier 
veloci^  (FIGURE  1  lb)  in  the  structure  demonstrates  that  in  regions  where  there  is  sig^iificant 
electron  transfer,  and  a  reduction  of  the  numbers  of  gamma  valley  carriers,  the  avoage  electron 
velocity  in  the  gallium  arsenide  continues  to  increase.  This  is  reasonable,  since  the  n»^ty  of  the 
subsidiary  valleys  is  considerably  smallo-  than  tiiat  of  the  gamma  valley  and  thereby  makes  a 
negligible  contribution  to  the  total  current  Thus,  even  when  electron  transfo-  occurs  the  current  is 
dominated  by  the  high  mobility  carriers. 
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(a)  (b) 

FIGURE  1 1  (a)  Average  election  eneigy  and  peiceniage  of  F-valley  elections 
along  the  GaAs  channel  for  y  =  0.3.  (b)  Average  election  velo^.  [9] 

Another  feature  of  importance  is  the  relative  concentration  of  carriers  in  die  subadiaty  valleys  of 
AlGaAs.  According  to  the  model  of  [9],  widi  die  exception  of  a  composition  of  10%  aluminium, 
in  which  case  most  of  the  carriers  remain  in  the  lower  energy  portions  of  the  conduction  band  or 
the  first  75%  of  the  device  length,  the  movement  of  the  energy  bands  closer  to  the  gamma  valley 
with  increased  aluminium  composititHi  iiiqilies  diat  the  percmtage  of  carriers  in  the  subsidiary 
valleys  increases.  Thus  there  are  fewer  light  mass  carriers  available  to  conduct  current  and  the 
mean  electron  velociQr  in  the  aluminium  gaUium  arsenide  is  expected  to  be  significantly  below  that 
of  gallium  arsenide.  FIGURE  12  displays  the  percmtage  of  carriers  in  the  subridiary  valley  and 
the  corresponding  average  velocity.  We  point  out  diat  the  velocity  for  all  but  the  x  =  0.1  mole 
fraction  is  approximately  an  order  of  magnitude  less  than  that  of  the  mean  velocity  in  gallium 
arsenide. 

In  [9],  the  authors  point  out  that  in  the  low  field  region  of  the  structure,  which  is  dominated  by 
gamma  valley  transport,  real  si»ce  transfer  from  the  gallium  arsenide  to  the  aluminium  gallium 
arsenide,  which  is  the  thermionic  emission  of  electrons  from  one  device  layer  to  anodier  because 
of  heating  of  the  carriers  by  an  electric  field,  is  approximately  balanced.  In  the  hi^  field  regime, 
which  occurs  near  the  end  of  the  structure,  gamma  valley  electrons  transfer  to  die  L  valley.  Tbe  L 
valley  electrons  in  the  GaAs  easily  undogo  real  space  transfer  into  the  L  valley  of  AlGaAs  due  to 
deformation  potential  scattering  and  a  relatively  low  offset  barrier  for  the  L  ^dl^  electrons.  The 
electrons  in  the  AlGaAs  undo-go  further  transfer  to  the  X  and  L  valleys.  The  X  valley  carriers  do 
not  transfer  efficiently  to  die  X  carriers  in  GaAs  due  to  an  unfavourable  offset  voltage. 
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FIGURE  12  (a)  Percentage  of  elections  in  AKSaAs  which  have  scattered  into  the  X  valley 
for  the  indicated  aluminium  compositions,  (b)  Average  velocity  in  AlGaAs.  [9]. 


The  study  of  [10]  while  similar  to  that  discussed  above  «i^hasised  the  role  of  the  transport  in  the 
AlGaAs  region  subject  to  different  doping  distributions.  In  particular,  they  examined  an 
AlGaAsA3aAs  MODFET,  with  a  gate  length  of  0.5  miorons,  ctmt^y  placed,  with  a  spacing  of 
0.3  microns  from  the  source  and  drain  metallisation  regions.  The  aluminium  percentage  was 
30%.  TWo  cases  were  considered.  In  the  first  an  n^  AlGaAs  reeion  500  A  deep  was  doped  to 
1.8  X  10  >*/cm3,  followed  by  an  undoped  AlGaAs  spacer  layer  50  A  de^,  and  a  2000  A  un^ped 
gallium  arsenide  layer  In  the  second  structure,  the  first  400  A  of  the  AlGaAs  are  undoped,  die 
last  100  A  is  doped  to  5  x  lO'^Aan^;  the  remaining  parts  of  the  structure  were  unchang^  The 
specific  doping  levels  were  chosen  to  assure  identical  gate  capacitance  and  threshold  voltage 
levels.  A  self-consistent  calculation  was  then  performed,  with  the  electrons  subject  to  the  same 
scattering  mechanisms  as  those  discussed  for  [9];  namely  polar  optical,  intervalley,  ionised 
in^urity,  and  electron-electron  scattering.  Fermi  statistics  is  also  included.  The  calculations  woe 
perform^  to  reveal  the  differences  in  die  high  field  transport  in  the  AlGaAs  for  the  two  structures, 
and  thereby  to  reveal  the  relative  merits  of  one  against  another.  For  a  composition  of  30% 
aluminium,  the  offset  voltage  was  taken  as  256  meV. 

Th?^  calculations  woe  performed  for  a  drain  bias  of  2.0  V  and  two  different  values  of  gate  bias, 
40.4  V.  A  classical  calculation  of  the  conchiction  band  jxofile  under  the  gate  contact  for  a  gate  bias 
of  0.4  V  is  shown  in  FIGURE  13  (D-H^iT  denotes  delta  doped  structure,  U-HOflT  denotes 
uniformly  doped  structure).  Note  that  the  higher  conduction  band  levels  for  the  D-HEMT  in  the 
vicinity  of  the  gate  represent  the  presence  of  fewer  electrons  than  for  the  U-HEMT.  The  channel 
electric  field  profile  for  this  calculation  fcH*  both  structures  at  a  gate  bias  of  0.4  V  is  also  displayed 
in  FIGURE  13.  Note  that  with  the  exception  of  a  small  region  near  the  drain  contact  the  field 
I»x)files  are  nearly  the  same,  signifying  that  comparisons  of  the  two  structures  are  relevant 
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FIGURE  13  (a)  Conduction  band  profile  akMig  a  line  perpendicular  to  the  middk  of  the  gate  contact  for  a  dtain 
bias  of  2.0  V  and  a  gate  bias  of  0.4  V.  The  zero  coordinate  is  the  Schottky  contact,  (b)  Channel  electric 
field  profile  along  the  device.  The  zero  coordinate  is  at  the  source.  [10]  < 

The  observations  of  FIGURE  14,  which  are  profiles  of  total  electron  concentration  along: 
(1)  select  regions  of  the  AlGaAs  (at  a  distance  of  5S0  A),  (2)  the  heterointerface  (at  a  distance  of 
1()0  A  from  physical  interface),  and  (3)  the  GaAs  interior  (at  a  distance  of  19(X)  A),  iiklicate  that 
undo-  reverse  bias  there  is  little  distinction  between  the  two  stnictures.  However,  at  a  bias  of 
0.4  V,  there  is  a  reduction  of  electrons  in  the  AlGaAs,  and  a  reduction  in  t^e  parasitic 
contributions.  A  supplementary  calculation  of  [  10]  shows  transconductance  levels  that  are  nearly 
die  same  under  reverse  bias,  but  with  substantial  in^rovements  under  forward  bias. 


OMMoaOim) 

(a) 


0<atonot(|ini) 
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FIGURE  14  (a)  Relative  electron  concentration  within  various  layers  for  a  drain  bias  of  2.0  V  and  a  gate 
bias  of -0.4  V.  The  zero  coordinate  is  at  the  source,  (b)  As  in  (a)  but  for  a  bias  of  0.4  V.  [10], 


-2II  - 


75  High  field  transport  and  the  Gunn  effect  in  AlGaAsIGaAs  structures 

H.  NEGATIVE  DIFFERENTIAL  RESISTANCE  THROUGH  REAL  SPACE 
TRANSFER 

Ii  th«  discussion  of  transport  in  the  MODFETS,  we  raised  the  issue  of  real  space  transfer.  This 
phenomenon  was  first  discussed  in  1979  [11],  and  occurs  as  a  function  of  bias  in  all  structures  in 
which  transport  is  parallel  to  a  heterointerface.  Consider  the  superlaitice  sketch  shown  in 
FIGURE  15,  and  imagine  carriers  travelling  parallel  to  the  interface.  If  the  structure  is  modulation 
doped  then  with  the  doping  only  in  the  wide  band  gap  material,  there  will  be  transpon  in  the 
narrow  as  well  as  wide  band  gap  regions,  with  the  mobility  of  the  narrow  band  gap  material  being 
higher  because  of  a  reduced  number  of  scattering  mechanisms.  Under  application  of  an  applied 
bias  electrons  travelling  in  the  quantum  well  can  acquire  a  sufficient  energy  far  above  their  thermal 
equilibrium  value.  Electron-electron  interactions  will  help  randomise  the  energy  gained  in  die  field 
direction,  and  the  transfer  rate  of  electrons  in  the  gallium  arsenide  and  the  aluminium  gallium 
arsenide  will  be  determined  by  the  thermionic  currents  from  the  wide  band  gap  material  to  the 
narrow  band  gap  material  and  vice  versa.  As  in  the  discussion  of  the  energy  gained  in  the 
AlGaAsAjaAs  MODFETs,  electrons  in  the  AlGaAs  will  not  be  heated  to  as  high  an  energy  as 
those  electrons  in  the  narrow  band  gap  region,  and  we  may  expect  that  the  transfer  rate  at  high 
fields  will  not  be  the  same.  Indeed  more  carriers  are  expected  to  transfer  from  the  GaAs  to  Ae 
AlGaAs,  with  the  consequent  reduction  in  local  mobility. 


RGURE  IS  Sketch  of  a  superlattice. 

Thus  we  can  imagine  a  device  in  which  the  source  and  drain  regions  permit  transport  parallel  to 
the  interface,  and  two  components  of  transport  emerge.  First;  at  sufficiently  high  values  of 
electric  Held  electrons  will  undergo  k-space  transfer  to  subsidiary  valleys.  For  sufficiently  long 
structures  this  is  known  to  lead  to  negative  differential  conductivity  and  tiie  consequences  thereof. 
Second:  real  space  transfer  introduces  another  component  to  the  negative  conductance,  where  it  is 
possible  for  gamma  valley  carriers  of  sufficient  energy  in  gallium  arsenide  to  transfer  to  gamma 
valley  carriers  in  a  wide  band  gap  material  with  a  consequent  region  of  negative  differential 
conductivity.  It  is  also  possible  for  L  valley  carriers  in  GaAs  to  transfer  to  L  valley  carriers  in 
AlGaAs  with  a  consequent  region  of  negative  differential  conductivity.  Several  initial 
experimental  results  are  summarised  in  [12].  But  perhaps  one  of  the  more  unusual  results 
associated  with  teal  space  transfer  has  been  the  construction  of  negative  differential  resistance 
transistors  NERFETs  [13].  These  are  reviewed  in  [14]. 


I.  CONCLUSION 

The  alloy  aluminium  gallium  arsenide  is  a  material  whose  transport  properties  under  conditions  of 
high  fields  offer  a  new  dimension  in  terms  of  the  design  of  electron  devices.  Those  designs  in 
which  transport  is  perpendicular  to  the  interfaces  are  based  in  terms  of  tiie  energetics  of  the 
carriers  entering  the  regions  of  interest,  or,  e.g.,  in  the  case  of  HBT  upon  the  reduction  of  hole 
injection  into  the  emitter,  and  the  enhancement  of  electron  injection  into  the  base.  These  devices 
do  not  depend  upon  the  ^jecific  high  field  properties  of  the  AlGaAs  for  their  operation.  Parallel 
transport,  however,  does  expose  the  high  field  transport  properties  of  AlGaAs.  Here  while  the 
high  field  transport  properties  associated  with  electron  transfer  in  k-space  are  present  in  all  devices 
whose  structure  size  is  large  enough  to  sense  the  NDM  region,  this  is  not  the  feature  emphasised. 
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Rather  the  effects  of  negative  conductance  through  real  space  transfer  are  emphasised  through  an 
understanding  of  its  effect  on  the  operation  of  such  devices  as  MODraTs,  or  in  the  construction 
of  newer  devices  such  as  the  NERFET. 
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1 1 .  TRANSIENTS  IN  QUANTUM  WIRES 


Dissipation  in  rectangular  quantum  wires  was  studies  through  Monte  Carlo 
simulations.  Optical  and  acoustic  phonons  were  considered  and  H  was 
demonstrated  that  hot-electron  cooling  is  determined  by  cascade  emission  of 
optical  phonons  followed  by  a  slow  second  stage  of  indastic  dectron-acousdc 
phonon  interactions  (as  well  as  by  nonequilibrium  hot  optical  phonons).  A  copy  of 
a  recent  paper  is  enclosed. 
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Mooie  Carte  ahnulationa  of  hot  noneqniHbrinm  dechoo  rdaxatka  far  lectangnlar  OaAa  <|itaalnm 

bM-olecizoQ  oooUng  dynamica  is  detennined  by  cascade  cmltdon  of  optical  and  ediifaiii 

ftmog  dqrendence  on  the  exdtatioo  energy.  The  second  (slow)  rdaiatioo  stage  ia  oonmlled  by 
aonagly  inelaatic  electron  iuenctions  with  accostic  pbonois  as  wdl  as  by  noneqnilibfiDm  Oxx) 
optical  phonons.  The  relaxation  times  obtaiaed  far  our  shnnlatioos  ate  in  good  agreemm  with  the 
remits  of  recent  Inminraoenre  e:q>6timents.  At  low  dectron  oonoenliatioos  ndrere  hot 
effects  ate  negti^de  the  cascade  emiadon  of  opticd  phonons  may  lead  to  tire  overcooUtig  of  the 
electron  system  to  tempaatme  below  the  lattice  tempenunre.  Them  elections  then  slowly  (doting 
tens  of  picoaeconds)  r^  to  equilibrinm  doe  to  tiie  fanenction  wiflt  acoustic  phonons.  At  oestam 
excitation  energies  strong  ittermbbsnd  electron  scattering  by  opticd  phonons  leads  to  dectron 
redistribntion  among  mbbands  and  intesmbband  population  hwwtsions.lftiie  electron  concentration 
exceeds  10^  cm~',  hot  phonon  effects  come  into  jday.  In  contrast  to  bulk  maieriais  and  ywnw" 
wells,  hot  phonon  effects  in  qiiantnm  wires  exhibit  strong  dqtendence  on  tire  initid  broadening  of 
the  enagy  disnibatian  of  die  electnons.  The  rtxy  initid  detaron  gas  relaxation  stage  far  qnanttm 
wires  is  faster  in  the  presence  of  hot  phonons,  wfaDe  for  t  >0.5  pa  the  hot  phonon  thesmalization  ihne 
defines  the  characteristic  electron  cooling  time. 


LINTROOUCnON 

Potentid  possibilities  of  utilizittj  nniqDe  properties  of 
<]uasi-ooeKfimaisiond  (ID)  semiconductor  structures  ^teak- 
like  density  of  sates.  Ugh  padcr^  densi^,  U^  operation 
frequencies)  for  the  devdofanem  of  a  new  generation  of 
electronic  and  optoelectronic  devices  have  stimulated  en¬ 
hanced  interest  in  investigating  die  ooneqdlUvinm  electroo 
relaxation  processes  in  them  sttnctnres.  Bietgy  (Bss^pation 
processes  trf  hot  dectron  gases  defirmsudi  characteristics  as 
device  Hug  speed,  efSdency,  gain,  transport  character- 
iatics,  noise,  etc.  Aldnugh  there  exist  a  nnmber  of  publica¬ 
tions  dealing  with  relaxation  processes  in  ID  dectron 
gases,'*'*  dm  relaxation  dynamics  of  photoexched  carriers 
under  highly  nonequilibrium  conditions  have  recdved  less 
•wtinrfaw-  Recent  experimentd  investigations  cleariy  demon¬ 
strate  different  behavion  of  photoexcited  efectron-hde  plas- 
mas  far  ID  and  2D  systems.’*'’  Thne-iesdved  huninescenoe 
measurements  hr  ID  quantum  wires  (QWb)**'’  indicaie  that 
nonequilibrinm  carrier  relaxation  to  the  lowest  ID  state  is 
rather  slow  compared  to  that  measured  in  quantum  weds. 
Thus,  with  current  interest  in  devdofring  a  new  generation  of 
devices  baaed  on  ID  aemkonductor  structures  it  is  inqxntant 
to  ■■■Mtutwwt  the  temporal  evolution  of  rdaxation  of  non- 
eqdlibiinm  caiiiers  following  inidal  exdtatioo.  TUs  evolu¬ 
tion,  especially  its  initial  stage  (first  few  picoseconds  follow¬ 
ing  polm  excitation),  is  of  great  importance  for  device 
qipUcations,  and,  hr  particular;  for  U^-qpeed  photonic  de¬ 
vices.  There  are  several  important  aspects  of  electron  relax¬ 
ation  in  ID  quantum  wires  to  be  considered  in  great  detail 


At  tlM»  OjiiMriiig  wmparliiiiM  nf  mnW  npntjtartmwir 

vices  (30-300  K),  often  the  only  important  energy  and  mo- 
mwanm  rdaxation  mechanism  is  the  dectron-phooan 
interaction.^  In  comraat  to  3D  or  2D  systems,  dectron- 
electron  pair  collisions  in  QWIs  do  not  affect  dectron  idax- 
ation,  whereas  dectroohphoaoii  interactions  remain  very 
mrong.*^'**'*  That  ia  why  at  low  dectron  concentration  ndien 
the  hnrisahhand  etocoon-dectron  interaction  is  weak,*^ 
doctpott  pbODOos  dctHuiiics  tibc  cocivo  ^dccttoo 

nlauiktt  dynamics.  It  qipeoB  ibttpboom 
from  the  same  m  fat  bulk  materials.^’’ Along  witii  election 

iliifogiiiifWb  miy  exist  fl 

2.a 

QBBn  UroB  IGQQSDC  BDO^lOwn  vC90BOHK 

QWU  is  essentially  famlastic  due  to  die  lack  of  translatinnal 
qrmmetry  and  the  resultant  .nftfttny  of  momentum  con¬ 
servation.  Tberefote,  far  considBtiug  the  dectron  ^mamir  b^ 
havlor  in  (2WIs  it  is  extremely  inqioRam  to  allow  for  the  red 
phonon  qiectinm  far  QWIs. 

When  dectrons  are  excited  wdl  above  the  bottom  of  the 

/■nnHui-riaw  limil  tiny  iiiJ«t  irf.  CMrMto  laiitMlnn  <rf  plmwnm 

and  drive  tire  phonon  qrstem  out  equilibcinaL  It  fe  now  a 
commonly  -neperd  nnrtmt  iMw«.pinti.<iini  (htu) 
{dionons  strongly  affect  dectron  transport  and  rdaxation  in 
bulk  materials  (see,  e.g..  Ret  2^.  Ifot  phonon  effects  also 
exdain  observations  of  very  slow  electron  coding  in  quan¬ 
tum  weUs  following  subpicosecond  pholoexdtation  of  hot 
electroos  (see,  e.g..  Ret  25).  The  hot  phonon  problem  far 
QWIs  has  been  addressed  previously  in  Re&  2,  6,  and  8. 
However,  the  kinetic  approach  used  in  Ret  2  does  not  pro- 
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vide  infonnaiioo  about  relaaation  dynamics,  udiile  in  Ref.  8 
die  pecnUarides  related  to  ID  nanire  of  noneqnilihrinm 
elecoon-pbonon  systeu  have  been  ovoiooked. 

In  contrast  to  30  or  2D  symems,  the  infnwihhanri 
electron-electron  pair  interaction  in  QWIt  leads  only  to  mo¬ 
mentum  exchange  between  interaction  tnrtitfingnishable 
electrons  and  does  not  contribute  to  the  relaxation  process. 
Hence,  at  least  during  initial  relaxation  stage,  the  deciron 
gas  cannot  be  described  by  a  Maxwellian  distribotion  func¬ 
tion  and  it  is  neceassiy  to  investigate  electron  relaxation 
without  any  a  priori  asaumptioos  about  the  etectron  distribu¬ 
tion  funfrion.^  Moroover,  the  ID  nainre  of  optical  phonons 
in  QWls  results  in  some  specific  pecuHaritiet  of  hot  phonon 
imtiHim  tint  should  strongly  modify  hot  phonon  effects  in 
QWIs.^ 

Ttl  r*T*»’  «h*  «t«nnl«ri«Mi  nf  Atr,  iwleiWion  of  hot  nen- 

juptniiriiiiB  electrons  has  been  carried  out  allowing  to  all 
the  specific  aspects  of  electron  dynamics  in  QWIs  mentioned 
above. 

a.  MODEL  AND  HETHOO 

Simulation  of  hot-electron  relaxation  has  been  per- 
tomed  in  a  rectangular  GaAs  QWIs  embedded  In  AlAs,  We 
empiny  an  ensemble  Moroe  Ctrio  tnrhnli|ne  specMOy  sailed 
to  ID  decnon  simulation.*^  A  twodimnskmaL  htifadtefy 
deep  square  posentid  wdi  confines  dectrons  in  die  QWI 
whb  a  mnltisnbband  energy  structure.  The  hot-electron  en¬ 
emy  distipstion  modd  includes  electron  interactions  widi 
confined  longhudinal-opticd  (LO),  locaHaed  imtoe  Qnler- 
toe)  opdcd  (SO)  phonons,'^'**  and  fanIk-Bice  acoustic 
ptaoBon^  as  wdi  as  nonequilibriam  opticd  phonon  popola- 
dons.  We  have  indnded  indasticity  of  dectroiHKOiistic  pho¬ 
non  aoMering  in  a  QWI  in  fiill  detail  ushig  the  technique 
proposed  in  Ref.  23.  The  initid  distribution  of  excited  dro- 
trons  among  tabbands  is  considered  to  be  defined  by  the 
density  of  states  to  a  given  excess  energy  hi  eadi  sabband. 

We  Stan  the  simulation  of  decnon  rdaxaikm  after  the 
inidd  exchation  by  a  shott  pulse  with  a  dmation  of  0.1  ps. 
We  have  not  simulated  electron  relaxation  in  coherent  regime 
(KlOO  ft)  uhich  lequhes  a  quantum  aserhaiticd  descrip¬ 
tion.  Insi^  we  have  focused  our  attention  on  the  thw 
range  t>0.1  ..ps  when  dectrons  cm  be  treated 
semidatajcally.*'”  We  are  interested  primarily  in  pecnliari- 
ties  of  the  dectron-phonon  intenction  in  Therefore, 
we  do  not  take  into  account  the  electron-hole  interaction. 
TUt  dtuation  could  take  place  sriiea  electrons  are  photoex- 
dted  fiom  netf-fflonoeneigetic  imporiQr  kvds.  The  initid 
state  to  electton  rdaxation  accounts  to  the  brnadeniiig  of 
the  electron  energy  distribution  due  to  tsro  effects:  0)  uncer¬ 
tainty  in  electron  initid  energy  due  to  the  dioit  electron  av¬ 
erage  lifetime  at  the  excited  levd  (A^10~^  eV  to 
Af~10~‘^  s);  (ii)  spectral  broadening  of  the  exciting  pulse 
srith  duration  of  the  Older  <rflO~“s.  In  accounting  to  fliese 
effects  we  assume  that  they  both  lead  to  a  Gaussian  distribu¬ 
tion  of  electron  eneigy  at  (*io.  wliidi  corresponds  to  the 
end  of  die  exdntion  pulse.  A  Gaussian  broadening  ftctor  is 
used  instead  of  a  Lorentdan  (typically  incoipanted  into  ide¬ 
alized  tbeoreticd  model^^^  to  prevem  die  unrealistically 
huge  qaead  of  electron  eneigies.  Hus  approadi  hu  been 


justified  to  the  electron-phonon  imeracrion  in  a  QWIs  dem¬ 
onstrating  diat  the  Golden  Rule  fesmaliam  m^  be  retained 
hy  convolviBg  a  Ganarian  hmsAiniiig  ftmrtimi  wwif  Iwjiig  ■ 

constant  broadening  factoc”  We  vary  the  excitation  energy 
which  coiroqwnds  to  fee  center  of  a  tx— «»«i  dianftu- 
tion,  as  svell  as.  Ac,  the  half-widdi  of  this  diatribntioa 
Noneqnilihrinm  phonons  have  been  by  calcu¬ 

lating  the  phonon  ocoqiation  number  venns  phonon  wave 
vector  (phonon  distribntion)  widdn  the  Moine  Cario  proce¬ 
dure.  In  aooordanoe  wife  the  ID  namre  of  opticri  phonons  in 
QWIs,  the  mcrement  of  phonon  occupation  after 

each  emisaion  (sign -f)  or  absorption  (sign -)  event  is  given 
by  the  term  ±(2n/Aq)(a/N),  when  Ag  is  die  step  of  die 
grid  in  q  apaet  used  to  teoord  the  AT,  hlatogtam,  a  Is  the 
electron  concentration  per  unit  length  of  a  QWI,  and  W  is  fee 

In  tllC  iilBIllfltkXL 

lit  Unma  rtio  atinnUrinM  «f  anHr  — <t  OH  imiMMptMK- 

rinm  electron-optical  phonon  systems,  fee  mesh  imecval  to 
the  phonon  occtqiaiion  number  Ag  is  not  a  crndal  parameter; 
givea  that  the  interval  is  much  Im  than  the  g-apace  region 
populated  by  noneqpiilibtinm  phonons  which  can  be  eaaQy 
estimated.  Hus  is  due  to  fee  feet  that  thephononteabaorp- 
tion  rate  depends  on  the  integrated  (average)  Of riipnncy  over 
die  emise  rmion  which  is  not  ond^  aemdtive  to  the  mesh 
iniervaL  However,  in  ID  systems  as  we  wiU  see  in  a  due 
the  mbiocptioii  d^wodt  onfyoo  Ab  IikhI  vsIob 
of  phonon  occupancy  A,  at  an  appropriate  g  value.  Hiere- 
for^  as  fee  meA  fanerval  becomes  amalhr,  bofe  the  local 
oocapancy  and  the  reabaorption  me  become  laiger;  This 
problem  is  paiticulariy  importut  when  «nn«i<fariiy  oear- 
monoenergelic  electron  excitation.  Hme  are,  of  oonrie, 
physical  Urnhs  on  the  magnitude  of  Ag.  These  limits  follow 
fiom  the  uncertain^  in  fee  phonon  longinidinal  wave  num¬ 
ber  due  to  the  finite  length  of  the  (2WL 

Wi  have  taken  a  QWI  of  lengfe  £.>10  tm,  so  that 
Ag>2ii/£.>6xl(>’  cm~*.  Hot  phonon  thcrmaliration  dne 
to  the  decay  of  optical  phonons  into  acoustic  phonons  is 
taken  into  account  by  recalculating  to  every  mesh  imer- 
val  at  the  end  of  ead  time  step,  b  has  been  demonstrated** 
thfeft  ffhf  pbonoo  th*— riwMb  ^  In 

stroctares  di^eadt  weakly  on  die  siaes  of  die  atnemae  and  is 
dose  to  the  bulk  value.  For  simulartnns  at  r>30  and  77  K. 
we  have  used  the  value  r^>7  ps,**  The  time  step  in  our 
simalationt  has  been  dmaea  to  be  smaller  fean  die  average 
time  between  two  events  of  electron  acatteriiig  by  optical 
pluyinnf  »iwl  Icm  *t««i«  fee  ptumnw  tlumiMliMtfaM  riiw 
We  have  not  taken  into  account  die  increase  in  die 
aconaticphononpopulationasaresukoftliedec^yofnon- 
equilibrinm  optical  phonons.  There  are  tsro  reasons  to  tins. 
I%st,  fee  boikhqi  of  nonequilibrinm  optical  phonons  oocnn 
only  in  a  very  nairow  r^ion  of  the  BriDonin  zone  (near  the 
zone  center),  so  that  over  the  entire  zone  the  average  oocn- 
pation  number  increases  ordy  negligibly.  Hds  is  tioe  to  sys¬ 
tems  oi  any  dimensionality  sfaioe  die  electrons 
srife  phonons  populate  only  the  center  r^lon  trf  die  Brillanin 
zme.  Second,  the  acoustic  phonons  in  QWb  embedded  in 
suirounding  materials  sridi  similar  dastic  propetties  (GaAs 
in  AlAs  in  our  case)  may  penetiate  dBongii  GaAs/AlAs  in- 
tetfeces  and  escape  ftom  to  QWL  Therefore,  we  have  ex- 
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FIC.  I.  TlBe  ««olacioa  of  te  qictd  elticw  onagy  la  •  QWI  wtt  cna 
«eci>aBl30x2aoA*ifter«mrtMrii»aw»ola«lce<amieiiww.(»)r-77K 
and  (b)  r>300  K.  Gave  1  la  (a)  coaw|wailatDelBcaonaaflMllwieaeny 
^-20  awV.  1 42  awV.  3.  «7  aeV.  4. 100  aocV.  la  (b)  1. 42  aaeV.  2. 67 
aaeV.3. 100  aaV. 


oeUent  tfaennal  coadiictivity  and  the  QWI  dioald  not  be 
heated  modi  mote  than  ite  whole  OaAa/AlAa  mctme. 
Given  that  the  sumninding  AlAa  is  mfUrinntiy  nuaiive,  the 
increase  in  temperainre  wooU  be  nqili^bie  even  if  die  QWI 
strongly  mdiaies  aconatic  phonons. 

m.  RESULTS  AND  D  SCUSSKM 
A.  Low  nlnctron  concutruiotie 

Let  ns  fitst  consider  electron  concentrations  less  than  10* 
cni~‘  ndiere  nanequilibrium  phonon  effects  can  be  n^Iecied. 

Calcnlationa  with  varions  escitation  energies  show  that 
in  the  dme  scale  of  10~*  s  electron  relaxation  exhibits  one  or 
two  disringnishabie  stages.  Hgnre  1  demonstrates  die  dec- 
tron  coolhig  dynamics  in  a  QWI  with  cross  section  150x230 
A*  for  lattice  temperamres  of  r»77  K  and  rB300  K,  as  wen 
as  for  different  electron  excitation  eneigies  ^  coomed  from 
the  bottom  of  the  lowest  condoction  aobband.  The  mean 
election  energy  plotted  on  the  vertical  axis  in  Rgs.  1(a)  and 
1(b)  is  calcnlated  reladve  to  the  bottom  of  the  first  sobband. 


Thns  it  consists  of  two  parts:  the  M»«**i<*  energy  conespond- 
ing  to  one  d^ree  (rf  freedom  in  a  QWI  ami  die  imersnbband 
sepandon  energy.  The  decnon  ncitaiian  enogy  ^  hu 
bera  varied  from  20  to  100  meV.  lUs  hrrplies  that  for  this 
panicnlar  cross  secdon  of  the  QWL  rip  so  the  three  lowest 
snbbands  can  be  occrrpied  by  electrons  at  the  initMi  rim*!, 
(H).  One  can  see  from  Rg.  1(a)  that  for  electron  excitation 
at  20  meV,  the  electron  gu  coding  is  slow  ('‘slow’’  stage). 
The  “fut”  stage  in  die  mean  electron  energy  dependence  on 
time  is  observed  when  electrons  are  excheddiove  the  opticd 
phonon  energy  (fiwu>  or  Aw,o>  where  Awu  and  Aw^  are 
energies  of  LO  and  SO  phonons,  reqncdvefy).  Electtona  ini* 
dally  (in  die  snbpicoaecond  dme  scale)  cod  down  »«««■§ 
Ihdr  energy  due  to  the  interaction  with  ruvjfrf  piirmn— 
Since  the  optied  phonon  abaoqidon  events  at  tMnp.— 
T-77  K  are  negli^y  rare,  the  electron  gu  ndaxadon  dy¬ 
namics  is  deSemined  by  die  »«««>— 1«««  of 
with  characteristic  timu  s  and  t,^-«10““  s 

^  dectron-LO  and  eleciron-SO  phonon  hneeaction.  respeo- 
dvdy).  It  is  worth  to  mention  that  ahrafrat  aoraeriaiHfaiiun 
carrier  relaxation  with  characteristic  dmu  of  the 

same  order  of  magninide  have  been  experimentally  obeenred 
in  dme-iesdved  piwitnimiiiiioiHw*  and  r«tti«Mrwiiimtw«Mu 
cenoe  measnrementa.*’^*’**  At  low  optical  excitation  levda 
the  dhaip  fine  of  band-edge  phoeohrmhuiacmce  ocean  rhrrlng 
d»  laser  excitation.****  This  implin  that  bd  catrien  low  the 
ni4)or  portion  of  their  excen  energy  daring  the  dme  mnch 
shorter  than  die  exchathai  pobe  of  25  pa.  The  andysis  of' 
low-iempexamre  cathodohminf  steme  spectra  saggasta  that 
carrier  capture  and  relaxation  to  die  bottom  ■<»*■■«««  in 
GaAa  QWIs  grown  on  nanpiaMr  nbatratea  ooenrs  in  a  anb- 
picoaecend  time  acale.**  At  a  lattice  teniperamw  of  r»300K 
the  electroo  coding  r^mamics  is  intinwinwH  atrongiy  by  op* 
tied  phonon  absorption  whkdi  rednou  the  dectron  gw  cod¬ 
ing  rate  [Rg.  1(b)]. 

Ihe  duration  of  die  “fut"  relaxation  stage  u  wdl  u  the 
entire  electron  gu  cooling  dynamics  for  exhibits 

a  strong  dependence  on  the  exchadon  energy.  As  rfiacussed 
below,  when  electrons  are  excited  just  above  the  LO 
energy  they  cod  down  to  the  bottom  of  the  first  ardiliand  on 
a  mbpicosecond  time  scale  [curve  2  in  Hg.  1(4  and  curve  1 
in  1(b)].  Electrons  eaut  optied  phonons  and  occupy  whm 
near  the  sobband  bottom.  Therefore,  the  mean  electron  en¬ 
ergy  drops  below  that  fori^*20  meV  [curve  1  faiHg.  1(a)]. 
For  a  lattice  temperature  of  r»300K  we  observe  anomahrea 
coding  ^mantics  when  eketrou  occur  bdow  die  dieinul 
eqoilibiiom  energy  [curve  1  in  Fig.  1(b)].  Overcoding  die 
electron  gu  occurs  if  the  dectitm  exdmdon  energy  frDs  hno 
the  range  iiwio<t„<Awio+kgTI2,  vidiere  kgT/2  is  die 
electron  kinetic  energy  at  a  given  temperature  T  correspond¬ 
ing  to  one  degree  of  freedom  in  a  QWL  At  lower  tempera- 
tniu(r««77K)  die  transient  electronovelcoding  disappears 
becanw  die  chosen  broadening  of  dectron  initid  energy  rSs- 
tiibntion  exceeds  die  electron  eqoiHfaiinm  energy 

*,r/2. 

The  “slow"  stage  of  dectron  rdaxadon  is  controlled  by 
tiiedectronii]teiactionwidiaconsticpiionooB.Onrca]cala- 
tkms  demonstiaie  thru  electron  gu  dxnnalixation  procen  in 
a  QWI  of  cron  section  of  150x250  A*  lasts  abow  1  u  at  a 
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lattice  tempenoire  r«30  K  (Fig.  1(a)].  and  30  pa  at  7-300 
K  (Rg.  1^}].  TUs  ime  dependa  atrangly  aoc  ooty  on  the 
iattioe  teaqMiatiiie  but  abo  on  the  craas  aoctkn  of  a  QWl  at 

Am  lim  prait»«<*  |il»ntinn  r— ^  Tho  «f  amm. 

dc  ^  fcsttoriiig  is  In  2  ftp  friffctwn 

wxjHng  ^namica  in  a  QWl  with  a  croai  aectian  40X40 
ooapned  with  the  oooiiiig  dynamka  in  a  1S0X2S0  A^  QWL 
The  dection  aiagf  leiaxation  due  to  the  imeraction  with 
aconadc  pbonone  it  aancb  faster  in  the  ddn  QWI  u  a  remit 
of  two  fKtocK  (i)  the  aconatic  phooan  acawrring  nte  is 
ranghfy  imeiaeiy  praporiiaoal  to  die  aoas  section  of  a  QWI, 
and  (ii)  the  inelasticity  of  the  etectroo-aooastic  phonon  haer- 
action  also  increases  with  the  decrease  of  dw  cioat  sectioa” 

Ihe  reiaaatioo  tunes  of  the  same  older  of  niagninide 
have  been  derived  from  the  time-resolved  photolnmmes- 
oenoe  measurements  at  low  eacitation  levela.”  Alter  initial 
bot-canier  relaxation  below  optical  phonon  taatf.  the  ftr- 
ther  evohdion  of  the  band-edge  i—  Une  shape  is 

chamcteriaed  by  the  lime  of  the  older  of  huntheds  of 
piroaennnds“  Therefore,  the  election-aeonstic  phonon  inter¬ 
action  might  be  responsible  for  the  dme  evolution  of  hnni* 
noscenre  spectra.  Rmgh  esdmases  yield  electron  dienaaliza- 
don  dme  due  to  niseraction  with  acoustic  I**—*—  aS  die 
Older  of  500  pa  for  the  strucnne  parsmetert  and  temperature 
(r-5K)ofRel  12. 

One  can  see  from  Hg.  1(a)  and  Fig.  1(b)  dmt  the  electron 
dieraial  eqnilibrinm  eneigy  to  7—300  K  is  larger  than  could 
be  eapectM  from  ^,7/2— 13  meV,  while  to  7-77  K  it 
practically  coincides  with  f  p7/2-3  J  meV.  Ihe  in 

thermal  egnilibrium  energict  conies  from  the  calcalstion  of 
the  dectron  mean  energy  in  QWIs  with  multimbband  enogy 
siinctare.  Appmaimately  one  third  of  electiQos  occupy  upper 
snbbands  in  the  equilibrium  state  at  a  lattice  temperature 
7-300  K  due  to  the  Bohzmann  distribotioo.  Ihe  mean  elec¬ 
tron  energy  indndes  the  electron  gas  Unedc  enogy  {kgT/2) 
of  free  morion  along  the  wire  and  die  energy  r^reaenting  the 
spatial  quandsadoo  (separation  between  mfabands)  in  two 


3lp 


HO.  3.  Tbae  avoimaa  at  ite  nrnipenrji  of  Uw  im  iSbema  fer  We  woe 
mean  pimrirri  aO  wftnHou  caoite  la  I;  (■)  r>77  K.  (b) 
r-300K. 


other  directioos.  In  the  extreme  Hmit  of  thick  QWIs,  when  a 
large  lannber  of  robbands  becomes  occupied,  the  electron 
mean  energy  tends  to  3kgTf2  coneaponding  to  the  3D  dec¬ 
tion  gas. 

Shnalation  of  bot-dectron  relaxation  dynamics  in  QWb 
demonstrates  that  fanermbband  dectron  primarily 

by  opdcd  phonons  leads  to  a  significant  cairier  redistribution 
among  mbbands  (Hg.  3).  When  electrons  are  exdted  wdl 
above  the  bottom  of  the  second  mbband  (^-100  meV) 
mnhiple  dection  tnnsitions  between  various  nw-tiih  due 
to  imeiactian  with  opticd  phonons  lead  to  a  nomnonotoaou 
dme  dependence  of  the  rdadve  occupancy  of  die  fint  (low¬ 
est)  mbband  [curves  3  in  Hgs.  3(a)  and  3(b)].  Aconadc  pho¬ 
non  scattering  is  also  reqxmsible  for  dectran  inteimbbond 
transitions.  For  the  case  where  the  crom  section  is  1S0X2S0 
A^  die  sqwrstion  between  the  first  and  the  second  mbband 
is  leas  dian  die  optical  phonon  energy,  so  that  electrons  can¬ 
not  be  scattered  from  die  bottom  of  die  second  mbband  by 
the  emissinn  of  qpdcal  phonons.  Accordingly,  at  low  tem- 
peratuies  electron^  be  “trapped"  in  die  second  mbband. 
They  are  skiudy  (with  characteristic  time  of  tens  and  hnn- 
dieds  of  picoseconds)  released  from  it  dne  to  in«eiwiHi«nrt 
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nan  te  aooad  labtaal  The  laaka  Mpvaaaa  k  r-?7  K. 


decmo  acanaring  by  acooatic  pbanana  (catvea  3  and  4  in 
Hg.  3(a)].  In  lUa  caae  (he  aecood  aabband  aervea  aa  a  hot* 
dectm  leaenrair  and  aignifiandy  alowa  dectno  cooling. 
Even  in  the  caae  when  moat  clecticaa  are  exdied  into  the 
liiat  anbband  (Vi«>20  meV),  a  amall  fiaction  of  them  (from 
the  Ugh-eneigy  tail  of  the  Oanaaian  exdtadon  energy  diatri- 
butioiO  are  initially  acanered  by  aconatic  pbooooa  10  the  aec- 
ond  anbband  and  then  return  to  die  fiixt  one  [cnrro  1  in  Fig. 
3(a)].  The  eneigy  of  the  pfattean  [carve  3  in  Rg.  1(a)]  viitn- 
ally  comddea  with  the  poaition  of  the  aeoond  anbband  widi 
reqiect  to  the  fiiat  anbband  bottom  (27  meV)  indicating  that 
clectrona  are  "napped”  diere.  At  biih  temperatnrea  [r»300 
K.  Fig.  3(b)]  eiecoona  “eacape”  fiom  die  npper  anbbanda 
dne  to  hneiaabband  abeotption  of  optical  phonona  aa  well  aa 
atronger  imennbband  acooatic  pbora  acatteiing,  and  reach 
an  etpdHhtinm  diatribotioo  among  anbbida  in  30  pa. 

Under  certain  excitation  conditioiia  intwaabband  dec* 
non  acatteiing  by  optical  phonona  any  lead  to  fanwaobband 
popalatton  invenkm.  We  obeerve  an  imeraubbaiid  popaladon 
invetaioo  when  two  condidona  are  mec  0)  aepmadon  be¬ 
tween  two  loweat  aabbands  in  the  QWI  it  ien  dim  minianmi 
optical  phonon  (LO  or  SO)  eneigy.  ao  that  electtona  cannot 
be  acanried  from  the  bottom  of  the  aeoond  aabband  by  the 
emtiaion  of  optical  phonona  (we  demonitrme  lesaka  for  a 
QWI  with  ciomaectioo  150X250  where  dda  condition  ia 

fidfiOed);  (ii)  eiectrona  are  excited  jaat  above  chaiacteriatic 
energy  where  cn  ia  the  energy  ol  the  bottom 

of  the  aecond  Mbband.  Dae  to  a  rignificant  difference  in  the 
nmdier  of  final  atalea  (peak-Ooe  denaity  of  ataiea  near  eadi 
aabband  bottom)  dectrona  from  bodi  the  fiat  and  the  aeoond 
anbbanda  are  acatiered  predominantly  into  die  aeoond  aab¬ 
band  bottom  after  the  emitainn  of  LO  phonona.  Thaa,  the 
of  eiectrona  at  the  bottom  of  die  aecond  aabband 
exceeda  the  annber  of  eiectrona  at  the  bottom  of  die  first 
aabband  and  a  atrong  imetsnbband  popalatton  inveratoo  oc¬ 
ean  near  the  center  of  the  Btilloain  xone  (A^O).  Rgnre  4 
picaenta  the  diatribatton  of  eiectrona  in  momentnm  apnoe  tor 


the  two  ioweat  anbbanda  3  pa  after  Elecoooa  in 

the  first  anbband  are  adll  hot  (wave  mmiben  4>2x  10^  cm~' 
on  carve  1  in  Fig.  4)  after  ««wi««inB  of  opiied  m»i 

they  relax  to  the  bottom  of  the  mbbmid  by  inenedng  with 
acooatto  phonons.  Eiectrona  in  the  aecond  anbbaad  occqiy 
atatoa  with  smaller  wave  veeaxs  near  die  mbband  bottom. 
TUa  popalatton  inversion  near  the  «««*«>  of  the  BiWUiin 
xone  (A-O)  lasu  abont  10  pa  at  a  hmioe  temperature  of 
7**77  K.  TUa  dme  is  by  onj  iniennb- 

band  electron  x*****'^  by  acoottir'  Inlecmbbmid 

electron  aratirring  by  acooatic  phonona  is  leaponsihie  tor 
deciron  ideaae  from  the  aecond  «n>*»nH  at  low  tempeta- 
tnres  where  optied  phonon  abaorptton  ia  vinaalty  fioaen  oat, 
while  intiambband  acooatic  phonon  leada  to  the 

thfarmaliialion  of  the  electron  tfiatribntion.  As  one  can  see 
from  Hg.  4  the  popalatton  mventon  at  aaaaU  electran  wave 
veettn  is  redneed  dae  to  pieaence  of  acme  fiaetton  of  dec- 
trons  near  the  bottom  of  the  first  aabband.  The  mnnber  ami 
eneigy  of  them  ekcirotia  depend  atrongiy  on  the 
r^inie.  Doe  to  the  deciron  eiMgy  tfis- 

iribation  some  eiectrona  ftom  the  Ugh-eneigy  tail  can «— <» 
two  optical  phonona  and  cool  down  to  the  bottom  of  the  fiist 
mbbmid.  Tbas,  m  the  eleciron  inittol  energy  in- 

creases,  the  occapatton  of  states  with  small  wave  vector  in 
the  first  mbband  also  increaaes,  and  the  effect  of  popnladoa 
invenion  decreases. 

Be  noi  pnonofi  BiiBra 

Dae  to  optied  phonon  qaandiation  and  die  reanhant  ID 
momeamm  conieivation  in  qaaniam  wires,  dectront  can 
emit  or  absorb  optied  phonons  with  wave  vectors  which  are 
atriedy  defined  ^  the  deciron  momeranm  and  the  phonon 
energy.  In  general,  the  phonon  wave  number  is  defined  by 
the  energy  and  momentum  conservation  eqnaltona  and  is 
given  by 

g-  V*^+*'*-2U'  cot  e .  (1) 

where  it  it  the  dectren  wave  number  betore  TaiaKritu, 
k'  *  is  the  deciron  wave  number  after 

absorption  (sign  +)  oremiaaion  (sign  -)  of  the  optied  pho¬ 
non  of  fieiineacy  aib>  ><><1  f  i*  (Bo  angle  between  electron 
wave  vectors  betore  and  after  acattering.  In  ID  atruennea 
there  are  Just  two  find  atates  tor  acanered  eiectrona:  forward 
araiierlug  with  cos  ^1  or  backward  afjfteriiig  with  cos  d 
*-l.  Conseqnendy;  there  are  two  posaiile  phonon  wave 
vecaora  avaflable  fbr  emitainn  (and  two  tor  ahaorptloiO  by 
my  single  electron: 

ga-|*+*'| .  (2) 

In  oontraat,  in  quanhim  wdb  (or  bofic  nuterials)  dne  to  ex- 
iatenoe  of  adtfithwd  degiee(t)  ci  fineedonu  cos  dean  take  aqy 
value  in  the  range  (- l,-l- 1),  so  thd  there  is  an  enttoe  range  of 
a  phonon  q  vafaim  from  |k— k'|  to  k+k'  mSUtile  tor  deo- 
tron  faneracdons. 

Tbeietore,  eiectrona  in  QWIa  having  qgxeciably  difier- 
ent  eneigiea  gf.iiHate  noneipiifibrinm  phonona  in  different 
nmrow  q-tpace  regions  wind  do  not  oveilap.  hi  nmi,  dieae 
phonona  can  be  reabaoibed  only  by  the  dectront  diat  have 
gendated  them,  nnitim  in  bulk  maieiiala  imd  in  ^lantnm 
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ateUs  aAm  ckctraot  caa  teahaab  ptnoona  Mated  by  fsdia 
dectrooa.  Conaeqnently.  electroot  wfakli  have  dMerent  eaer- 
gict  ttM>  Uyipii  tBiMW|Q6Bt  VD* 

abaorptiaa  of  optical  phoaona.  Tlaia.  eledxoaa  and  Ifae 
pboom  aaaociaiBd  with  then  (with  appropriate  warn  vec- 
lon)  are  iaolated  from  other  eiectroa-pboaao  peka  if  elec- 
trooB  are  in  diflBreni  ataies.  Tbns.  for  ID  atraetnea  we  eiimi- 
naie  one  croaa-oanelaiiaa  effect  which  ia  aiwaya  paeaett  in 
the  nonegmlihrium  eiectreii-phooon  ayitea  in  aataiiala 
and  ia  tpianBnn  welia  and  which  ia  important  ediea  conaki- 
ering  electric  noise.^ 

Another  conaequence  of  the  ID  naiare  of  the  electron* 
phonon  BMerection  in  QWb  ia  that  the  leabaonxion  probabil¬ 
ity  for  each  aingie  electroa  in  the  QWI  does  not  depend  on 
the  iniegraied  phonon  occupation  number  but  oniy  on  the 
local  occupancy  at  a  oeitaia  q,  TIub  leabaorprion  probabiliqr 
decreaaea  aa  the  phonon  diiiribmion  apreada  over  q  apace 
(given  that  the  initialed  occupancy  ia  by  die  con- 
ceniraiion  of  eaciied  electrona  and  tenuina  conatam).  The 
apeead  of  the  nonegnilihrinm  phonon  popnlaiion  in  g  qiace 
reauha  from  the  broadening  of  the  eleciron  energy  diatribo- 
don.  Aa  a  leanh,  the  reabaoiption  rate  and  hot  phonon  effMtt 
depend  airongly  on  the  enecD'  diairibntion  of  eadirid  dec- 
trona. 

Let  na  firai  oonaider  a  aimplified  picmre  in  which  there  ia 
only  one  energy  aubband  and  SO  phonona  are  ne^ected;  that 
ia,  only  LO  and  acouadc  phonona  are  preaetf  far  the  QWL 
Thia  ajwiptMMMt  picture  ailowa  na  to  refine  d»  pme  ID  effect 
of  the  broadening  of  the  electron  energy  diatiibntion  on  the 
liniUhip  of  hot  pl«««twia  »n«t,  on  nnnUt^  dy. 

"— !<•«  Rfuie  3  ifinatratea  electron  coding  dynandca  in  a 
1S0X2S0  A*  QWI  at  r-30  K  fdlowing  iiddal  electron  ex¬ 
citation  at  an  energy  4.5  timea  the  LO  phonon  energy  for  two 
differem  Oanaaian  eiecnon  diatribndon  half  widdia:  30  and  4 
meV.  For  compariaon,  we  plot  the  eleciron  rrlaxaiiondynam- 
ica  widMut  nonequilibrira  optical  phonorm.  When  hot 


phonona  are  negiecied  the  cooling  dymica  diapi^  two 
dtniiqiiialiahle  atapea:  the  faat  ataga  (wMi  ■Oniiwwiiwwl  «hi. 

ladon)  due  to  the  caacade  emiaaion  of  optical  phonorm,  and 
tfaeaecood.alowaiageofdectronthenuiBimiondneioin- 
reracdona  with  acouadc  phonona.  It  aanat  be  nored.  bowover; 
that  in  the  given  time  acale  of  10  pa,  phonon  acat- 

leiing  doea  not  viaibly  infinence  ite  ify- 

namica  in  dua  QWI  with  a  rather  large  croaa  aecdon  of  ISO 
X250  k*.  Aa  one  can  aee  from  Fig.  5  the  very  fandal 
relaxation  aiage  (t<0J  pa)  ia  faaier  in  the  r"irTr  of  hot 
phermna.  The  higher  imnegmlilinwin  phnann  pripiltim,. 

created  (4  me V),  the  frater  ia  the  very  inidal  reiaxadon  atage. 
Thia  effect  can  be  mideritood  if  one  firatconaiden  the  aem- 

pemme  depwiHwife  eS  die  lelMaritw  »«at>  At  tfmpffa- 

tnrea  both  the  emiaaion  and  abaotpdon  niea  are  highei:  TUa 
leada  to  faat  eneigy  rediaiiifandon  of  excited  eiectrana.  The 
cooling  ate  of  electroru  which  *■»>■»  gprirai  piwnf  m. 
creaaea  and  that  of  electrona  which  riiaorb  phonona  decreaaea 
becanae  of  the  a'*'*  eneoy  dependence  of  ID  denaify  of 
atatea  and  acattering  latea.  The  increaae.  however;  ia  frater 
than  the  decreaae  due  to  the  aam*  Theniaie, 

the  total  ID  electron  gaa  cooling  me  >««•»»*— when  the 
*W«iw  «n»»gy  —HinaMt**  iW 

of  eprt***!  At  vety  fafgh  reaapnMma  tUa  h^pona 

on  a  very  abort  dare  acak  vriiile  eieciran  aaerihndon  at  low 
temperamme  adn  reaaaina  unchanged.  Hence,  the  vary  haWai 
electron  oooliDg  rate  fat  QWb  incieaaea  when  hacreaaing  the 
laidoe  teaoperamre,  provided  that  electrom  are  exched  weO 
above  optical  phonon  encagy  and  ihwmal  en- 

eagy.  dint  fat  bulk  materiala,  edme  the  eaaiaaion  md 
abeorpdon  ratea  increare  with  tBagf  the  reiaxadon  tree  ia 
hi|)ieratlowlaitioeiempetainiea.Thereaiionldbei»ieaa- 
peramre  dqiendence  of  the  initial  reiaxadon  rare  in  2D  aya- 
tema.)  To  obaerve  an  appreciable  ttanperomre  effect  on  the 
reiaxadon  rate  it  janeceaaaiydiatphonmorrupadon  number 
be  greater  than  1.  Under  phonon  eqnilibrinmauch  occupation 
mmabert  could  even  be  unachievable  in  a  aoUd  atato.  How- 
e^  due  to  alrong  buildup  of  noneqoiHbrinm  pi««K»««  at 

to  certain  pborton  modea  may  be  oonaidendily  higher  titan  i. 
Thia  ia  why  the  initial  rrlaxatioo  ia  frater  to  higher  nonegni- 

Htirhti  tfanSp  Iv  BIDOWer 

electron  energy  tfianibutiona  (Rg.  5). 

One  can  aee  from  Fig.  S  that  the  onaet  of  hot  pbanana 
leada  to  a  rnhatantial  redncdon  of  the  elecaon  gaa 
rate  to  rX).S  pa  doe  to  atrong  reabaotpdon  of  noneqnOib- 
rinm  pbormna.  The  onaet  of  hot  phonona  Occam  aooner  if  the 
electron  energy  diairibution  ia  narrower  (4  meV).  Hence, 
electron  cooling  ia  riower  to  narrow  eleciron  diatiibuHona. 
The  effea  of  narrowing  of  the  eleciron  eneigy  dteibndoo  ia 
aimilar  to  that  of  fatcreaaing  the  eleciron  oonceanadon  and, 
aa  we  have  already  diacnaaed,  ilia  a  pnrrfy  ID  effect.  Aahaa 
been  denronatraied  In  a  previooa  anbaecdon,  far  QWb  with 
tmallcroaBaectiona(40X40A^theaoonadcphononacal(er- 
ing  rate  b  hitfier  and  thia  acattering  b  mnch  more  indaadc 
than  far  QWb  widi  large  croai  aecdona  (150X250  A^. 
Therefore,  in  a  QWI  vrith  a  40X40  A^  croaa  aecdon,  aoooadc 
phonon  acattering  b  a  very  effective  eneigy  tfiaaipation 
mechaniwn  and  it  b  reaponaible  to  fiat  reiaxadon  and 
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pkoaoa.  ElaeM  cooetaMiM  a  aqari  UP  Ok*'.  oAar  f««am  mA 
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mearing  oat  liw  effiect  of  die  initial  broadeniag  of  the  dec- 
tran  ^stabatioa  oo  the  cooling  fate. 

We  have  alao  cooiideied  the  realiitic  caae  where  the 
maltiiiibbe»l  atmctaie  d  (be  QWI  ia  taken  kMo  aocoont 
along  with  all  poadbk  optical  phonon  modea  QjO  and  SO). 
Hgore  6  ahowa  the  deciron  coding  dynatnica  in  dda  leafiatk 
atractme  for  4  and  30  meV  electron  excitation  linewidtha. 
ThedependenceofhotpbooonbaildDpontiieeiecnandia- 
trihotion  broadening  ia  waahed  out  altnoat  compieaely  fat  tfaia 
reaiiatic  anucttiie  due  to  varioua  fancaanbband  and  faneianb- 
band  (TBnaiiiooa  aaaiated  by  the  LO  and  the  two  SO  modea. 
The  leaaoo  ia  that  nooeqnUibriam  phonon  peafca  in  9  apace 
in  tlda  caae  overtap  and  fonn  a  con^texbro^diatributiM  in 
9  apace  viitnally  independent  of  initial  election  diatribatioa 
The  main  effect  which  cornea  into  pUy  within  tiria  tealiatic 
model  ia  the  dqiendenoe  of  the  number  the  upper  aub- 
benda  iiivoived  in  election  cooUng  on  the  initial  eleciiQo 
cnenpr  tfiatributioa  F^nre  7  demonatratea  tito  time  evolution 
of  the  anbband  filling  by  elecnona.  In  die  caae  of  a  broad 
dectron  itiitid  datribution  (30  meV)  tine  are  more  anb- 
banda  occupied  by  dectrona  acaiiered  from  the  highenecgy 
tafl.  Therefore,  the  retnin  of  dectrona  to  the  ffaat  anbbnnd  ia 
alower  than  for  a  nairow  dectron  diatiibution  (4  meV).  Hot 
phonona  lead  to  anonger  intenabband  electron  redianibatioo 
and  alower  remm  to  the  loweat  anbband.  By  comparing  Figa. 
6and7ooenoticeathattfae<fiflerentoccuputionof  aubbanda 
for  4  and  30  meV  excitation  Hnewidlha  viiinally  doea  not 
affect  the  mean-dectron  em^  wUdi  coincide  for  both 
excitation  regiinea  after  3  pa  foUowing  excitation.  lUa,  at 
the  bat  glaaoe,  atnu^  behavior  ia  related  to  the  fact  that  the 
dectron  kinetic  cnagy  rdaied  to  ID  fine  motion  fat  each 
anbband  ia  higher  fbr  4  meV  excitarton  Ifaiewidtfa  due  to 
atronger  hot  phonon  effecta.  Ihia  dflerenoe  in  kfaMdc  ener- 
giea  la  co«Pivv»«»*»v*  by  the  higher  occupetion  of  the  upper 
Wftiip«nii«  in  the  caae  of  a  30  me  V  hnevrtdtli.  Conaeqnentiy, 
the  4  meV  curve  in  Rg.  6  containa  a  larger  part  of  kiiietic 


nay.  TlaeevolnfceotntwlanwewmadBedneatOoewi)  tofc. 
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energy  and  a  amalkr  pait  of  “potentid**  eneigy  than  the  30 
meV  curve. 

IV.  SUMMARY 

Sinmlation  leveala  comidex  depraideiice  of  hot>electnm 
gaa  cooUng  tiynamica  on  excfaation  energy,  lattice  Mn^era- 
tnie,  and  atmctnre  paiameieta  of  QWL  Efecnon  relaxatioo  fat 
sobpicotooood  If  oooiipUsd  by  tatsficiioo 

yiffc  qpicil  wbevets  tboBiiimkMiof  die 

is  by  SMCfliitDy  fa**^*— ^ 

f IfCtfCTt-irOTitir  g— 

two  ia  much  foater  in  a  thitmer  QWI  due  to  hi^Kr  aconatic- 
phonon  acanmng  laie  and  attonger  farlaaticity  of  dectron- 
acooatic  phonon  faueraction.  The  relaxation  timeaobtajned  fat 
our  aimidatioiia  are  in  good  agreement  with  the  leauita  of 
optied  meaamementa.*’*'**^ 

Vhriatioo  of  the  initid  dectron  energy  aubataiUiaHy 
dtangea  the  entire  picture  of  hot-eiectran  relaxation  due  to 
the  iaemetion  of  eiectrona  with  varioua  phonona  fat  (2WIa. 
Caicniationa  demonattate  potentid  poaaibHitiea  of  two  ef> 
fecta:  electrao  gaa  overcooliag  and  dynamic  fanmubbanri 
population  faweiaion.  Both  effectaexhihitatreng  dependence 
on  the  lattice  temperature. 

Population  inwraionia  well  pmwxinned  at  tow  tempera- 
tarea,  whereaa  dectron  gaa  overcooliiig  henelita  from  fal^ 
temperaturea.  At  tow  temperamrea  eiectrona  can  be 
“ngped”  in  the  upper  aubbainda  bdow  d»  optied  phonon 
emgy  and  atay  thm  for  tjuite  a  tong  time  defined  by  faner- 
mbband  electron  ecouatic  phonon  itneraction. 

Vfe  have  found  that  hot  phonon  efiecta  fat  (JWb  are  well 
prononnoed  for  electron  conoenirationa  eipid  to  or  Itigher 
than  10*  cm~'  and  depend  atrongly  on  the  energy  diatribn* 
tion  of  excited  elecnona.  Hot  phonon  effecta  become  wedtor 
aa  die  broadening  of  die  exdM  electron  enetgy  tfiatiibnlion 
"VTWBBM  TUa  reault  ia  fat  complete  cooliaat  to  the  caae  of 
bulk  matetiala  and  quatnnm  wdia  where  the  energy  diatiibn* 
tion  of  excited  eiectrona  viitually  doea  not  affect  dm  boOdqi 
of  nooequilibiinm  phonona  and  die  reabsoiption  rate. 
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12.  RECOMMENDATIONS 

The  density  matrix  in  the  coordinate  representation  is  an  extraordinary 
means  of  studyiiig  quantum  transport  in  nanoscale  devices.  While  not  discussed  in 
this  docummt  this  approach  is  beii^  used  to  study  the  formation  of  barriers  at 
metal/semiconductor  interfiices.  Presently,  the  denaty  nutrix  in  the  coordituUe 
representation,  along  with  Green’s  fiinction  iq)proadies  are  the  only  quantum 
device  simulations  expUdtly  accounting  for  disapation.  The  density  matrix  is  also 
the  only  quantum  transport  procedure  presently  crq>able  of  handling  both  dectron 
and  hole  tran^rt.  Its  ability  to  directly  compute  current  from  the  nonequilibrium 
quantum  distribution  function  is  a  agnificant  advance  over  those  current 
algorithms  that  employ  the  Tsu-Esaki  formulation.  Further,  the  presence  of  menu- 
driven  algorithms  fur  both  workstations  and  PCs  suggests  that  the  approach  taken 
at  SRA  will  provide  a  significant  design  tool. 
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